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Corticosteroids  are  often  prescribed  for  chronic  conditions.  Diseases  such  as 

rheumatoid  arthritis  and  asthma  require  long-term  drug  therapy  and,  regardless  of  the 

disease  or  route  of  administration,  both  the  desired  and  undesired  effects  must  be 

weighed.  In  some  cases,  systemic  activity  is  necessary  and  side  effects  are  unavoidable. 

Pharmacokinetic/Pharmacodynamic  (PK/PD)  models  may  be  constructed  to  assess  and 

predict  the  degree  of  systemic  effects  of  corticosteroids.  Blood  lymphocytes  were  chosen 

as  surrogate  markers  for  systemic  activity  since  their  decline  is  directly  related  to  possible 

immunosuppression  in  chronic  corticosteroid  therapy.  Total  lymphocytes  and  selected 

subtypes  were  monitored  in  a  clinical  study.  Lymphocyte  counts  are  suppressed  by 

corticosteroids,  endogenous  and  exogenous,  while  endogenous  Cortisol  itself  is 

suppressed  by  the  exogenous  corticosteroid  administered.  All  of  these  interactions  were 

included  in  a  PK/PD  model.  Since  nonlinear  binding  of  Cortisol  to  plasma  proteins 

influences  both  its  elimination  and  distribution,  a  new  pharmacokinetic  model  based  on 


IX 


the  flux  of  unbound  Cortisol  was  developed  to  describe  these  phenomena  and 
incorporated  into  the  lymphocyte  module  of  the  PK/PD  model. 

The  PK/PD  analysis  of  total  lymphocytes  and  selected  subtypes  revealed  that  total 
lymphocytes  were  suppressed  by  10%  after  9  mg  oral  budesonide.  Suppression  of  CD4 
and  CD  19  was  greater  and  that  of  CD8  was  minimal.  Based  on  these  results,  there  is  no 
apparent  benefit  in  monitoring  subtypes  over  total  lymphocytes  solely  for  obtaining  an 
adequate  surrogate  marker.  PK/PD  simulations  of  total  lymphocytes  using  parameters 
obtained  from  clinical  data  showed  that  the  degree  of  systemic  effects  after  one  dose  per 
day  depends  on  administration  time.  Effects  on  lymphocytes  are  least  for  4  AM 
administration  and  greatest  for  4  PM.  If  Cortisol  is  monitored  as  the  surrogate  marker,  the 
opposite  trend  is  observed.  Lymphocyte  reduction  is  inversely  related  to  Cortisol 
suppression.  Since  Cortisol  is  suppressed  by  exogenous  corticosteroids,  there  is  less  total 
corticosteroid  present  and  no  net  effects  on  lymphocytes  are  predicted  for  bioavailable 
budesonide  doses  of  0.1  to  0.8  mg  or  less.  This  buffering  capacity  is  most  apparent  when 
the  potency  of  the  drug  on  Cortisol  is  three  or  more  times  greater  than  the  potency  on 
lymphocytes. 


CHAPTER  1 
COMBINING  PHARMACOKINETICS  AND  PHARMACODYNAMICS 

Introduction 

This  research  focused  on  various  aspects  of  the  theory,  development,  and 
application  of  the  pharmacokinetics  and  pharmacokinetic/pharmacodynamic  (PK/PD) 
modeling  of  corticosteroids.  Anti-inflammatory  corticosteroids  are  used  to  treat  a  variety 
of  chronic  conditions  including  rheumatoid  arthritis,  Crohn's  disease,  and  asthma.  Since 
chronic  diseases  require  long-term  drug  therapy,  the  degree  of  both  the  desired  effects 
and  the  undesired  side  effects  must  be  weighed.  Potential  side  effects  include  glaucoma, 
hypertension,  hyperglycemia,  Cortisol  suppression,  osteoporosis,  and  peptic  ulcers 
(Caldwell  and  Furst  1991;  Howland  1996;  Sorkness  1998;  Lipworth  1999;  Tripathi, 
Parapuram  et  al.  1999).  The  most  important  therapeutic  actions  are  the  anti-inflammatory, 
immunosuppressive,  and  lympholitic  effects  of  these  drugs. 

The  use  of  inhaled  corticosteroids  for  asthma  allows,  in  effect,  local  delivery  of 
the  drug  to  the  lung.  The  goal  of  inhaled  corticosteroids  is  to  produce  long-lasting 
therapeutic  effects  at  the  pulmonary  target  site  and  to  minimize  systemic  effects  by  rapid 
clearance  of  the  absorbed  drug  (Derendorf,  Hochhaus  et  al.  1998).  However,  despite  this 
local  delivery,  the  entire  inhaled  dose  reaching  the  lung  (as  well  as  the  fraction  of  the 
swallowed  portion  of  the  inhaled  dose  not  eliminated  during  first-pass  metabolism)  will 
ultimately  enter  systemic  circulation.  Similarly,  delayed  release  oral  formulations  may 
allow  for  local  delivery  of  the  drug  for  the  treatment  of  inflammatory  bowel  disease. 


In  the  treatment  of  rheumatoid  arthritis  with  oral  corticosteroids,  systemic  activity 
is  necessary  and  undesired  side  effects  are  often  unavoidable.  This  is  due  to  the  fact  that, 
unlike  many  classes  of  drugs,  both  desired  and  undesired  effects  of  corticosteroids  are 
mediated  by  the  same  drug-receptor  interaction  (Hochhaus,  Mollmann  et  al.  1997). 

Mechanisms  of  Action 

Corticosteroids  contribute  to  a  broad  array  of  physiological  and  biochemical 
processes  such  as  carbohydrate,  protein,  and  lipid  metabolism  as  well  as  the  electrolyte 
and  water  balance  (Kragballe  1989;  Jusko  1994).  Receptor/gene-mediated,  receptor/non- 
gene-mediated,  and  nonspecific  effects  are  attributed  to  corticosteroids.  The  main 
mechanism  of  action  of  corticosteroids  is  generally  accepted  to  be  via  interaction  with 
specific  receptors,  which  results  in  an  increase  or  decrease  in  gene  transcription.  These 
are  characterized  by  slow  onset  and  slow  dissipation  of  the  effects.  Examples  of  such 
receptor/gene-mediated  effects  include  anti-inflammatory  activity  and  the  induction  of 
the  hepatic  enzyme  tyrosine  aminotransferase  (Xu,  Sun  et  al.  1995). 

The  glucocorticoid  receptor  is  found  mainly  as  a  water-soluble  protein  within  the 
cell  rather  than  as  an  integral  part  of  the  cell  membrane.  However,  membrane  binding 
sites  have  been  found  in  several  tissues  with  a  chemical  similarity  to  the  cytosolic 
corticosteroid  receptor.  Both  appear  to  cooperate  for  the  modulation  of  cellular  events. 
The  majority  of  mammalian  cells  have  been  found  to  contain  corticosteroid  receptors,  in 
agreement  with  the  importance  of  corticosteroids  for  the  regulation  of  energy  metabolism 
and  immune  and  anti-inflammatory  response  (Hochhaus,  Mollmann  et  al.  1997). 

In  addition  to  the  classical  receptor/gene-mediated  mechanism,  a  nucleus- 
independent  cytosolic  receptor-mediated  effect  has  also  been  proposed.  These  are 


characterized  by  an  immediate  onset  after  drug  administration  and  include  steroid- 
induced  Cortisol  suppression  and  the  redistribution  of  lymphocytes  (Balow,  Hurley  et  al. 
1975;  Fauci  1975;  Matsuse,  Shimoda  et  al.  1999;  Meibohm,  Derendorf  et  al.  1999).  The 
obstruction  of  calcium  influx  and  the  interference  with  c-AMP-dependent  protein  kinase 
are  two  examples  of  the  modes  of  action  for  non-gene-mediated  corticosteroid  effects. 

Clinical  efficacy  of  high  doses  of  corticosteroids  may  not  always  be  explained  by 
receptor  interactions  because  these  doses  result  in  plasma  and  tissue  concentrations  far 
exceeding  those  for  receptor  saturation.  High  doses  are  indicated  for  spinal  cord  injury, 
arthritis,  and  for  the  inhibition  of  rejection  during  organ  transplant  (van  den  Brink,  van 
Wijk  et  al.  1994;  Ferraris,  Tambutti  et  al.  2000).  High  doses  of  flunisolide  have  been 
proposed  for  treatment  of  acute  asthma  attacks  (Rodrigo  and  Rodrigo  1998).  Despite 
these  common  clinical  applications,  not  much  is  known  about  the  mode  of  action  for 
these  nonspecific  effects. 

Since  effects  are  mediated  via  a  single  receptor  type,  corticosteroids  with  different 
receptor  binding  affinities  do  not  differ  in  their  qualitative  or  quantitative  effects. 
However,  the  concentrations  necessary  to  achieve  the  desired  effects  will  be  different.  In 
general,  lower  receptor  affinity  may  be  overcome  by  increasing  the  dose,  an  indication 
that  high  pharmacodynamic  activity  is  not  necessarily  the  most  critical  factor  in 
corticosteroid  therapy.  It  has  been  suggested  that  the  pharmacokinetic  properties  of  the 
drug  may  be  the  most  important  considerations  for  optimizing  the  therapeutic  index 
(Hochhaus,  Mollmann  et  al.  1997). 

Pharmacokinetics 

The  treatment  of  asthma  has  been  significantly  improved  as  a  result  of  the 


introduction  of  corticosteroids  such  as  fluticasone  propionate  (FP)  and  budesonide  (BUD) 
having  improved  therapeutic  ratios  compared  to  other  drugs  such  as  hydrocortisone  (HC), 
predisolone  (P),  and  dexamethasone  (DEX).  Although  there  is  a  marked  increase  in  the 
relative  receptor  binding  affinity  (RRA),  this  improvement  is  mainly  due  to  optimized 
pharmacokinetic  properties  (Derendorf,  Hochhaus  et  al.  1998).  Pharmacokinetic  and 
pharmacodynamic  parameters  of  some  selected  corticosteroids  are  contained  in  the  table 
below.  The  corticosteroids  BUD  and  FP  are  characterized  by  a  high  clearance  (CL)  and  a 
low  oral  bioavailability  (Forai).  These  are  especially  useful  for  treating  asthma  since  the 
swallowed  portion  of  the  dose  undergoes  significant  first-pass  metabolism  and  the  drug 
entering  systemic  circulation  by  oral  or  pulmonary  absorption  is  rapidly  cleared. 
However,  modern  corticosteroids  differ  in  volume  of  distribution  (Vd)  and,  therefore, 
half-life  (ti/2).  Fluticasone  propionate  has  one  of  the  largest  steady-state  volumes  of 
distribution. 

Table.  Pharmacokinetic  and  pharmacodynamic  properties  of  selected  corticosteroids. 

Drug        RRA       fu       CL  (L/h)       VdsslL)      FnraJ  (%) 

HC  9       20  18 

P  16       25  6 

DEX  100      23  17 

FLU  180      20  58 

TA  233      29  37 

BUD         935      12  84 

FP  1800      10  69 

RRA=Relative  receptor  affinity  (compared  to  dexamethasone,  defined  as  100), 
fu=fraction  unbound  in  plasma,  CL=total  body  clearance,  Vdss=steady-state  volume  of 
distribution,  Forai=bioavailability  after  oral  administration. 
HC=hydrocortisone,  Prednisolone,  DEX=dexamethasone,  FLU=flunisoIide, 
TA=triamcinolone  acetonide,  BUD=budesonide,  FP=fluticasone  propionate. 
Parameter  values  were  obtained  from  previous  publications  (Rohdewald,  Rehder  et  al. 
1987;  Derendorf,  Mollmann  et  al.  1991;  Rohatagi,  Hochhaus  et  al.  1995;  Derendorf  1997; 
Derendorf,  Mollmann  et  al.  1997;  Hochhaus,  Mollmann  et  al.  1997;  Rohatagi,  Barth  et  al.' 
1997;  Derendorf,  Hochhaus  et  al.  1998;  Meibohm,  Hochhaus  et  al.  1999). 
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The  degrees  of  plasma  protein  and  tissue  binding  are  also  important 
considerations.  Binding  is  important  from  a  pharmacokinetic  perspective  (since  CL  may 
be  influenced  by  plasma  protein  binding  and  Vdss  is  determined  by  both  plasma  protein 
and  tissue  binding)  and  from  a  pharmacodynamic  perspective  (since  only  the  unbound 
fraction  of  the  drug  is  pharmacologically  active).  Most  corticosteroids  have  a  high  deg 
of  binding  to  plasma  proteins  with  a  small  unbound  fraction  in  plasma  (fu).  For  high 
extraction  drugs,  fu  is  a  factor  determining  oral  bioavailablity.  In  general,  overall  tissue 
binding  is  difficult  to  determine.  However,  tissue  binding  or  the  fraction  unbound  in 
tissues  (fut)  may  be  estimated  indirectly  using  the  steady  state  volume  of  distribution 
(Vdss)  since  this  term  is  dependent  on  the  ratio  of  fu:fut.  Due  to  the  importance  of  binding 
properties,  these  will  be  discussed  in  depth  and  correlated  to  other  physicochemical 
properties  such  as  solubility  and  lipophilicity.  Discrepancies  between  previously 
published  values  of  fu  (flunisolide)  and  nonlinear  plasma  protein  binding  (hydrocortisone) 
will  also  be  addressed. 

Pharmacokinetic/Pharmacodynamic  Models 
Corticosteroid  dosing  for  different  therapeutic  outcomes  is  usually  determined 
empirically  and  is  subject  to  considerable  debate.  The  goal  in  designing  a  dosing  regimen 
for  corticosteroids  (or  any  drug,  for  that  matter)  is  to  provide  a  dose  sufficient  to  achieve 
the  desired  therapeutic  effect  while  keeping  side  effects  to  a  minimum.  Even  for 
local/topical  administration,  undesired  systemic  (and  local)  side  effects  are  possible  if  not 
completely  likely.  Although  equally  applicable  to  local  effects,  the  following  discussion 
is  geared  mainly  towards  systemic  activity,  whether  desired  or  adverse. 

In  order  to  both  assess  and  predict  the  degree  of  systemic  effects  of 


corticosteroids,  mathematical  relationships  may  be  constructed.  These  relationships, 
usually  referred  to  as  models,  are  often  complex  since  they  must  include  both 
pharmacokinetics,  the  time  course  of  the  drug  in  the  body,  and  pharmacodynamics,  the 
effect  observed  at  certain  drug  concentrations.  More  formally,  pharmacokinetics 
describes  the  fate  of  drugs  in  the  body  over  a  period  of  time,  including  their  absorption, 
distribution,  localization  in  tissues,  biotransformation,  and  elimination  or  excretion.  A 
more  complete  definition  of  pharmacodynamics  is  the  study  of  the  biochemical  and 
physiological  effects  of  drugs  and  the  mechanisms  of  their  actions.  By  combining  the 
two,  pharmacokinetic/pharmacodynamic  (PK/PD)  models  enable  the  effects  of  a  drug  to 
be  described  as  a  function  of  time.  In  addition  to  describing  clinical  data,  PK/PD  models 
may  be  used  to  predict  the  effects  of  corticosteroids  having  different  kinetic  and  dynamic 
properties  and  to  identify  important  considerations  for  future  studies,  thereby  enabling  the 
design  of  more  informative  and  cost-effective  clinical  trials.  Other  goals  in  PK/PD 
modeling  include  testing  competing  hypotheses  for  differences  observed  in  the  efficacy 
of  drugs,  predicting  responses  under  new  conditions  such  as  various  doses  and 
administration  times,  and  estimating  inaccessible  model  parameters  such  as  the 
amplification  factor  in  a  biochemical  cascade  yielding  the  observed  effect  after  receptor 
activation.  It  should  be  noted  that  most  therapeutic  effects  of  corticosteroids  are  receptor- 
mediated  and  there  seems  to  be  no  substance-specific  difference  in  the  post-receptor 
reaction  cascade  (Derendorf,  Mollmann  et  al.  1997).  Therefore,  the  extent  and  duration  of 
the  effects  and  side  effects  depend  only  on  receptor  occupancy  determined  by  the 
availability  of  the  drug  at  the  receptor  site  and  its  affinity  to  the  receptor  (Hochhaus, 
Mollmann  et  al.  1997),  making  corticosteroids  ideal  candidates  for  PK/PD  modeling. 


As  stated  by  Dirac  (Dirac  and  Dalitz  1995),  the  physical  laws  necessary  for  the 
mathematical  theory  of  a  large  part  of  physics  and  the  whole  of  chemistry  are  completely 
known  and  the  difficulty  is  only  that  the  exact  application  of  these  laws  leads  to  equations 
too  complex  to  be  solved.  This  is  especially  true  when  attempting  to  construct 
mathematical  models  to  describe  biological  systems  or,  in  the  present  case,  the  fate  and 
action  of  drugs  in  the  body.  Due  to  the  complexity  of  some  kinetic  and  dynamic 
phenomena,  some  proposed  models  are  entirely  theoretical  without  close  correlation  to 
physiological/biochemical  events  occurring  in  vivo.  Examples  of  this  include  describing 
the  complex  distribution  and  elimination  of  corticosteroids  using  compartmental  models 
consisting  of  poly-exponential  equations  and  describing  their  effects  in  terms  of  cosine 
functions,  harmonic  oscillations,  and  Fourier  series  (Chakraborty,  Krzyzanski  et  al. 
1999).  The  present  study  was  not  approached  with  this  philosophy  and  every  attempt  was 
made  to  apply  models  based  on  the  physiological  elimination,  distribution,  and  mode  of 
action  as  well  as  the  physicochemical  properties  of  the  drug. 

Physiolofiical  Models 

The  goal  of  the  present  study  was  not  to  construct  mathematical  models  that 
simply  describe  data  sets  but  rather  to  adapt  or  develop  physiologically  based  models  that 
are  also  sound  mathematically.  By  doing  so,  it  is  not  only  possible  to  describe  kinetic  and 
dynamic  data  but  also  to  obtain  a  better  understanding  and  appreciation  of  the  complex 
events  occurring  in  vivo,  thereby  bridging  the  gap  between  mathematical  description  and 
physiological  reality.  The  PK/PD  models  used  in  this  study  may  best  be  described  as 
semi-empirical  in  that  some  aspects  of  the  models  are  based  on  physiological 
observations  while  others  are  based  on  the  straightforward  application  of  established 


physical  principles  to  biological  systems.  Examples  of  such  physiological  PK/PD  models 
for  corticosteroids  include  those  for  the  circadian  rhythm  of  Cortisol  (Rohatagi,  Hochhaus 
et  al.  1995),  the  reduction  of  blood  lymphocytes  (Meibohm,  Derendorf  et  al.  1999),  and 
pulmonary  targeting  (Hochhaus,  Mollmann  et  al.  1997)  or  deposition  (Gonda  1981). 

According  to  Gibaldi  (Gibaldi  and  Perrier  1982),  the  kinds  of  information 
required  by  the  model  can  be  classified  as  anatomical  (organ  and  tissue  volumes), 
physiological  (blood  flow  rates  and  enzyme  reaction  parameters),  thermodynamic  (drug- 
protein  binding  isotherms),  and  transport  (membrane  permeability).  The  PK/PD  models 
in  the  present  study  allow  the  inclusion  of  parameters  such  as  plasma  protein  binding, 
tissue  partitioning,  and  receptor  binding  affinities  obtained  in  vitro.  In  addition, 
absorption  and  distribution  may  be  described  by  diffusion  processes,  including  membrane 
permeability,  and  elimination  may  be  described  in  terms  of  physiological  models  such  as 
the  Wilkinson-Shand  expression  for  hepatic  clearance,  including  enzyme  activity 
(Wilkinson  and  Shand  1974). 

Surrogate  Markers 

Clinical  trials  are  the  standard  scientific  method  for  evaluating  a  new  drug, 
device,  or  procedure  for  the  treatment  or  prevention  of  disease.  The  phase  three  clinical 
trial  is  designed  to  evaluate  the  clinical  benefit  and  possible  side  effects  of  a  new  drug 
which  often  requires  many  participants  to  be  followed  for  a  long  period  of  time.  The 
primary  endpoint  should  be  clinically  relevant  to  the  patient  and  examples  include  death, 
loss  of  vision,  symptomatic  events  of  acquired  immunodeficiency  syndrome  (AIDS),  the 
need  for  ventilatory  support,  and  other  events  causing  a  reduction  in  quality  of  life 
(Fleming,  Prentice  et  al.  1994).  There  has  recently  been  great  interest  in  using  surrogate 


endpoints  such  as  tumor  shrinkage,  changes  in  cholesterol  levels,  blood  pressure,  white 
blood  cell  counts,  or  other  laboratory  measures  to  reduce  the  cost  and  duration  of  clinical 
trials  (Fleming  1994;  Fleming,  Prentice  et  al.  1994).  Surrogate  markers  are  objective 
observations  that  reliably  serve  as  early  indicators  for  clinically  meaningful  endpoints 
and  there  should  be  a  demonstrable  linkage  between  the  surrogate  and  the  ultimate 
endpoint.  Surrogate  markers  are  only  useful  when  meaningful  correlations  with  clinical 
events  are  well-established.  However,  it  is  a  misconception  that  if  a  measurable 
observation  correlates  to  the  desired  clinical  outcome,  then  it  may  be  used  as  a  valid 
surrogate  endpoint.  Proper  justification  for  such  a  replacement  requires  that  the  effect  of 
the  drug  on  the  surrogate  end  point  may  be  used  to  predict  the  clinical  outcome.  Criteria 
have  been  developed  to  validate  surrogate  endpoints  in  clinical  trials  (Boissel,  Collet  et 
al.  1992).  These  criteria  essentially  require  that  the  surrogate  must  be  a  correlate  of  the 
true  clinical  outcome  and  fully  capture  the  net  effect  of  the  drug  on  the  clinical  outcome. 
Thus,  the  greatest  potential  for  validating  a  surrogate  is  that  the  surrogate  is  in  the  only 
contributory  pathway  of  the  disease  process  and  that  the  entire  effect  on  the  true  clinical 
outcome  is  mediated  through  its  effect  on  the  surrogate. 

Although  this  strict  definition  leads  to  the  identification  of  the  best  possible 
surrogate  marker,  in  practice  the  mechanisms  of  drugs  are  often  complex  or  poorly 
understood  and  the  validity  of  a  surrogate  endpoint  may  rarely  be  rigorously  established 
using  these  criteria.  Thus,  it  may  be  argued  that  such  a  strict  definition  of  surrogates 
limits  the  information  to  be  gained  in  a  clinical  study  and  that  surrogate  markers  should 
be  thought  of  in  broader  terms.  Surrogates  can  not  only  predict  efficacy  but  can  also  be 
used  to  characterize  toxicity  or  undesired  effects  and  these  may  not  necessarily  correlate 
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to  the  same  pathophysiologic  process  as  the  effect  on  the  true  clinical  outcome.  Such 
surrogates  may  be  termed  biochemical  markers  or  "biomarkers"  (Colburn  1995;  Lee, 
Hulse  et  al.  1995;  Colburn  1997;  Colburn  2000).  Identifying  biomarkers  enhances  the 
ability  to  find  more  effective  and  less  toxic  therapies  more  accurately  and  efficiently. 
Surrogates  for  Systemic  Effects  of  Corticosteroids 

Surrogate  markers  are  particularly  important  in  the  study  of  corticosteroids  since 
many  corticosteroid  effects,  desired  or  otherwise,  are  either  difficult  to  quantify  or 
indeterminate  during  short-term  observation.  Two  such  markers  have  been  identified  fi 
the  systemic  effects  of  corticosteroids:  endogenous  Cortisol  (Rohatagi,  Hochhaus  et  al. 
1995)  and  blood  lymphocytes  (Meibohm,  Derendorf  et  al.  1999^  These  are  easily 
measured  in  blood  samples. 
Cortisol  Suppression 

Cortisol  is  very  sensitive  to  the  presence  of  exogenous  corticosteroids  and  the 
release  of  serum  Cortisol  is  often  dramatically  suppressed  after  drug  administration. 

Cortisol  suppression  has  been  used  as  an  indication  of  the  amount  of  corticosteroid 

\ 

entering  systemic  circulation  after  inhalation  and  has  been  suggested  as  a  surrogate 
marker  for  the  effects  on  the  hypothalamic-pituitary-adrenal  axis  (Meibohm,  Hochhaus  et 
al.  1997;  Meibohm,  Hochhaus  et  al.  1999).  However,  Cortisol  is  not  necessarily  a 
preferred  marker  for  systemic  activity.  As  stated  previously,  there  should  be  a  close 
correlation  between  a  surrogate  marker  and  the  relevant  effects  of  the  drug.  Although  the 
lack  of  Cortisol  can  alter  the  body's  regulation  of  protein,  carbohydrate,  and  fat 
metabolism,  blood  sugar  levels,  and  inflammatory  responses  as  observed  in  Addison's 
disease,  when  Cortisol  suppression  is  induced  by  the  administration  of  an  exogenous 
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corticosteroid,  the  presence  of  the  exogenous  corticosteroid  may  compensate  for  the  lack 
of  endogenous  Cortisol  in  systemic  circulation.  In  this  case,  Cortisol  is  merely  replaced  by 
another  corticosteroid,  albeit  a  more  potent  one  in  most  circumstances. 
Lymphocyte  Reduction 

In  this  study,  the  decline  in  blood  lymphocytes  was  chosen  as  the  surrogate 
marker  for  the  systemic  activity  of  corticosteroids.  This  is  an  appropriate  choice  since  the 
decline  in  lymphocytes  is  directly  related  to  the  suppression  of  the  immune  system  often 
observed  in  chronic  corticosteroid  therapy,  whether  it  be  an  undesired  side  effect  (Cupps, 
Edgar  et  al.  1984;  Pytsky  1991;  Wilckens  and  De  Rijk  1997)  or  the  desired  clinical 
outcome  (Wofsy  1990;  Dupont,  Huygen  et  al.  1985).  Unless  stated  otherwise,  the 
discussion  of  lymphocytes  and  lymphocyte  subsets  below  is  a  compilation  of  information 
obtained  from  general  immunology  texts  and  review  articles  (Hildemann  1984;  Paul 
1989;  Eisen  1990;  Kuby  1997;  Peakman  and  Vergani  1997). 

The  two  key  features  of  lymphocytes  are  specificity  and  memory.  Lymphocytes 
are  involved  in  the  specific  immune  defense  system  and  can  distinguish  foreign  cells  and 
antigens,  having  mechanisms  for  selecting  a  precisely  defined  target.  The  second  feature 
of  lymphocytes  is  memory,  allowing  immunization  and  resistance  to  reinfection  with  the 
same  microorganism.  Lymphocytes  constitute  between  25%  and  40%  of  all  white  blood 
cells.  Less  than  1%  of  all  lymphocytes  circulate  in  the  blood,  most  being  found  in  tissues. 
Lymphocyte  Subsets 

The  effects  of  corticosteroids  on  lymphocyte  subsets  were  also  investigated.  Since 
different  subtypes  play  different  roles  in  regulating  the  body's  immune  response,  there  is 
some  interest  in  the  effects  of  corticosteroids  on  subtypes,  particularly  helper  T  cells 
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(CD4)  and  suppressor  T  cells  (CD8)  and  their  ratio  (Gerblich,  Urda  et  al.  1985; 
Tornatore,  Venuto  et  al.  1998;  Gergely  1999).  It  is  generally  accepted  that  proper 
regulation  of  the  immune  system  requires  CD4  and  CD8  cells  be  present  in  a  2:1  ratio. 

Lymphocytes  are  divided  into  two  main  classes,  bone-marrow-derived  (B)  and 
thymus-derived  (T)  lymphocytes.  Natural  killer  (NK)  cells,  a  small  population  of 
lymphocytes,  do  not  express  characteristics  of  either  T  or  B  cells.  Although 
morphologically  similar,  T  and  B  cells  are  distinguishable  by  their  different  cell-surface 
components  and  are  classified  on  the  basis  of  these  surface  markers.  The  NK  cells  may 
also  be  distinguished  by  these  markers. 

The  cluster  designation  (CD)  is  based  on  surface  antigens  (designated  by 
numbers)  that  can  be  recognized  by  specific  antibodies.  Mature  T  cells  have  CD3/CD4 
(helper  cells)  or  CD3/CD8  (suppressor  or  cytotoxic  cells),  B  cells  have  CD19-CD22,  and 
NK  cells  have  CD1 6  and/or  CD56.  Of  the  mature  lymphocytes  in  circulation,  20-30%  are 
B  cells,  65-75%  are  T  cells,  and  less  than  5%  are  NK  or  null  cells. 

The  B  lymphocytes  mature  in  the  marrow  and  migrate  directly  to  the  spleen  and 
the  lymph  nodes.  They  are  then  ready  to  react  to  foreign  substances  (antigens).  The 
principal  function  of  B  cells  is  to  synthesize  and  secrete  antibodies  for  a  rapid  response  to 
extracellular  and  mucosal  microorganisms,  including  viruses  (if  part  of  their  life  cycle  is 
spent  in  the  extracellular  fluid).  When  these  cells  are  stimulated  by  encountering  a 
foreign  antigen,  they  go  through  several  cycles  of  cell  division  and  then  differentiate  into 
specialized  antibody-secreting  plasma  cells.  Each  B  cell  is  programmed  to  make  one 
specific  antibody.  The  secreted  antibodies  typically  interact  with  circulating  antigens 
such  as  bacteria  and  toxic  molecules. 
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The  T  lymphocytes  migrate  from  the  bone  marrow  to  the  thymus  where  they 
mature  before  migrating  to  the  peripheral  lymphoid  organs  (spleen  and  lymph  nodes). 
Several  subsets  of  T  cells  exist  and  they  serve  two  roles.  The  first  is  as  the  coordinator  of 
other  immune  responses,  accomplished  by  their  production  of  a  wide  variety  of  cytokines 
and  surface  cell  signals  (helper  and  suppressor  T  cells).  The  second  is  as  a  responder  to 
long-term  intracellular  infections  and  cancerous  cells  (cytotoxic  T  cells).  The  T 
lymphocytes  do  not  secrete  a  single  major  protein  such  as  immunoglobulin  but  carry  out 
their  function  primarily  by  direct  contact  with  other  cells. 

Helper  T  cells  (CD4)  recognize  foreign  peptides  on  cell  surfaces  and  secrete 
factors  that  stimulate  other  cells  involved  in  the  immune  response.  For  example,  CD4 
cells  stimulate  to  B  cells  to  bring  about  proliferation  and  the  secretion  of  antibodies.  The 
CD4  cells  also  activate  cytotoxic  T  cells  and  natural  killer  cells.  The  helper  T  cells  are 
central  players  in  most  immune  responses  and  are  the  targets  of  human 
immunodeficiency  virus  (HIV),  the  retrovirus  that  causes  acquired  human 
immonodeficiency  (AIDS).  The  CD4  subset  is  a  common  surrogate  marker  for  the 
progression  of  this  disease  since  the  loss  of  helper  T  cells  causes  failure  of  the  immune 
system. 

Suppressor  T  cells  (CD8)  can  switch  off  antibody  production  by  suppressing  the 
action  of  other  cells  in  the  immune  system,  notably  B  cells  and  helper  T  cells.  Suppressor 
cells  are  important  for  terminating  a  response  to  a  foreign  pathogen  but  may  also  prevent 
inappropriate  immune  responses  to  self-antigens  (autoimmunity). 

Cytotoxic  T  lymphocytes  also  have  the  CD8  surface  marker  and  have  the  simplest 
mode  of  action  of  all  T  cells.  They  recognize  and  kill  cells  that  have  foreign  peptides  on 
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their  surfaces,  for  example,  fragments  of  viral  proteins  on  the  surface  of  a  virus-infected 
cell.  Tumor  cells  are  also  recognized,  having  proteins  that  do  not  normally  constitute  the 
cell.  Cytotoxic  T  cells  kill  infected  or  tumor  cells  by  a  lysis  mechanism  or  by  initiating 
programmed  cell  death  (apoptosis)  in  the  target  cell  and  are  responsible  for  the  rejection 
of  tissue  and  organ  grafts. 

Null  cells  or  natural  killer  (NK)  cells  are  large  granular  lymphocytes  derived  from 
stem  cells  in  bone  marrow  and  are  able  to  kill  cells  through  a  mechanism  known  as 
antibody  dependent  cellular  cytotoxicity,  similar  to  cytotoxic  T  lymphocytes.  The  NK 
cells  are  spontaneously  cytotoxic  to  a  variety  of  targets,  including  certain  cancer  cells. 
Combined  Effects  of  Endogenous  and  Exogenous  Corticosteroids 

As  will  be  discussed,  the  degree  of  systemic  activity  as  indicated  by  a  decrease  in 
lymphocyte  counts  depends  on  the  pharmacokinetic  disposition  of  the  drug,  the  potency 
of  the  drug,  and  the  interaction  of  the  drug  with  endogenous  compounds,  namely  Cortisol. 
Diurnal  variations  were  previously  reported  for  lymphocytes  and  subsets  (Signore, 
Cugini  et  al.  1985)  and  later  correlated  to  the  circadian  rhythm  of  Cortisol  (Milad,  Ludwig 
et  al.  1994;  Meibohm,  Derendorf  et  al.  1999).  Lymphocyte  counts  are  influenced  by  both 
endogenous  Cortisol  and  the  exogenous  administered  drug.  In  addition,  exogenous 
corticosteroids  suppress  the  release  of  endogenous  Cortisol.  Thus,  to  adequately  describe 
the  effects  of  corticosteroids  on  lymphocytes  with  a  PK/PD  model,  both  Cortisol 
concentrations  and  Cortisol  suppression  by  an  exogenous  corticosteroid  must  be 
considered.  A  schematic  diagram  of  this  model  is  shown  in  Figure  1-1. 

The  ultimate  goal  of  this  research  was  to  adapt  and  apply  a  previously  developed 
physiological  model  for  describing  the  combined  effects  of  endogenous  and  exogenous 
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A  PK/PD  model  for  the  combined  effects  of  corticosteroids 
on  blood  lymphocytes 


ECSs  suppress  the  release 
Ik.  of  endogenous  Cortisol 


Exogenous  Corticosteroids  (ECSs) 


Direct  Effect:  ECSs  inhibit 
the  influx  of  lymphocytes 
into  the  blood  from  the 
extravascular  space 


Lymphocytes  exhibit  a 
circadian  rhythm  inversely 
proportional  to  that  of  Cortisol 

Indirect  Effect:  Cortisol 
suppression  by  ECSs  results 
in  less  Cortisol  present  to 
influence  lymphocyte  influx 


Figure  1-1 .  A  schematic  diagram  of  the  PK/PD  model  for  describing  the  combined 

effects  of  Cortisol  and  an  exogenous  corticosteroid  on  blood  lymphocytes. 
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corticosteroids  on  blood  lymphocytes  (Meibohm,  Derendorf  et  al.  1999)  and  subtypes.  As 
will  be  discussed,  the  degree  of  immunosuppression  after  the  administration  of 
corticosteroids  may  be  assessed  (and  later  predicted  using  the  PK/PD  model)  by 
monitoring  blood  lymphocytes,  total  and  subtypes,  satisfying  many  of  the  criteria  for  the 
validation  of  surrogate  markers  stated  previously. 

Budesonide 

Most  of  the  research  presented  in  this  dissertation  involved  budesonide  (BUD),  a 
highly  potent  corticosteroid  indicated  for  the  treatment  of  lung  diseases  (Newhouse, 
Knight  et  al.  2000)  and  other  inflammatory  conditions  such  as  Chrone's  syndrome 
(Friend  1998)  and  psoriasis  (de  Jong,  Ferrier  et  al.  1995). 

As  shown  in  Figure  1-2,  budesonide  contains  an  asymmetric  16a,  17a-acetal 
group  and  exists  as  a  1 :1  mixture  of  the  two  epimers,  22R  and  22S,  that  have  the  same 
qualitative  pharmacologic  effect.  The  22R  epimer  is  two  to  three  times  more  potent  than 
22S  (Ryrfeldt,  Edsbacker  et  al.  1984). 


CH9OH 


Figure  1-2.  Chemical  structure  of  budesonide. 
The  clearance  of  BUD  has  been  reported  to  be  83.7  L/h,  a  value  close  to  liver 
blood  indicating  that  BUD  is  a  high  extraction  drug,  with  a  half  life  of  2.8  h,  a  steady 


17 

state  volume  of  distribution  of  183  L,  and  an  oral  bioavailability  of  10.7%  (Ryrfeldt, 
Andersson  et  al.  1982;  Thorsson,  Edsbacker  et  al.  1994).  The  plasma  protein  binding  of 
BUD  is  88%  without  affinity  to  transcortin  (Ryrfeldt,  Andersson  et  al.  1982).  The  major 
metabolites  are  6p-hydroxy  budesonide  and  16a-hydroxy  prednisolone  and  there  appears 
to  be  no  extra-hepatic  biotransformation  to  these  metabolites  (Edsbacker,  Jonsson  et  al. 
1983).  However,  it  has  been  hypothesized  that  BUD  undergoes  reversible  conjugation 
with  fatty  acids  in  vivo  to  yield  palmitates  and  oleates  via  an  ester  linkage  (Tunek,  Sjodin 
et  al.  1997;  Miller-Larsson,  Mattsson  et  al.  1998).  The  conjugates  could  then  serve  as  a 
reservoir,  providing  sustained  release  of  BUD.  If  these  conjugates  are  present  in 
sufficient  quantities,  BUD  may  exhibit  improved  lung  retention  after  inhalation  and  an 
improved  pharmacological  profile  with  prolonged  anti-inflammatory  effects  in  the  lung 
and  a  greater  local-to-systemic  effect  ratio  (Wieslander,  Delander  et  al.  1998). 

Additional  Therapeutic  Corticosteroids 
Other  corticosteroids  such  as  triamcinolone  acetonide  (TA),  fluticasone 
propionate  (FP),  and  flunisolide  (FLU)  were  investigated  to  a  lesser  degree.  However,  a 
considerable  amount  of  effort  was  invested  in  describing  the  pharmacokinetics  of 
hydrocortisone  (HC)  or  endogenous  Cortisol.  The  reason  for  this  is  that  hydrocortisone 
exhibits  nonlinear  binding  to  plasma  proteins  and  it  was  believed  that  this  could  effect 
both  the  clearance  and  distribution  of  this  drug.  Various  compartmental  models  were 
investigated  and  a  new  model  was  proposed  based  on  the  distribution  of  free  drug 
between  compartments.  This  model  was  designed  to  describe  both  nonlinear  clearance 
and  nonlinear  distribution  (a  new  concept)  attributed  to  the  concentration-dependent  free 
fraction  of  hydrocortisone  in  systemic  circulation. 


18 


Hypotheses 
The  following  hypotheses  were  tested  using  either  experimental  methods, 
compartmental  analysis  of  clinical  kinetic  data,  PK/PD  modeling  procedures  for  clinical 
kinetic  and  dynamic  data,  or  PK/PD  simulations: 

1 .  Corticosteroid  binding  to  plasma  proteins  and  tissues  may  be  determined  using  in 
vitro  experimental  procedures.  These  binding  phenomena  follow  the  same  trends 
established  for  solubility  and  lipophilicity. 

2.  The  nonlinear  binding  of  hydrocortisone  to  plasma  proteins  is  responsible  for  both 
nonlinear  elimination  and  nonlinear  distribution.  For  large  doses,  no  singular  CL  or 
Vd  term  is  sufficient  for  adequately  describing  the  concentration-time  profile.  The 
pharmacokinetics  of  hydrocortisone  may  be  described  for  any  given  dose  by  using  a 
free  concentration  gradient-based  distribution  model  in  conjunction  with  the 
Wilkinson-Shand  expression  for  hepatic  clearance,  which  is  a  function  of  the 
unbound  fraction  in  plasma.  This  model  may  be  integrated  numerically  and  used  to 
describe  and  predict  nonlinear  pharmacokinetics  with  various  input  parameters. 

3.  The  reduction/redistribution  of  blood  lymphocytes  may  be  used  as  a  surrogate  marker 
for  the  systemic  activity  of  corticosteroids.  The  pronounced  first-pass  metabolism  of 
budesonide  (BUD)  precludes  the  realization  of  substantial  systemic  effects  for  normal 
oral  doses.  However,  differential  effects  may  be  observed  for  lymphocyte  subtypes. 
Systemic  effects  may  be  quantified  by  comparing  the  area  under  the  effect-time 
profiles  after  placebo  to  that  after  budesonide  administration. 

4.  Administration  time,  potency,  pharmacokinetic  disposition,  the  circadian  rhythm  of 
Cortisol,  and  Cortisol  suppression  affect  the  degree  of  systemic  activity  (as  measured 
by  the  reduction  in  blood  lymphocytes)  observed  after  corticosteroid  administration. 
Minimal  systemic  effects  are  observed  when  the  administered  dose  entering  systemic 
circulation  is  counterbalanced  by  the  steroid-induced  Cortisol  suppression. 
Theoretically,  an  exogenous  corticosteroid  may  be  administered  at  a  dose  such  that  no 
net  systemic  effects  are  observed. 


CHAPTER  2 
BINDING  TO  TISSUES  AND  PLASMA  PROTEINS 

Background 

Physiological  pharmacokinetic  models  have  been  used  to  quantitatively  describe 
the  distribution  and  elimination  of  drugs  (Chen  and  Gross  1979).  These  models  differ 
from  conventional  compartmental  models  in  that  many  (or  all)  of  the  parameters  have  an 
anatomical,  physiological,  or  physicochemical  basis.  One  of  the  most  important 
parameters  in  these  models  is  the  partition  coefficient  (Kp).  Partition  coefficients  are 
sometimes  referred  to  as  distribution  coefficients  and  usually  defined  as  the  ratio  of  drug 
concentration  in  tissue  to  that  in  plasma  or  blood,  which  is  used  to  estimate  the 
disposition  of  the  drug  (Poulin  and  Theil  2000).  The  value  of  Kp  is  directly  proportional 
to  the  volume  of  distribution  at  steady-state  and  it  has  been  suggested  that  Kp  influences 
clearance  in  addition  to  distribution  (Lam,  Chen  et  al.  1982).  The  proposed  relationship 
between  Kp  and  clearance  was  not  considered  in  the  present  study  and  may  be  subject  to 
debate. 

In  addition  to  describing  drug  distribution,  the  degree  of  tissue  binding  is 
important  because,  for  a  given  dose,  it  determines  the  amount  of  free  drug  in  the  tissue 
available  for  pharmacological  activity  as  only  unbound  drug  may  interact  with  receptors. 
One  of  the  best  examples  of  this  is  the  local  activity  in  the  lung  after  corticosteroid 
inhalation.  It  has  been  suggested  that  lung  tissue  partitioning  is  a  critical  pharmacokinetic 
parameter  and  the  term  "pulmonary  volume  of  distribution"  has  been  coined,  relating  the 
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amount  of  drug  deposited  in  the  lung  to  the  pharmacologically  active  free  drug 
concentration  within  lung  tissue  immediately  after  inhalation  (Hochhaus,  Mollmann  et  al. 
1997).  The  pulmonary  volume  of  distribution  (Vdpul)  is  a  necessary  parameter  in  the 
pharmacokinetic/pharmacodynamic  (PK/PD)  model  used  to  assess  the  pulmonary 
selectivity  of  inhaled  corticosteroids  under  ongoing  development  in  our  research  group. 
This  model  is  shown  in  Figure  2-1 .  The  free  concentration  of  drug  in  the  lung,  Cfflung), 
after  inhalation  is 

CAlung)  =  -^-  = £m> 

Vd*      Vlung.Kp{lung)  ^2A 

where  D,nh  is  the  inhaled  dose  deposited  in  the  lung,  V,ung  is  the  volume  of  the  lung,  and 
Kp  is  the  partition  coefficient  of  the  drug  into  lung  tissue.  Human  lung  volume  is  1.17  L 
(Davies  and  Morris  1993). 

There  is  no  general  consensus  on  the  most  appropriate  method  for  the 
determination  of  Kp  values;  Both  in  vivo  and  in  vitro  procedures  have  been  used  and  good 
correlations  were  observed  between  the  two  methods  (Kato,  Hirate  et  al.  1987).  Although 
algorithms  have  been  proposed  for  assessing  Kp  values  after  various  routes  of 
administration  (Chen  and  Gross  1979),  the  most  straightforward  experimental  design 
involves  the  constant-rate  infusion  of  the  drug  until  steady-state  is  reached,  the  time  at 
which  free  levels  in  the  plasma  and  tissue  are  equal.  In  vivo  partition  studies  have  been 
performed  for  loteprednol  etabonate  (LE)  in  rats  (Hochhaus,  Chen  et  al.  1992)  and  for 
prednisolone  (P)  in  laboratory  rabbits  and  rats  (  Khalafallah  and  Jusko  1984a  and  1984b; 
Mishina  and  Jusko  1994a  and  1994b).  F^^nheLE  study  in  rats,  significantly  higher  levels 
^yer£foundinliyer  and  kidney  compared  to  other  tissues. 
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Figure  2-1.  Pharmacokinetic/Pharmacodynamic  model  for  assessing  pulmonary 
selectivity  of  inhaled  corticosteroids. 
kdiSS-the  dissolution  rate  of  drug  deposited  in  the  lung  after  inhalation,  Vdpu|-the 
pulmonary  volume  of  distribution,  ka>pui-the  absorption  rate  of  drug  from  the  lung  into 
sustemic  circulation,  ka,Gi-the  absorption  rate  of  the  swallowed  portion  of  the  dose  from 
the  GI  tract  (gastrointestinal  tract)  into  systemic  circulation,  CL-total  body  clearance  of 
the  drug,  Vd-the  volume  of  distribution,  ke-the  first  order  elimination  rate  constant. 

Roughly  20%  of  an  inhaled  dose  is  deposited  in  the  lung  and  the  remaining 
portion  is  swallowed  and  subject  to  first-pass  metabolism  prior  to  entering  systemic 
circulation.  The  drug  deposited  in  the  lung  dissolves  and  the  resulting  free  concentration 
in  the  lung  (initially  treated  as  an  isolated  one-compartment  system)  is  determined  by  the 
pulmonary  volume  of  distribution,  dependent  on  the  degree  of  binding  to  tissue.  Free 
drug  may  interact  with  corticosteroid  receptors  in  the  lung  or  be  absorbed  into  systemic 
circulation.  In  animal  studies,  pulmonary  selectivity  has  been  defined  as  the  difference  in 
receptor  occupancy  between  lung  and  systemic  tissues  (Hochhaus,  Mollmann  et  al. 
1997). 
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For  the  steady-state  infusion  experiment  in  rabbits,  the  total  P  concentration  in  lung  tissue 
was  2.9  times  greater  than  the  unbound  concentration  in  plasma  (Khalafallah  and  Jusko 
1984a). 

Few  in  vivo  studies  have  been  reported  for  human  tissue  partitioning  of 
corticosteroids  and  these  involved  the  inhalation  of  a  drug  prior  to  lung  resection  surgery. 
After  budesonide  inhalation,  the  concentration  in  lung  tissue  was  8  times  greater  than  that 
inj)lasma  (total  rather  than  unboundconcentrations)  (Van  den  Bosch,  Westermann  et  al. 
1993).  After  the  inhalation  of  fluticasone  propionate,  central  lung  concentrationsjyere  3 
t0  Ujmes  greater  than  periphsraLiung  concentratix^fe-v^wey^a^irn^^^in  times 
greater  than  serum  levels  (Esmailpour,  Hogger  et  al.  1997). 

In  vitro  studies  have  been  performed  using  both  human  and  rat  tissues  and  it  was 
found  that  partitioning  of  aromatic  compounds  into  rat  adipose  tissue  is  predictive  of 
partitioning  into  human  adipose  tissue  (Pierce,  Dills  et  al.  1996).  In  order  to  reach  steady- 
state  conditions  during  in  vitro  experiments,  sufficient  incubation  times  must  be  allowed 
for  complete  distribution  of  drug  into  the  tissues  and  equilibrium  between  the  tissue  and 
the  surrounding  drug  solution.  These  partitioning  experiments  are  commonly  performed 
by  placing  a  tissue  sample  in  a  buffer  solution  (pH  7.4)  containing  the  drug  and 
incubating  the  tissue/solution  for  18  to  24  h  at  37°C  (Jepson,  Hoover  et  al.  1994). 

Related  in  vitro  experiments  have  been  described  for  human  lung  tissue  (Hogger 
and  Rohdewald  1998).  Flunisolide  (FLU),  beclomethasone  dipropionate  (BDP), 
beclomethasone  monopropionate  (BMP),  hydrocortisone  (HC),  BUD,  and  FP  were 
considered  in  this  experiment.  In  this  study,  human  lung  tissues  samples  were  placed  in 
corticosteroid  solutions  (0.3  ug/mL)  for  1  h.  The  tissues  were  then  placed  in  human 
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plasma  for  re-equilibration  of  drug  between  tissue  and  plasma.  The  results  were  reported 
qualitatively  and  the  uptake  of  drug  into  tissue  exhibited  the  following  rank  order: 

FP  >  BDP  >  BMP^>  BUD  >  FLU  >  HC. 
From  viewing  the  diagrams  in  this  publication,  the  uptake  of  FP  appears  to  be  roughly  2.5 
times  greater  than  that  of  BUD.  Tissue  concentrations  after  re-equilibrium  in  human 
plasma  exhibited  the  following  rank  order: 

BDP  =  BMP  <  FP  <  BUD  =  HC  <  FLU. 


A  possible  drawback  in  the  design  of  these  experiments  is  that  the  solutions  used  for  the 
initial  uptake  of  FP  and  BDP  into  lung  tissue  exceeded  the  aqueous  solubilities  of  these 
two  drugs.  It  is  not  possible  to  directly  compute  and  compare  Kp  values  from  these  results 
to  those  in  other  publications  and  those  in  the  present  study  because  total  rather  than 
unbound  plasma  concentrations  were  reported,  tissue  concentrations  were  reported  in 
units  of  ug  drug/g  of  tissue  (preventing  Kp  from  being  expressed  as  a  unitless  parameter), 
and  the  experiments  were  conducted  at  ambient  room  temperature  rather  than  at  37°C. 
In  addition  to  predicting  the  distribution  and  free  fractions  of  drugs  in  tissues, 
partition  coefficients  are  of  key  interest  in  the  design  of  new  therapeutic  agents.  It  has 
been  shown  that  the  partition  coefficient  of  drugs  in  a  two-phase  system  of  octanol  and 
water  (normally  expressed  as  the  log  of  the  partition  coefficient,  logP)  is  related  to 
biological  activity  (Dunn  and  Hansch  1974),  uptake  into  tissues  (Cornford  1982),  and 
solubility  (Valvani,  Yalkowsky  et  al.  1981).  Traditionally,  the  relationship  to  activity  has 
been  explained  in  terms  of  thermodynamics.  One  the  driving  forces  for  the  drug-receptor 
interaction  is  the  so-called  hydrophobic  interaction.  In  the  systems  involving  water  as  the 
polar  phase,  the  hydrophobic  effect  is  always  the  driving  force  that  governs  the 
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distribution  process  irrespective  of  the  interacting  or  noninteracting  nature  of  the 
substances  studied  (Ruelle  2000).  For  lipophilic  drugs  in  an  aqueous  environment, 
complexation  with  a  receptor  results  in  a  favorable  change  in  the  Gibb's  free  energy.  This 
favorable  change  (a  negative  change  since  energy  is  released)  is  often  due  to  an  increase 
in  entropy  rather  than  a  decrease  in  enthalpy. 

Estimations  of  Kp  values  of  new  drug  candidates  using  conventional  in  vitro 
and/or  in  vivo  methods  is  time  and  cost  intensive.  For  this  reason,  octanol-water  or  olive 
oil-water  is  often  used  to  calculate  partition  coefficients.  It  has  been  shown  that  the 
solubility  in  phospholipids  is  a  function  of  solubility  in  octanol-water  (Poulin  and 
Krishnan  1995)  and  tissue  partitioning  equations  have  been  developed  in  which  the  only 
compound-dependent  input  parameters  are  the  lipophilicity  parameters,  such  as  olive  oil- 
water  or  octanol-water  partitioning,  and/or  unbound  fraction  in  plasma,  fu,  determined  in 
vitro  (Poulin,  Schoenlein  et  al.  2001). 

Partition  coefficients  are  used  in  the  selection  of  compounds  for  further 
development.  Biological  activity  was  related  to  lipid-water  partition  coefficients  as  early 
as  1901  and  most  modern  quantitative  structure-activity  relationships  (QSAR)  are  based 
on  octanol-water  partition  coefficients  (Van  Valkenburg  and  American  Chemical  Society. 
Division  of  Pesticide  Chemistry.  1972;  Hansch  and  Leo  1979).  The  "Hansch  approach" 
to  understanding  drug  action  is  based  on  the  assumption  that  physicochemical  properties 
governing  drug  transport  and  drug-receptor  interaction  can  be  factored  into  electronic, 
hydrophobic,  and  steric  components.  One  of  the  most  widely  used  equations  for  relating 
biological  activity  (BA)  to  these  physicochemical  properties  within  a  series  of 
structurally  related  compounds  is 
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BA  ~  log 


FC 


=  -a ■  (logP)2  +b  \ogP  +  pcr  +  SEs  +c  Eq.  2-2 


where  EC50  is  the  drug  concentration  producing  50%  of  the  maximal  response,  P  is  the 
octanol-water  partition  coefficient  (the  hydrophobic  factor),  a  is  an  electronic  factor,  Es  is 
related  to  steric  interactions,  and  a,  b,  p,  5,  and  c  are  constants  (Van  Valkenburg  and 
American  Chemical  Society.  Division  of  Pesticide  Chemistry.  1972;  Dunn  and  Hansch 
1 974).  The  expression  above  is  written  as  a  quadratic  equation  in  terms  of  logP  to 
account  for  any  nonlinear  dependence  of  drug  activity  on  the  hydrophobic  character  of 
the  molecule.  This  is  not  always  the  case  and  a  simplified  equation, 


BA 


EC 


=  a-logP  +  b,  Eq.  2-3 


is  often  sufficient  for  QSARs  (Van  Valkenburg  and  American  Chemical  Society. 
Division  of  Pesticide  Chemistry.  1972;  Hansch  and  Leo  1979).  In  the  case  of  lipophilic 
molecules  (such  as  most  therapeutic  corticosteroids),  relative  lipid  content  (logP)  has 
been  shown  to  be  the  sole  mechanistic  determinant  of  tissue-plasma  partition  coefficients 
(Haddad,  Poulin  et  al.  2000).  To  aid  in  assessing  the  hydrophobic  factors  in  these 
expressions,  various  empirical  and  semi-empirical  algorithms  have  been  proposed  (Flynn 
1971;  Bodor  and  Huang  1992a;  Bodor  and  Huang  1992b;  Poulin  and  Theil  2000;  Ruelle 
2000). 

As  stated,  the  degree  of  plasma  protein  binding,  often  expressed  as  the  fraction 
unbound  (fu),  impacts  the  pharmacokinetic  disposition  of  drugs  by  modulating  the 
primary  pharmacokinetic  parameters,  clearance  and  volume  of  distribution.  In  addition,  fu 
is  important  from  a  pharmacodynamic  perspective  because  only  free  drug  is 
pharmacologically  active.  Thus,  fu  is  necessary  in  most  PK/PD  modeling  procedures. 
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Objectives 
The  steady  state  tissue  partitioning  of  selected  corticosteroids  was  initially 
assessed  in  vitro  using  lung  tissue  (rat)  to  obtain  an  estimate  of  the  pulmonary  volume  of 
distribution,  an  important  consideration  in  the  pharmacokinetic/pharmacodynamic 
(PK/PD)  modeling  of  inhaled  corticosteroids.  Budesonide  (BUD)  and  triamcinolone 
acetonide  (TA)  were  considered  because  these  two  corticosteroids  were  used  in  previous 
animal  studies  in  our  group  and  the  data  obtained  in  the  present  study  may  be 
incorporated  into  the  pharmacokinetic/pharmacodynamic  (PK/PD)  model  used  in  the 
previous  studies. 

Values  of  logP  were  calculated  using  a  semi-empirical  quantum  mechanical 
program  and  compared  to  the  experimentally  determined  tissue  and  plasma  binding.  This 
was  done  in  order  to  determine  the  correlation  between  binding/partitioning  phenomena 
and  other  physicochemical  properties  (stemming  from  molecular  size,  shape,  structure, 
and  functional  groups)  such  as  aqueous  solubility. 

Plasma  protein  binding  was  determined  for  various  corticosteroids  in  rat  and 
human  plasma  in  order  to  allow  total  drug  concentrations  to  be  converted  into 
pharmacologically  active  free  levels  in  the  PK/PD  model. 

In  previous  animal  studies  and  applications  of  the  PK/PD  model  for  assessing 
pulmonary  selectivity,  receptor  binding  was  used  as  a  surrogate  marker  for  local  (lung) 
and  systemic  (various  tissues)  effects.  In  these  receptor  binding  studies,  receptor 
occupancies  were  reported  under  the  assumption  that  all  binding  represented  binding  to 
the  corticosteroid  receptor  and  that  non-specific  binding  was  minimal.  Thus,  non-specifi 
binding  in  lung  cytosol  was  assessed  to  validate  this  assumption. 
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Tissue  Partitioning 
Methods 

Sprague-Dawley  male  rats  were  sacrificed  by  decapitation.  Tissues  (lung,  liver, 
kidney,  spleen,  brain,  heart,  and  thigh  muscle)  were  immediately  removed  and  excess 
blood  was  removed  by  blotting  the  tissues  with  paper  towels.  Portions  of  these  tissues 
(roughly  0.5  g)  were  accurately  weighed  and  tissue  volumes  were  determined  for  these 
tissues  (n=8)  by  displacement  using  phosphate-buffered  saline  solution  (PBS:  24  g 
dibasic  sodium  phosphate,  4  g  monobasic  sodium  phosphate,  4  g  sodium  chloride  in  4  L 
distilled  water,  pH  7.4)  in  a  10  ml  graduated  cylinder.  The  degree  of  drug  partitioning 
was  determined  by  placing  the  weighed  tissue  samples  in  15  ml  Falcon  tubes  and  adding 
10  mL  corticosteroid  solution  in  PBS  of  known  concentration.  For  BUD,  the  number  of 
replications  was  n=14  and  for  TA,  n=6.  Tubes  were  placed  in  a  rotating  rack  in  an  oven 
preheated  to  37°C.  Due  to  the  fact  that  the  liver  is  a  metabolizing  organ,  partitioning  into 
liver  tissue  was  assessed  at  4°C  to  minimize  enzyme  activity.  Tissue  samples  were 
incubated  at  constant  temperature  with  rotation  for  18  h.  Tissues  were  removed  and 
homogenized  in  5  mL  distilled  water.  Portions  (1  mL)  of  the  tissue  homogenates, 
solutions  in  which  the  tissues  were  incubated,  and  corticosteroid  solutions  (incubated  for 
18  h  at  37°C  without  the  addition  of  any  tissue  samples)  were  extracted  with  two  2-mL 
volumes  of  ethyl  acetate.  Extractions  were  evaporated  under  nitrogen  gas  (Meyer  N-Evap 
analytical  evaporator,  model  no.  1 12,  Organomation  Associates,  Inc.)  and  reconstituted  in 
0.5  mL  acetonitrile/water  (50:50).  Concentrations  were  determined  using  reverse-phase 
HPLC  with  a  C-18  column  (Zorbax  SB-C18,  4.6  mm  x  7.5  cm,  PN866953.902)  and  a 
mobile  phase  of  acetonitrile/water/trifluroacetic  acid  (50/49.07/0.03).  Samples  (50  pL) 
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were  injected  via  an  auto  sampler  (Perkin-Elmer  auto  sampler  model  ISS-100).  The  flow 
rate  of  mobile  phase  was  lmL/min  (Milton  Roy  multiple  solvent  delivery  system 
CM4000).  Detection  was  set  at  254  nm  (Milton  Roy  programmable  wavelength  detector 
SM4000)  and  the  results  were  recorded  using  a  Hewlett  Packard  integrator  (model  HP 
3394A).  Calibration  standards  (1,  2,  5,  10,  25,  50  ug/mL)  were  prepared  from  a  stock 
solution  of  the  drug  in  methanol. 

Tissue  masses  were  converted  to  volumes  using  previously  determined  densities 
for  each  of  the  selected  tissues  except  muscle.  Muscle  density  was  not  determined 
experimentally  and  a  density  of  1 .0  g/mL  was  used  for  skeletal  muscle,  as  in  other  studies 
(Khalafallah  and  Jusko  1984b).  This  enabled  the  tissue  partition  coefficients,  Kp,  to  be 
determined  as  a  unitless  parameter.  The  Kp  values  were  calculated  directly  from 
experimentally  determined  tissue  and  incubation  solution  concentrations  and  defined  as 

where  CtiSSue  is  the  drug  concentration  in  the  tissue  and  CSO|Ution  is  the  final  concentration 
in  the  incubation  solution.  The  fraction  of  drug  bound  to  non- water  components  in  each 
tissue,  reported  as  a  percent  (%fbt),  was  calculated  using 

%A=ioo-^A  Eq;2.5 

where  Kp  is  the  tissue  partition  coefficient  and  fw  is  the  fractional  water  content  in  tissue, 
averaging  0.761  for  the  tissues  considered  in  this  study  (Khalafallah  and  Jusko  1984b). 
The  percent  of  total  drug  that  is  unbound  (%fut)  and  pharmacologically  active  drug  in  the 
tissue  is 

%fu,=  100 -ft,,.  Eq.2-6 
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For  comparison,  octanol-water  partition  coefficients  (expressed  as  logP)  were 
predicted  for  various  corticosteroids  using  a  semi-empirical  method  relating  logP  to  the 
structure,  size,  shape,  and  functional  groups  of  a  molecule  (Bodor  and  Huang  1992a).  In 
this  user-friendly  program,  the  chemical  structure  of  the  molecule  is  entered  using  a 
graphical  interface  and  logP  values  are  generated  automatically  after  optimization  of  the 
ground-state  structure  of  the  molecule,  based  on  semi-empirical  quantum  mechanics  and 
molecular  orbital  theory. 
Results  and  Discussion 

The  results  of  the  density  experiments  are  shown  in  Table  2-1.  Densities  were 
very  close  to  unity  for  most  tissues.  Lung  density  was  slightly  less  and  this  value  is  likely 
due  to  the  air  content  in  the  tissue  immediately  after  excision. 

Calculated  Kp  values  for  BUD  and  TA  are  shown  in  Table  2-2  and  the  calculated 
fb,  values  (expressed  as  a  percentage)  are  shown  in  Table  2-3.  Previously  published 
results  for  prednisolone  (P)  (Khalafallah  and  Jusko  1984b)  are  also  included  for 
comparison.  In  the  experiments  conducted  in  this  study,  smallest  Kp  values  were 
observed  for  muscle  tissue,  largest  values  were  observed  for  liver  (BUD)  and  kidney 
(TA).  Similar  results  have  been  reported  for  other  corticosteroids  in  vivo  with  a  lesser 
degree  of  partitioning  into  muscle  for  P  (Khalafallah  and  Jusko  1984b)  and  a  greater 
degree  of  partitioning  into  liver  and  kidney  for  loteprednol  etabonate  (Hochhaus,  Chen  et 
al.  1992).  The  larger  standard  deviation  for  liver  tissue  binding,  observed  in  the  present  as 
well  as  other  published  studies,  is  likely  due  to  the  loss  of  drug  through  enzymatic 
activity  (Khalafallah  and  Jusko  1984b).  No  distinct  trends  were  observed  for  partitioning 
of  BUD,  TA,  and  P  into  other  tissues. 
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Table  2-1 .  Densities  of  selected  tissues. 


Tissue 

Density  (e/mL) 

Lung 

0.72  ±0.12 

Liver 

1.08  ±0.06 

Kidney 

1.07  ±0.21 

Spleen 

0.94  ±0.17 

Brain 

0.98  ±0.15 

Heart 

1.07  ±0.16 

Muscle 

* 

*  Muscle  density  was  not  determined  experimentally. 


Table  2-2.  Partition  coefficients  (Kc)  for  budesonide  (BUD),  triamcinolone 
acetonide  (TA),  and  prednisolone  (P)  in  selected  tissues. 
Tissue  KP(BUD)  KpjTA)  Kp  (P) 

Lung  8.60  ±2.95       3.01  ±  1.08  2.86 

Liver  12.20  ±7.16      7.82  ±2.64      0.38-4.47 

Kidney         8.49  ±4.01       6.02  ±1.65  2.91 

Spleen         7.72  ±  3.30       4.19  ±1.45  1.16 

Brain  5.50  ±3.21       3.96  ±1.11 

Heart  9.26  ±3.06       4.71  ±  1.35  2.92 

Muscle        2.11  ±0.44  —  1.54 

Values  of  Kp  for  BUD  and  TA  were  calculated  according  to  Eq.  2-4  with  n=14  and  n=6, 
respecively.  Values  for  P  are  previously  reported  results  from  other  researchers 
(Khalafallah  and  Jusko  1984b). 


Table  2-3.  Tissue  binding  (%fbt)  for  budesonide  (BUD),  triamcinolone 
acetonide  (TA),  and  prednisolone  (P)  in  selected  tissues- 
Tissue         %fbL(BUD) 

Lung         /  91.2  ±  28.1 

Liver         I  93.8  ±55.0 

Kidney       XLO  ±  43. oj 

Spleen         90.1  ±38.6 

Brain         C^6.2±503^> 

Heart  9L8T310 

Muscle        63.9  ±  6.8 
Values  of  %fb,  were  calculated  using  Eq.  2-5  and  based  on  the  values  of  partition 
coefficients  (Kp,  Table  2-2)  determined  experimentally  and  previous  published  values  of 
the  fractional  water  content  in  each  tissue  (Khalafallah  and  Jusko  1984b). 


%fb,(TA) 

%fb,  (P) 

74.7  ±  26.8 

72.6 

90.3  ±30.5 

58.7-84.3 

87.4  ±24.0 

72.8 

81.8  ±28.3 

34.9 

80.8  ±22.6 

59.3 

83.8  ±24.0 

72.9 

... 

50.3 
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Interestingly,  the  partitioning  of  TA  into  lung  tissue  was  somewhat  less  than  the 
partitioning  into  other  tissues.  This  trend  was  not  observed  for  BUD  and  it  is  possible  that 
the  larger  Kp  for  BUD  in  lung  tissue  is  be  due  to  the  formation  of  fatty  acid  conjugates,  as 
was  hypothesized  in  other  publications  (Miller-Larsson,  Mattsson  et  al.  1998; 
Wieslander,  Delander  et  al.  1998).  However,  the  conjugates  were  not  assayed  in  the 
present  study  and  the  importance  of  these  conjugates  may  not  be  stated  conclusively.  In 
comparison  to  TA,  a  larger  uptake  of  BUD  was  also  observed  in  vitro  in  previous  human 
lung  tissue  studies  (Hogger  and  Rohdewald  1998). 

Predicted  octanol-water  logP  values  for  various  corticosteroids  are  shown  in 
Table  2-4.  Consistent  with  the  results  obtained  in  the  in  vitro  tissue  partitioning 
experiments,  larger  logP  values  were  observed  for  BUD.  With  the  exception  of  BMP,  the 
predicted  logP  values  followed  the  same  trend  as  drug  uptake  into  human  lung  tissue 
(Hogger  and  Rohdewald  1998). 

Table  2-4.  Predicted  logP  values  for  selected  corticosteroids. 

Corticosteroid        LogP 

TA  1.03 

BMP  2.42 

Fluticasone  2.63 

BUD  3.24 

FP  3.89 

Values  of  LogP  for  octanol-water  partitioning  were  not  determined  experimentally.  The 
reported  values  were  calculated  using  a  previously  developed  semi-empirical  quantum 
mechanical  method  (Bodor  and  Huang  1992a). 

Binding  to  Plasma  Proteins 
Methods 

Whole  blood  from  Sprague-Dawley  rats  was  collected  into  heparinized  tubes 
immediately  after  decapitation  or  cardiac  puncture.  Tubes  containing  whole  blood  were 
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centrifuged  at  3000  rpm  (DYNAC  II  centrifuge,  Clay  Adams/a  division  of  Becton, 
Dickson,  and  Company)  for  20  min  at  room  temperature  to  precipitate  blood  cells. 
Plasma  was  removed  with  a  pipette  and  placed  in  borosilicate  tubes  and  refrigerated  at 
4°C  while  preparations  were  made  to  continue  the  experiment  (not  more  than  two  hours). 

Corticosteroid  solutions  in  rat  plasma  were  prepared  using  a  stock  solution  of  the 
drug  in  methanol  to  yield  a  final  plasma  solution  of  drug  in  0.1%  methanol.  Dilutions 
were  also  performed  with  phosphate  buffered  saline  (PBS)  solution  to  provide  control 
solutions.  All  solutions  were  warmed  to  37°C  before  filtration.  Free  drug  was  separated 
from  the  protein-bound  drug  by  ultrafiltration  (Centrifree  micropartition  device  for  in 
vitro  diagnostic  use,  30  Kdal  cutoff,  Amicon  Bioseparations,  Millipore  Corporation). 
Three  1-mL  portions  of  the  plasma  solution  and  the  controls  were  pipetted  into  Centrifree 
micropartition  devices.  These  were  placed  into  a  centrifuge  (Beckman  Centrifuge  model 
JA-20)  that  had  been  pre-warmed  to  37°C  and  centrifuged  for  10  min  at  4000  rpm  (a 
force  of  1935  g).  Less  than  20%  of  the  total  plasma  volume  was  filtered  to  prevent 
dramatic  changes  in  free  concentrations.  An  optimal  filtration  time  of  10  min  was 
previously  determined  by  filtration  rate  experiments  as  shown  in  Figure  2-2.  Portions  of 
all  solutions  were  retained  unfiltered  in  order  to  assess  total  concentrations.  All  binding 
experiments  were  performed  in  triplicate. 

Calibration  standards,  plasma  solutions  (filtered  and  unfiltered),  and  PBS 
solutions  (filtered  and  unfiltered)  were  subjected  to  the  same  preparation  procedure 
before  analysis.  The  procedure  was  as  follows: 

Aliquots  (100  uL)  of  each  solution  to  be  analyzed  were  placed  in  heavy  walled 
borosilicate  culture  tubes.  Dexamethasone  (DEX,  20  ug/mL  in  methanol)  was  used  as  an 
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Plasma  filtered  vs  time 
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Figure  2-2.  Filtration  rate  of  human  plasma  through  Centrifree  partition  devices. 

internal  standard  (IS)  for  TA  and  FLU  samples  and  TA  (10  ng/mL  in  methanol)  was  used 
as  an  IS  in  BUD  samples.  The  IS  (100  uL)  and  800  uL  PBS  were  added  to  the  samples  to 
produce  a  total  volume  of  1  mL.  The  TA  and  FLU  samples  were  extracted  with  2  mL  of 
ethyl  acetate,  capping  and  shaking  for  15  min  in  an  Eberbach  shaker  set  on  low.  Tubes 
were  then  centrifuged  at  3000  rpm  for  10  min.  Ethyl  acetate  was  taken  off  the  top  of  each 
tube  with  disposable  glass  pipettes  and  placed  into  borosilicate  culture  tubes.  The 
extraction  procedure  was  performed  twice.  Samples  were  evaporated  under  nitrogen  gas 
until  dry  and  residues  were  reconstituted  in  100  uL  acetonitrile/water  (50:50). 

The  DEX  and  TA/FLU  contents  of  each  sample  were  determined  using  reverse- 
phase  HPLC  analysis.  The  calibration  curves  for  TA  and  FLU  were  constructed  using 
calibration  standards  in  methanol  (1,  4,  10,  20,  50  ug/mL).  The  mobile  phase  used  was 
30%  acetonitrile  and  0.03%  trifluoroacetic  acid  in  distilled  water.  The  HPLC  system  used 
for  analysis  was  the  same  as  that  described  in  the  previous  section  regarding  drug 
concentrations  in  tissue  extractions. 
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The  BUD  samples  were  analyzed  separately  by  a  liquid  chromatography-mass 
spectrometry  (LC-MS)  assay  following  solid  phase  extraction  (Lindberg,  Paulson  et  al. 
1987;  Kronkvist,  Gustavsson  et  al.  1998;  Krishnaswami,  Mollmann  et  al.  2000). 
Retention  of  the  corticosteroids  by  the  membrane  in  the  filtration  device  was  determined 
by  comparing  concentrations  of  filtered  and  unfiltered  solutions  of  the  drug  in  PBS 
solution. 

A  similar  method  for  determining  the  protein  binding  of  TA,  BUD,  and  FLU  in 
human  plasma  was  performed  using  human  plasma  purchased  from  Civitan  Regional 
Blood  Center  (Gainesville,  FL).  Experiments  were  performed  within  24  h  of  obtaining 
plasma  from  the  blood  center.  Plasma  was  kept  under  refrigeration  and  never  frozen. 

Nonspecific  binding  of  TA  and  BUD  in  rat  lung  cytosol  was  also  determined  by 
ultrafiltration.  The  cytosol  was  prepared  according  to  previously  published  methods 
(Hochhaus,  Rohdewald  et  al.  1983).  The  binding  of  TA  was  determined  in  undiluted 
cytosol  and  that  of  BUD  was  determined  in  diluted  cytosol,  cytosohwater  (1:3). 

The  tissue  and  plasma  protein  binding  data  may  be  combined  to  estimate  the 
volume  of  distribution  at  steady  state  (Vdss).  This  volume  term  is  proportional  to  the  ratio 
of  the  unbound  fraction  of  drug  in  plasma  to  that  in  tissue, 

r"ss     vp^yi     .  Eq.  2-7 

J  ul 

where  Vp  is  the  volume  of  plasma  (~3  L),  V,  is  the  volume  of  tissue  water  (-38  L),  fu  is 
the  unbound  fraction  in  plasma,  and  fut  is  the  unbound  fraction  in  tissue.  The  fraction  of 
drug  bound  to  plasma  proteins  (fb)  is  related  to  the  free  (unbound)  fraction  (fu)  by 

fu=X~fb  Eq.  2-8 

The  unbound  fraction  in  tissue  may  be  calculated  using  Eq.  2-4  through  Eq.  2-6. 
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Values  for  Vdss  calculated  using  the  fut  determined  in  the  present  study  would 
most  certainly  underestimate  the  actual  volumes  of  distribution  at  steady-state.  The 
reason  for  this  is  that,  in  general,  corticosteroids  exhibit  pronounced  partitioning  into 
fatty  tissues.  Since  adipose  tissue  partitioning  was  not  assessed,  absolute  values  of  Vdss 
may  not  be  calculated  accurately.  Instead,  the  ratio  of  tissue-to-plasma  protein  binding 
was  calculated  and  compared  to  previously  reported  Vdss  values.  Values  of  the  Vdss  (and 
tissue-to-plasma  protein  binding  ratios)  and  other  relevant  parameters  for  BUD  and  TA 
were  correlated.  The  parameter  values  for  BUD  were  normalized  to  those  for  TA  in  order 
to  facilitate  comparisons. 
Results  and  Discussion 

The  results  of  plasma  protein  binding  in  rat  and  human  plasma,  expressed  as  a 

percent  of  the  fraction  bound  (fb),  are  shown  in  Table  2-5.  Nonspecific  binding  in  cytosol 

is  given  in  Table  2-6.  There  was  no  significant  binding  to  the  filter  or  the  filtration  tube 

when  PBS  solutions  containing  the  corticosteroids  were  filtered. 

Table  2-5.  Plasma  protein  binding  (%  fh)  of  selected  corticosteroids. 
Species  TA  BUD  FLU 

Rat  87.2  ±5.2 

Human      68.8  ±1.8    89.2  ±2.3    71.5  ±3.9 

Table  2-6.  Nonspecific  binding  in  rat  lung  cvtosol  of  TA  and  BUD. 
Corticosteroid       Cvtosol  dilution  %  f j, 

TA  Undiluted  18.57  ±3  94 

BUD  1:3  16.12  +  7^56 

Plasma  protein  bindings  are  in  close  agreement  with  previously  published  values 
of  71%  for  TA  (Rohatagi,  Hochhaus  et  al.  1995)  and  88%  for  BUD  (Ryrfeldt,  Andersson 
et  al.  1982).  The  fb  of  71.5%  for  FLU  is  more  plausible  than  the  20%  bound  reported  in 
diluted  human  plasma  (Tomlinson,  Runkel  et  al.  1982).  Nonspecific  binding  of 


36 

corticosteroids  of  less  than  20%  in  cytosol  is  not  predicted  to  be  of  any  consequence  in 
the  interpretation  of  previously  acquired  receptor  binding  data.  The  reason  for  this  is  that 
at  least  80%  of  the  drug  exists  as  free  drug  inside  the  cell  and  is  able  to  interact  with  the 
receptor.  If  significant  binding  occurred,  the  receptor  occupancies  determined  in  previous 
studies  would  need  to  be  corrected  for  the  non-specific  binding.  This  is  not  the  case  and  it 
is  likely  that  any  perturbation  due  to  non-specific  binding  in  the  receptor  binding  assay  is 
within  typical  experimental  error. 

Summary  and  Conclusions 

For  BUD  and  TA,  there  is  a  very  close  correlation  between  the  tissue  partition 
coefficient,  solubility,  octanol-water  partitioning,  and  receptor  affinity.  These  and  other 
parameters  (normalized  to  the  values  for  TA)  are  shown  in  Figure  2-3.  This  suggests  that 
these  parameters  be  strongly  influenced  by  the  hydrophobic  character  (logP)  of  the 
molecules.  Thus,  although  electronic  and  steric  factors  can  not  be  completely  neglected, 
these  are  not  determining  factors  for  tissue  distribution  and  solubility.  Trends  for  receptor 
binding  and  solubility  have  already  been  correlated  to  lipophilicity  (Hogger  and 
Rohdewald  1998).  However,  these  have  been  established  on  a  qualitative  rather  than 
quantitative  basis.  From  viewing  the  values  tabulated  in  Figure  2-3,  there  is  a  roughly  a 
3-fold  difference  in  the  tissue  partitioning/binding,  solubility,  lipophilicity,  and  receptor 
affinity  values  between  BUD  and  TA.  The  larger  difference  in  the  RRA  ratio  suggests 
that  either  TA  has  some  steric  hindrance  in  interacting  with  the  receptor  or  that  BUD  has 
a  more  favorable  electronic  interaction,  yielding  a  more  negative  change  in  enthalpy  (and 
therefore  Gibb's  free  energy)  upon  receptor  binding. 

The  difference  between  fu  values  for  BUD  and  TA  is  not  as  great  as  that  between 
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Normalized  Parameters  for  TA  and  BUD 
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Figure  2-3.  Normalized  partitioning,  binding,  and  distribution  parameters  for 
triamcinolone  acetonide  (TA)  and  budesonide  (BUD). 
RRA-relative  receptor  affinity,  Kp-experimentally  determined  tissue  partition 
coefficient  (lung),  logP-estimated  octanol-water  partition  coefficient,  1 /Sol-inverse 
of  the  aqueous  solubility,  1/fut-inverse  of  the  experimentally  determined  unbound 
fraction  in  tissue  (lung),  1/fu-inverse  of  the  experimentally  determined  unbound 
fraction  in  plasma,  Vdss-the  volume  of  distribution  at  steady-state.  Parameter 
values  not  determined  in  this  study  were  obtained  from  tabulated  data  in  previous 
publications  (Rohatagi,  Hochhaus  et  al.  1995;  Hochhaus,  Mollmann  et  al.  1997; 
Derendorf,  Hochhaus  et  al.  1998;  Hogger  and  Rohdewald  1998). 


38 

logP  and,  therefore,  protein  binding  may  not  be  explained  solely  on  the  basis  of  the 
hydrophobic  interaction  or  lipophilicity.  There  is,  however,  a  very  close  correlation 
between  the  degree  of  plasma  protein  binding  and  the  experimental  tissue  binding-to- 
protein  binding  ratio  as  well  as  the  Vdss  values  obtained  from  clinical  data.  This  suggests 
that  despite  the  inability  to  calculate  absolute  values  of  Vdss  directly  from  the  data 
acquired  in  the  present  study,  the  in  vitro  tissue  experiments  provided  an  excellent 
estimate  of  the  relative  overall  tissue  partitioning  of  BUD  and  TA. 

It  should  be  noted  that  the  in  vitro  procedures  conducted  in  this  study  only  reflect 
the  permeability  and/or  affinity  of  the  drug  to  a  particular  tissue.  Perfusion  was  not  a 
factor.  This  is  of  no  consequence  when  these  parameters  are  used  in  the  PK/PD  model  for 
assessing  pulmonary  selectivity  because  the  binding  in  the  lung  immediately  after 
inhalation  does  not  depend  on  perfusion.  The  Kp  value  for  lung  tissue  may  be  used  to 
convert  the  amount  of  an  inhaled  dose  deposited  in  the  lung  to  the  respective  free 
concentration.  By  knowing  both  the  concentration  and  the  affinity  to  the  corticosteroid 
receptor,  initial  receptor  occupancy  may  be  calculated  for  a  given  dose.  The  decline  in 
receptor  occupancy  is  proportional  to  the  decline  in  drug  concentration  in  the  lung  as  the 
drug  is  absorbed  into  systemic  circulation.  As  stated,  there  is  a  3-fold  difference  in  both 
Kp  and  RRA  between  BUD  and  TA.  For  the  same  dose,  initial  free  concentrations  in  the 
lung  will  be  3  times  lower  for  BUD  compared  to  TA.  However,  the  RRA  for  BUD  is 
roughly  3  times  greater  than  that  for  TA.  Thus,  for  a  given  dose,  the  same  initial  receptor 
occupancy  will  be  observed  for  both  drugs  and  any  differences  in  pulmonary  selectivity 
between  BUD  and  TA  are  due  to  the  degree  of  lung  deposition,  retention  of  drug  in  the 
lung,  the  bioavailability  of  the  swallowed  fraction  of  the  dose,  and  the  half-life  of  the 
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drug  once  it  is  absorbed  into  systemic  circulation.  Sample  calculations  are  provided  in 
Figure  2-4  for  a  10  ug  dose  and  the  inhibition  of  T  cell  proliferation.  In  this  example,  the 
initial  effects  for  BUD  and  TA  were  almost  identical,  at  71%  and  70%,  despite  the 
difference  in  potency  (RRA).  The  dependence  of  pulmonary  selectivity  on 
pharmacokinetics  (namely  clearance)  has  been  addressed  in  previous  studies  (Hochhaus, 
Mollmann  et  al.  1997)  and  the  results  of  the  present  study  further  support  this  hypothesis. 


Pulmonary  Effects  of  BUD  and  TA 

Anti-Inflammatory  Response:  Inhibition  of  T  Cell  Proliferation 
Inhaled  Dose  (ug)  Cone  (ug/L)  Effect  (%) 
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The  pharmacologically  free  concentration  in  the  lung  is  determined  by  the 

pulmonary  volume  of  distribution  for  the  drug,  budesonide  (BUD)  or  triamcinolone 

acetonide  (TA).  Lung  volume  (human)  is  1.17  L;  Kp  is  the  partition  coefficient  of 

the  drug.  Effects  (anti-inflammatory  response)  were  calculated  using  an  Emax  model. 


Figure  2-4.  Pulmonary  effects  for  the  same  dose  of  budesonide  CBUD^)  and 

triamcinolone  acetonide  (TA)  on  the  inhibition  of  T  cell  proliferation. 
Calculations  were  performed  using  experimentally  determined  partition  coefficients  for 
BUD  and  TA  into  lung  tissue  of  8.6  and  3.0,  respectively,  from  Table  2-2,  a  lung  volume 
of  1 . 1 7  L  (Davies  and  Morris  1993),  and  previously  reported  EC50  values  for  the 
inhibition  of  T  cell  proliferation  (Hogger  and  Rohdewald  1998).  A  value  of  100%  was 
used  as  the  maximal  effect,  Emax. 


CHAPTER  3 
PHARMACOKINETICS  OF  HYDROCORTISONE 


Background  Information  and  Theory 
The  suppression  of  endogenous  Cortisol  has  been  used  as  a  biomarker  (discussed 
in  Chapter  1)  in  pharmacokinetic/pharmacodynamic  (PK7PD)  studies  of  exogenous 
corticosteroids  (Derendorf,  Mollmann  et  al.  1997).  Cortisol  is  also  an  important 
component  in  the  PK/PD  model  describing  the  reduction  of  blood  lymphocytes  observed 
after  the  administration  of  corticosteroids  because  the  number  of  lymphocytes  circulating 
in  the  blood  is  modulated  by  the  total  plasma  concentration  of  corticosteroid,  endogenous 
plus  exogenous  (Meibohm,  Derendorf  et  al.  1999).  Thus,  in  order  to  describe  the  effects 
of  exogenous  corticosteroids  on  lymphocytes  adequately  using  a  PK/PD  model,  Cortisol 
and  Cortisol  suppression  must  be  included  in  the  model.  In  addition  to  the  usual 
pharmacokinetic  properties,  the  Cortisol  module  of  the  PK/PD  model  must  describe  its 
circadian  rhythm.  This  has  been  accomplished  using  a  variable  production  (or  release) 
rate  of  Cortisol  (Rohatagi,  Hochhaus  et  al.  1995).  The  degree  of  Cortisol  binding  to 
plasma  proteins  must  also  be  included  because  not  only  does  this  determine  the  free 
fraction  of  Cortisol  that  is  pharmacologically  active  (on  lymphocytes),  as  will  be 
discussed,  it  may  also  influence  the  pharmacokinetic  disposition  of  Cortisol. 

Both  endogenous  Cortisol  and  exogenous  hydrocortisone  were  considered  this 
facet  of  the  present  PK/PD  study.  Determining  the  effects  of  nonlinear  binding  on 
pharmacokinetics  using  only  endogenous  Cortisol  data  (after  placebo)  is  problematic  due 
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to  the  variable  release  rate  of  Cortisol.  However,  by  suppressing  the  release  of 
endogenous  Cortisol  with  dexamethasone  prior  to  hydrocortisone  administration,  the 
pharmacokinetics  of  Cortisol/hydrocortisone  may  be  assessed  in  the  absence  of  any 
diurnal  variations  (Derendorf,  Mollmann  et  al.  1991).  Therefore,  pharmacokinetic  models 
for  hydrocortisone  will  be  discussed  prior  to  addressing  the  more  complex  expressions 
needed  to  describe  and  predict  endogenous  Cortisol  levels. 
Binding  to  Plasma  Proteins 

Hydrocortisone,  or  endogenous  Cortisol,  exhibits  saturable,  nonlinear  binding  to 
transcortin  and  linear  binding  to  albumin  (Ballard  1979a  and  1979b).  This  is  due  to  the 
high  affinity  of  hydrocortisone  to  transcortin,  a  plasma  protein  present  in  much  lower 
concentrations  than  albumin,  and  the  low  affinity  of  hydrocortisone  to  albumin.  Thus,  the 
pharmacologically  free  concentration  is  not  constant  but  concentration-dependent. 
Assuming  that  hydrocortisone  binds  only  to  albumin  and  transcortin  and  that  there  is  only 
one  binding  site  per  protein,  the  total  hydrocortisone  concentration  may  be  written  as 

CTo>ai  =   1+        —     +  KAlb  ■  QAlb  ■  Cf  +Cf  Eq.  3- 1 

where  Cf  is  the  concentration  of  free  (unbound)  hydrocortisone,  KTc  is  the  affinity 
constant  for  transcortin  (3.0  x  107  L/mol),  KAlb  is  the  affinity  constant  for  albumin  (5.0  x 
10  L/mol),  and  QTc  and  QAib  are  the  normal  physiological  concentrations  of  transcortin 
and  albumin,  0.7  umol/L  and  550  umol/L,  respectively  (Ballard  1979b).  This  expression 
may  be  rearranged  to  give  a  quadratic  equation  in  terms  of  Cf,  the  free  concentration. 
After  solving  the  quadratic  for  Cf,  the  fraction  unbound,  fu,  may  then  be  determined  by 


fu   = 


r      •  Eq.  3-2 

*"  Total 
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In  order  to  determine  the  concentrations  at  which  transcortin  saturation  occurs,  values  of 
fu  were  calculated  for  a  large  concentration  range  using  these  formulae  and  the  nonlinear 
binding  is  shown  in  Figure  3-1 .  Total  and  free  concentrations  for  a  smaller  range  of 
concentrations  are  provided  in  a  separate  plot  in  which  nonlinear  binding  is  compared  to 
a  constant  binding  of  95%. 


Unbound  fraction  of  Hydrocortisone 

0.3 


0.25 

0.2 

f„  0.15 

0.1 
0.05 

0 


10     100    1000   10000  100000 
Cp  (ng/mL) 


Figure  3-1 .  Nonlinear  binding  of  hydrocortisone  to  plasma  proteins. 
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Values  of  the  unbound  fraction  (fu)  were  calculated  from  total  concentrations  and  average 
plasma  protein  concentrations  and  binding  affinities  according  to  Eq.  3-1  and  Eq.  3-2.  For 
comparison,  free  concentrations  were  also  predicted  using  a  linear  binding  of  95%. 


For  concentrations  less  than  100  ng/mL,  the  fraction  unbound  is  relatively 
constant,  roughly  0.05  (5%).  This  is  apparent  in  the  second  plot  in  which  free 
concentrations  were  determined  with  both  nonlinear  and  linear  binding:  The  predicted 
free  levels  are  equal  for  total  hydrocortisone  concentrations  below  100  ng/mL.  However, 
for  larger  concentrations,  the  nonlinearity  is  pronounced  and  the  concentration  of  free 
drug  rises  sharply  with  increasing  total  concentrations,  reaching  a  plateau  of  just  over 
0.25  (25%)  for  concentrations  above  10,000  ng/mL  (10  ug/mL),  the  concentration  at 
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which  transcortin  binding  sites  are  saturated  and  hydrocortisone  exhibits  linear  binding  to 
albumin.  Thus,  for  a  concentration  between  0.1  and  10  ug/mL,  plasma  protein  binding  is 
nonlinear  and  the  pharmacologically  active  free  fraction  is  concentration-dependent. 
Pharmacokinetic  Considerations 

It  has  been  suggested  that  nonlinear  protein  binding  leads  to  dose-dependent 
pharmacokinetics  of  hydrocortisone  (Toothaker  and  Welling  1982;  Derendorf,  Mollmann 
et  al.  1991).  The  same  is  true  for  prednisolone,  which  binds  to  both  albumin  and 
transcortin  (Boekenoogen,  Szefler  et  al.  1983;  Rohdewald,  Rehder  et  al.  1987;  Rohatagi, 
Hochhaus  et  al.  1995).  Nonlinear  binding  may  modulate  the  pharmacokinetics  of 
hydrocortisone  in  two  ways. 

Clearance.  Perhaps  the  most  obvious  way  by  which  pharmacokinetics  are 
dependent  on  plasma  protein  binding  is  by  affecting  the  clearance  and,  in  fact,  increased 
clearance  values  have  been  observed  with  higher  doses.  The  total  body  clearance  of 
hydrocortisone  after  a  20  mg  IV  bolus  dose  (n=8)  was  reported  as  18.2  ±  4.2  L/h 
(Derendorf,  Mollmann  et  al.  1991).  Noncompartmental  analysis  of  unpublished  data  after 
100  mg  hydrocortisone  IV  (n=5)  produced  a  larger  total  body  clearance,  21.9  ±  3.0  L/h. 
These  clearance  values  indicate  that  hydrocortisone  is  a  low  extraction  drug  with  a 
variable  hepatic  extraction  ratio.  For  low  extraction  drugs,  the  Wilkinson-Shand  equation 
for  hepatic  clearance  (Wilkinson  and  Shand  1974), 

may  be  simplified  to 

CL  =  CLm  ■  fu  Eq.  3-4 
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where  QH  is  hepatic  blood  flow,  fu  is  the  unbound  fraction  of  the  drug,  and  CLint  is  the 
intrinsic  clearance,  a  measure  of  enzyme  activity.  According  to  Eq.  3-4,  it  is  obvious  that 
if  fu  increases,  then  CL  will  increase  proportionally.  Of  course,  determining  CL  for 
hydrocortisone  presents  a  complex  problem  because  fu  is  not  constant  for  concentrations 
between  0.1  and  10  ug/mL.  To  illustrate  the  dependence  of  total  CL  on  fu,  clearance 
values  were  predicted  for  a  large  range  of  total  hydrocortisone  concentrations.  Values  of 
CL  were  calculated  using  Eq.  3-3  with  QH  of  90  L  and  CLint  of  365  L/h  (estimated  using  a 
total  CL  of  1 8.2  L/h  and  a  free  fraction  of  0.05,  the  value  at  low  concentrations).  The 
predicted  clearance  values  are  plotted  versus  hydrocortisone  concentrations  in  Figure  3-2. 
In  order  to  illustrate  the  changes  in  clearance  at  lower  concentrations,  a  second  plot  is 
provided  with  concentration  expanded  logarithmically.  The  predicted  clearance  at  10,000 
ng/mL  (10  ug/mL)  is  about  four  times  greater  than  that  at  100  ng/mL  (0.1  u-g/mL). 
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Figure  3-2.  Predicted  clearance  of  hydrocortisone  at  various  concentrations 
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Clearance  values  were  calculated  using  Eq.  3-3.  The  value  of  intrinsic  clearance  (CLint) 
was  estimated  from  the  published  CL  value  of  1 8.2  L/h  (Derendorf  Mollmann  et  al 
1991)  and  fu  of  0.05. 

Volume  of  distribution.  The  second  manner  in  which  nonlinear  binding  may 
modulate  the  pharmacokinetics  of  hydrocortisone  is  by  affecting  the  volume  of 


45 


distribution.  By  definition,  the  volume  of  distribution  is  the  volume  that  relates  the 
amount  of  drug  in  the  body  to  the  concentration  found  in  the  plasma.  Assessing  the 
changes  in  this  pharmacokinetic  parameter  is  somewhat  complex  because  most  drugs, 
including  hydrocortisone,  exhibit  a  distribution  phase  after  IV  bolus  administration  and, 
consequently,  several  volume  terms  may  be  defined.  Two  of  these,  the  volume  of 
distribution  of  the  central  compartment  (Vc)  and  the  steady-state  volume  of  distribution 
(Vdss)  are  discussed  below. 

Both  perfusion  and  permeability  control  drug  distribution.  Since  most  therapeutic 
corticosteroids  are  lipophilic  and  unionized  such  that  membranes  to  not  present  a  barrier 
for  distribution,  the  distribution  of  these  drugs  is  perfusion-limited  rather  than 
permeability-limited.  If  drug  distribution  is  related  to  blood  flow,  then  highly  perfused 
organs  and  tissues  such  as  the  liver  and  kidney  should  be  in  rapid  distribution  equilibrium 
with  the  blood  (or  plasma).  The  plasma  and  well-perfused  tissues  are  often  treated  as  a 
homogenous  unit  generally  referred  to  as  the  central  compartment.  This  is  shown  in  the 
distribution  diagram,  Figure  3-3.  Kinetic  homogeneity  does  not  necessarily  mean  that  the 
drug  concentrations  in  all  tissues  of  the  central  compartment  are  the  same  at  a  given  time 
point.  However,  it  is  assumed  that  any  changes  in  the  plasma  levels  will  be  reflected 
(quantitatively)  by  changes  in  central  compartment  tissue  levels  (Gibaldi  and  Perrier 
1982).  In  contrast,  drug  levels  in  poorly  perfused  tissues  require  a  finite  time  to  reach 
equilibrium  with  those  in  the  central  compartment.  Drug  levels  in  poorly  perfused  tissues 
first  increase,  reach  a  maximum,  and  then  decline.  After  some  time,  pseudo-equilibrium 
is  attained  between  the  highly  perfused  tissues  and  fluids  of  the  central  compartment  and 
the  poorly  perfused  tissues.  Distribution  of  drug  into  the  poorly  perfused  tissues  may 
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occur  at  different  rates  and,  for  this  reason,  poorly  refused  tissues  are  often  grouped 
together  into  a  single  peripheral  compartment.  As  shown  in  Figure  3-3,  the  hypothetical 
peripheral  compartment  is,  at  best,  a  hybrid  of  several  physiologic  units  (Gibaldi  and 
Perrier  1982). 


Central 


Peripheral 


IV  Bolus 
Dose 


Perfusion-limited  drug  distribution  may  be  described  in 
terms  of  physiologic  units  based  on  the  degree  of  perfusion 


Figure  3-3.  Perfusion-limited  drug  distribution. 


The  volume  of  distribution  of  the  central  compartment.  From  a  purely 
physiological  viewpoint,  the  initial  volume  of  distribution  after  an  IV  bolus  injection  is 
the  volume  of  plasma  (3.22  L  for  a  70  kg  man)  or,  if  the  drug  binds  to  erythrocytes,  the 
volume  of  whole  blood  (5.25  L  for  a  70  kg  man)  (Gibaldi  and  McNamara  1977;  Davies 
and  Morris  1993).  In  practice,  however,  there  is  normally  rapid  distribution  into  well- 
perfused  tissues  and  initial  distribution  volumes  are  larger  than  3.22  L  when  determined 
experimentally  with  a  practical  lower  limit  of  7-7.5  L  (Gibaldi  and  McNamara  1978; 
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Rowland  and  Tozer  1995).  (The  apparent  volume  of  plasma  proteins,  about  7.5  L  for 
albumin,  is  a  better  estimate  of  the  volume  of  distribution  for  drugs  highly  bound  to 
plasma  proteins  but  not  bound  to  tissues.  This  volume  term  is  larger  since  plasma 
proteins,  and,  therefore,  protein-bound  drugs,  equilibrate  between  plasma  and  other 
extracellular  fluids.  However,  if  the  drug  is  restricted  to  the  plasma  due  to  a  very  high 
molecular  weight,  the  volume  of  distribution  tends  to  approximate  the  volume  of  plasma, 
about  3  L,  Rowland  and  Tozer  1995).  The  volume  of  distribution  of  the  central 
compartment,  Vc,  relates  the  drug  concentration  immediately  after  an  IV  bolus  injection 
to  the  observed  plasma  concentration.  Of  course,  it  must  be  assumed  that  no  elimination 
has  taken  place  and  that  distribution  is  minimal.  A  Vc  of  23.9  ±  4.8  was  previously 
reported  for  hydrocortisone  after  a  20  mg  IV  bolus  dose  (Derendorf,  Mollmann  et  al. 
1991). 

Although  often  overlooked,  the  experimentally  determined  Vc  may  depend  on  the 
degree  of  binding  to  plasma  proteins.  The  reason  for  this  is  that,  despite  the  rapidity  of 
the  early  distribution  phase,  the  extent  of  distribution  into  the  well-perfused  tissues 
constituting  the  central  compartment  is  determined  by  the  ability  of  the  drug  to  distribute 
into  tissues,  the  degree  of  binding  in  the  well-perfused  tissues,  and  fu  (because  only 
unbound  drug  can  cross  membranes).  For  example,  if  two  drugs  differ  only  in  the  degree 
of  protein  binding,  the  drug  with  the  larger  fu,  i.e.,  the  drug  having  a  greater  free  fraction 
for  the  same  dose,  will  have  a  larger  Vc  after  IV  bolus.  In  the  present  case  of 
hydrocortisone  with  saturable  protein  binding,  a  larger  dose  yields  a  larger  f„  and  Vc  is 
dose-dependent.  Thus,  an  increase  in  the  dose  does  not  necessarily  produce  a  proportional 
increase  in  the  initial  concentration  observed  after  IV  bolus  administration. 
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The  volume  of  distribution  at  steady-state.  Possibly  the  best  approach  to 
assessing  the  effects  of  nonlinear  binding  on  drug  distribution  is  to  consider  the  steady 
state  volume  of  distribution,  Vdss.  By  definition,  Vdss  is  the  volume  of  distribution 
(during  distribution)  at  which  the  free  concentration  in  the  central  compartment  is  equal 
to  the  free  concentration  in  the  peripheral  compartment.  In  the  case  of  linear  kinetics,  the 
value  of  Vdss  is  between  Vc  and  Vdp,  the  volume  of  distribution  in  the  terminal  phase 
(Gibaldi  1969).  For  most  drugs,  Vdss  and  Vdp  differ  by  less  than  10%  (Gibaldi,  Levy  et 
al.  1978).  Although  this  volume  term  is  short-lived  in  its  ability  to  relate  the  amount  of 
drug  in  the  body  to  the  concentration  observed  in  the  plasma  after  a  single  dose,  under  the 
assumption  that  passive  diffusion  is  the  primary  means  of  drug  transport,  Vdss  is  the 
apparent  volume  of  distribution  during  a  constant  IV  infusion  at  steady  state.  The  Vdss  is 
easily  calculated  from  concentration-time  data  by  noncompartmental  analysis  and  may 
also  be  defined  in  terms  of  physiological  volumes  as 

v*n  =  v,+vt~  Eti-3-5 

J  ui 

where  Vp  is  the  plasma  volume  (central  compartment),  Vt  is  the  volume  of  tissue  water 
into  which  the  drug  distributes  (peripheral  compartment),  fu  is  the  fraction  unbound  in  the 
plasma,  and  fut  is  the  fraction  unbound  in  the  tissue.  From  this  expression,  Vdss  is  directly 
proportional  to  the  unbound  fraction  in  plasma  and  inversely  proportional  to  the  unbound 
fraction  in  peripheral  tissues.  The  Vt  depends  on  the  ability  of  the  drug  to  cross 
membranes  and  the  degree  to  which  drug  molecules  gain  access  to  intracellular  fluids 
(Rowland  and  Tozer  1995).  For  lipophilic  drugs  that  cross  membranes  easily,  Vt  is  often 
taken  to  be  total  body  water  minus  the  plasma  volume,  38.71  L  for  a  70  kg  man.  Without 
this  assumption,  fut  and  Vt  may  not  be  considered  independently  and  only  their  ratio  may 
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be  used  in  Eq.  3-5.  To  illustrate  how  varying  the  unbound  fraction  influences  the  steady 
state  volume  of  distribution,  Vdss  was  calculated  over  a  wide  range  of  fu  values  after 
defining  fut  as  0.3,  estimated  from  the  published  Vdss  of  33.7  L  for  hydrocortisone 
(Derendorf,  Mollmann  et  al.  1991)  and  assuming  fu  is  0.05.  See  Figure  3-4.  (This  is  a 
general  illustration  of  Vdss-dependency  on  fu  with  fu  varying  between  0.001  and  1.0  rather 
than  the  range  of  hydrocortisone  of  0.05  and  0.25  shown  in  Figure  3-1). 

If  fu  is  not  a  constant,  no  singular  value  of  Vdss  may  be  used  to  characterize  the 
distribution  of  the  drug.  Thus,  the  volume  of  distribution  during  a  constant  rate  infusion 
for  a  drug  with  nonlinear  protein  binding  depends  on  the  infusion  rate.  Under  the 
assumptions  used  to  generate  Figure  3-3,  the  Vdss  of  hydrocortisone  (0.05  <  fu  <  0.25) 
may  range  between  9  and  40  L. 
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Fifiure  3-4.  Steady  state  volume  of  distribution  with  changing  f„. 
Steady-state  volumes  of  distribution  (Vdss)  were  calculated  according  to  Eq.  3-5  using 
fut=0.3  and  Vt=38.71  L.  Values  corresponding  to  free  fractions  (fu)  of  0.05  and  0.25  are 
highlighted  to  indicate  the  range  of  possible  Vdss  values  for  hydrocortisone,  based  on 
nonlinear  binding  to  plasma  proteins. 
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Distribution  Kinetics  and  Compartmental  Models 

The  most  common  approach  to  the  pharmacokinetic  characterization  of  a  drug  is 
to  represent  the  body  as  a  system  of  compartments  and  to  assume  that  the  transfer  of  drug 
between  the  compartments  and  drug  elimination  may  be  described  using  first-order  rate 
constants  (linear  kinetics)  (Gibaldi  and  Perrier  1982).  The  simplest  model  depicts  the 
body  as  a  single,  homogeneous  unit.  Thus,  this  model  is  commonly  referred  to  as  the  one- 
compartment  body  model.  Portraying  the  body  as  a  single  compartment  may  be  useful  for 
laying  the  foundation  of  pharmacokinetics,  but  this  is  usually  an  inaccurate  representation 
for  the  events  that  follow  drug  administration.  The  basic  assumption  of  the  one- 
compartment  distribution  model  is  that  equilibration  of  drug  between  blood  and  tissues 
occurs  spontaneously  (instantaneously)  (Rowland  and  Tozer  1995).  In  reality,  however, 
distribution  takes  time.  As  mentioned  previously,  the  time  required  for  distribution 
depends  on  tissue  perfusion,  permeability  characteristics  of  the  drug,  and  the  partitioning 
between  tissue  and  blood.  To  account  for  this  time  factor  in  drug  distribution,  multi- 
compartments models  are  often  used,  the  two-compartment  model  being  the  most  favored 
one  based  on  physiologic  considerations  (Rowland  and  Tozer  1995). 

The  two-compartment  model.  Although  a  one-compartment  model  has  been  used 
to  describe  the  concentration-time  profile  of  hydrocortisone  after  oral  administration 
(Toothaker,  Sundaresan  et  al.  1982),  a  two-compartment  model  was  required  to  describe 
clinical  data  after  IV  bolus  administration  (Toothaker  and  Welling  1982;  Derendorf, 
Mollmann  et  al.  1991).  This  was  due  to  the  pronounced  distribution  phase  observed  after 
IV  administration.  For  a  single  IV  bolus  dose,  the  two -compartment  model  may  be 
written  as  a  simple  biexponential  curve, 
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-fit 


C(t)  =  A-e-°'+B-e^  Eq.  3.6 

where  a  and  (3  are  macroconstants  characterized  by  the  first-order  rate  constants  for  drug 
transfer  out  of  and  into  the  central  compartment,  k,2  and  k2i,  respectively,  and  the  first- 
order  rate  constant  for  drug  elimination  from  the  central  compartment,  k10, 


a 


P 


_  (*10  +  ku  +  k2l)  +  {(k]0  +  ku  +  *2, )2  - 4 •  kw  ■  k2X }~2 

2 

_  (*10  +  kn  +k2X)- {(k]0  +kn+k2l)2-4- kw  •  k2i }^ 


and  A  and  B  are  pre-exponential  factors  defined  in  terms  of  the  dose,  Vc,  and  k2i, 
D    (k2X-a) 


A  = 


B  = 


Vc     {p-a) 


Vc     (a-fi) 


A  diagram  of  this  model  is  shown  in  Figure  3-5. 


Compartment  1 


Compartment  2 


■10 


Y  (  1st  order  elimination  rate  constant) 


Figure  3-5.  Two-compartment  pharmacokinetic  model. 
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The  two-compartment  model  adequately  described  the  concentration-time  profiles 
of  hydrocortisone  for  the  dose  and  dosage  forms  considered  (Derendorf,  Mollmann  et  al. 
1991).  However,  because  of  the  nonlinear  processes  taking  place,  parameters  obtained 
from  fitting  data  to  this  model  must  be  regarded  as  dose-specific  and  may  not  be  used  to 
predict  the  disposition  of  hydrocortisone  for  other  doses.  This  is  especially  true  if  the 
distribution-  or  a-phase  is  to  be  accurately  described  since  saturable  protein  binding  is 
most  apparent  at  the  high  concentrations  achieved  immediately  after  injection. 

Compartmental  models  and  nonlinearitv.  After  considering  the  potential  dose- 
dependency  of  both  CL  and  Vd,  new  models  were  tested  in  which  Eq.  3-3  was  used  to 
describe  the  clearance  and  Eq.  3-5  was  used  to  describe  Vd  as  each  of  these  equations 
allows  fu  to  be  treated  as  a  separate  variable.  These  models  were  based  on  a  one- 
compartment  system, 

C(0  =  £■*-*«', 

where  D  is  the  dose,  Vd  is  the  volume  of  distribution,  and  ke  is  the  rate  of  elimination, 

Qn-CLmx.fu 

k  _CL_QH+CLmi.fu 

K-y- 7—.  Eq.3-7 

<i         v  +V  •  — 

P  1         r 

J  ul 

Using  the  one-compartment  model  allowed  distribution  to  be  described  in  terms  of  the 
nonlinear  fu.  (The  two-compartment  model  is  formally  defined  only  in  terms  of  Vc,  the 
volume  of  the  central  compartment).  Although  Eq.  3-5  for  Vdss  is  only  true  at  steady-state 
for  pharmacokinetic  models  with  distribution,  under  the  assumption  that  a  dynamic 
equilibrium  is  established  immediately  after  drug  administration,  Vdss  is  equal  to  Vd  in  a 
one  compartment  system.  By  incorporating  the  physiologically  based  expressions  for  CL 
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and  Vd,  intrinsic  clearance  and  overall  tissue  binding  for  hydrocortisone  may  be 
determined  from  concentration-time  data. 

This  model  allowed  for  any  deviation  from  linear  pharmacokinetics  due  to 
changes  in  fu,  the  data  were  well-described  in  the  elimination  phase,  and  this  model 
offered  a  slight  improvement  over  the  usual  one-compartment  model  having  constant 
values  of  CL  and  Vd.  However,  early  time  points  were  severely  underestimated  due  to  the 
lack  of  any  distribution  phase.  Several  other  models  were  constructed  by  starting  with  a 
one-compartment  system  yet  allowing  Vd  to  expand.  These  pseudo-one-compartment 
models  dramatically  improved  fitting  of  the  data  and  although  they  were  interesting 
conceptually,  they  were  somewhat  empirical,  based  more  on  intuition  than  on  sound 
physical  or  mathematical  principles. 

The  two-compartment  model  was  revisited  with  the  first  order  elimination  rate 
constant,  kio,  defined  by  the  Wilkinson-Shand  equation  for  hepatic  clearance  and  the 
physiological  expression  for  volume  of  distribution,  Eq.  3-3  and  Eq.  3-5,  respectively. 
This  model  performed  well  for  individual  data  sets  but  the  same  distribution  rate 
constants,  k)2  and  k2i,  could  not  be  used  to  describe  hydrocortisone  after  different  doses 
(data  after  20  and  100  mg  IV  bolus  doses). 

Although  some  of  these  models  may  have  worked  in  practice,  they  were  all 
flawed  in  theory  since  nonlinear  conditions  were  imposed  on  functions  derived  under  the 
assumption  that  all  processes,  namely  distribution  and  elimination,  are  strictly  linear.  It 
was  concluded  that  Eq.  3-3  was  useful  for  describing  nonlinear  elimination  but  that  the 
distribution  characteristics  were  too  complex  to  be  fully  described  using  traditional 
pharmacokinetic  models.  Ultimately,  a  new  model  was  derived  based  on  the  distribution 
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of  free  drug  by  passive  diffusion  that  allowed  for  nonlinear  processes. 

Purpose 
The  objectives  of  this  project  were  as  follows: 

1 .  To  investigate  the  possible  nonlinear  pharmacokinetics  of  hydrocortisone; 

2.  To  determine  whether  or  not  such  nonlinear  processes  are  important  for  describing 
the  concentration-time  profile  of  endogenous  Cortisol  in  the  physiological 
concentration  range;  and 

3.  To  construct  a  physiological  pharmacokinetic  model  that  could  produce  good  fits  of 
concentration-time  profiles  as  well  as  provide  some  insight  as  to  how  nonlinear 
binding  affects  these  profiles. 

Methods 

Data  from  a  previously  published  clinical  study  of  hydrocortisone  (single  IV 
bolus  doses  of  20  mg  administered  to  eight  volunteers)  were  obtained  (Derendorf, 
Mollmann  et  al.  1991).  Unpublished  data  after  single  IV  bolus  doses  of  100  mg  to  five 
volunteers  were  also  obtained. 
Preliminary  Calculations 

Before  compartmental  analysis,  preliminary  calculations  were  performed  to 
determine  the  variability  in  the  unbound  fraction  (fuj)  and  other  parameters  such  as  the 
volume  of  distribution  using  concentration-time  data  after  a  20  mg  IV  bolus  dose  (n=8). 
(In  the  remainder  of  this  chapter,  the  unbound  fraction  of  drug  in  plasma  is  denoted  ful, 
designating  the  central  compartment,  compartment  1 ,  as  a  subscript.  The  overall  unbound 
fraction  in  tissues  is  denoted  fu2  rather  than  fut,  indicating  the  second,  or  peripheral, 
compartment.  This  notation  is  used  to  distinguish  between  compartment-  and  time- 
dependent  parameters,  t  being  the  most  appropriate  symbol  for  time). 

In  order  to  test  the  hypothesis  that  nonlinear  binding  greatly  affects  the  initial 
distribution  phase  of  hydrocortisone,  values  of  fui  were  calculated  for  each  time  point  for 
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each  subject,  averaged,  and  plotted  for  each  time  point  with  standard  deviations.  The 
volume  of  distribution  was  also  calculated  for  each  time  point,  averaged,  and  plotted.  By 
definition,  the  volume  of  distribution  relates  the  amount  of  drug  in  the  body  at  any  given 
time  (Xt)  to  the  concentration  observed  in  the  plasma  (Cp), 

which  may  be  rearranged  to  give 

y.'-frfa-  Eq.3-8 

where  V,  is  the  apparent  volume  of  distribution  at  time  t.  By  relating  the  amount  of  drug 
eliminated  at  time  t  to  the  area  under  the  curve  from  time-zero  to  time  t  (AUCt),  Eq.3-8 
may  be  expressed  as 

Y_D.jAUC.-AUC,) 

Cp(t)-AUCn  Eq'3'9 

where  AUCoo  is  the  total  AUC  (extrapolated  to  infinity)  (Chiou  1972;  Niazi  1976a;  Niazi 

1976b;  McNamara,  Slattery  et  al.  1979). 

Compartmental  Models  with  Linear  and  Nonlinear  Pharmacokinetics 

Individual  and  average  hydrocortisone-time  data  after  IV  administration,  20  mg 
hydrocortisone  to  eight  patients  and  100  mg  to  five  patients,  were  fitted  to  different 
pharmacokinetic  models.  The  two-compartment  model,  Eq.  3-6  using  microconstants, 
was  compared  to  a  new  two-compartment  model  having  free  concentration  gradient- 
based  distribution.  (When  the  conventional  two-compartment  model  was  used,  the  dose, 
20  mg  or  100  mg,  was  factored  out  of  the  pre-exponential  parameters,  A  and  B,  and 
entered  as  a  constant  for  each  data  set).  A  complete  derivation  and  discussion  of  the 
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gradient-based  distribution  model  is  provided  in  Appendix  A.  In  integrated  form  (suitable 
only  for  linear  binding  and  linear  pharmacokinetics),  this  model  is 

Cx(t)  =  A-e-a'+B-e-<"  Eq.  3_10 

where 


D- 

A  = 


v 

V        V2 


Vr(p-a)  E*3-U 


and 


D- 


V2 


B=      vAa-P)      •  Eq-3'12 

In  these  equations,  D  is  the  dose,  K  is  the  diffusion  coefficient,  Vj  is  the  volume  of  the 
central  compartment  (formally  defined  as  the  plasma  volume,  3.22  L),  and  fu2  is  the  free 
fraction  of  drug  in  the  periphery.  The  rate  constants  a  and  [3  are  functions  of  the  free 
fraction  in  the  central  compartment  (ful),  fu2,  V,,  and  V2,  the  volume  of  the  peripheral 
compartment  (formally  defined  as  38.71  L).  From  viewing  the  integrated  form  of  this 
new  model,  it  is  apparent  that  there  is  a  close  structural  similarity  to  the  conventional 
two-compartment  model.  In  fact  k2i,  the  rate  constant  in  the  conventional  two- 
compartment  model  describing  the  transfer  of  drug  into  the  central  compartment  from  the 
peripheral  compartment  is  present  in  the  gradient-based  distribution  model.  However,  in 
the  latter,  it  is  defined  in  terms  of  the  diffusion  coefficient,  tissue  binding,  and  the  volume 
of  the  peripheral  compartment.  From  Eq.  3-6  and  Eq.  3-10  through  Eq.  3-12, 


k2x 

v2 


Likewise, 
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*12 


K-fui 


Vy 


This  similarity  regarding  k,2  is  found  in  the  expressions  included  in  Appendix  A. 

Due  to  the  fact  that  ful  is  not  constant  in  the  present  application  to  hydrocortisone 
data,  no  integrated  solution  for  the  gradient-based  distribution  model  is  possible  in  closed 
form.  The  differential  form  was  used  for  all  hydrocortisone  data  fitting  procedures.  The 
differential  equations  are 


dCx 

dt 


r,       ") 


^—-C,  Eq.3-13 


and 


dC7      K-f.  K-f, 


i         7} c2-  Eq.  3-14 


dt  V2         '         V2 

where  d  is  the  concentration  in  the  central  compartment  and  C2  is  the  concentration  in 
the  peripheral  compartment.  Values  of  ful  were  calculated  at  a  given  concentration  using 
Eq.  3-1  and  Eq.  3-2.  The  values  of  K  and  fu2  were  variables  determined  in  fitting  the  data 
to  the  model.  The  elimination  rate  from  the  central  compartment  (kio)  was  variable  and 
defined  in  terms  of  Eq.  3-3  and  V]  as 

QH-CLmrfu] 

k    _CL_QH+CLmrfuX 

10  ~T v Eq-3-15 

The  intrinsic  clearance,  CL,nt,  was  determined  in  fitting  the  data  to  the  model  and  QH  was 
fixed  at  an  average  hepatic  blood  flow  of  90  L  (Rowland  and  Tozer  1995). 

The  gradient-based  distribution  model  was  incorporated  into  an  Excel™ 
spreadsheet  to  allow  simulations  to  be  performed.  The  model  was  entered  in  differential 
form  and  is  integrated  numerically  over  2000  intervals.  This  allows  nonlinear  processes, 
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either  clearance  or  distribution,  to  be  simulated.  Two  sets  of  input  parameters  may  be 
entered  to  generate  two  concentration-time  profiles  for  each  simulation.  This 
comprehensive  spreadsheet  includes  (1)  variable  simulation  times,  (2)  on/off  functions 
for  linear  and  nonlinear  binding  to  plasma  proteins  (albumin  and  transcortin)  with 
variable  plasma  concentrations  and  binding  affinities,  (3)  on/off  functions  for  hepatic 
and/or  renal  clearance  with  variable  hepatic  and/or  renal  function,  (4)  plots  for  plasma 
and  tissue  concentrations  (total  and  free),  (5)  plots  for  clearance,  volume  of  distribution 
(calculated  at  each  time  point  using  a  noncompartmental  algorithm),  and  the  unbound 
fraction  in  plasma,  and  (6)  a  complete  noncompartmental  analysis  of  the  simulated 
concentration-times  profiles.  The  current  version  of  this  spreadsheet  is  limited  to  IV 
bolus  administration  and  tissue  binding  is  assumed  to  be  linear  for  all  doses. 
Simulations  with  the  Gradient-Based  Distribution  Model 

After  fitting  clinical  data  (hydrocortisone,  20  mg  IV  to  eight  volunteers)  to  the 
gradient-based  distribution  model  in  the  SCIENTIST  program,  the  fitted  parameters  were 
entered  into  the  numerically  integrated  model  in  the  Excel™  spreadsheet  in  order  to 
conduct  simulations  specifically  for  hydrocortisone.  Of  particular  interest  was  the  time- 
dependent  volume  of  distribution.  Simulated  Vd-time  profiles  were  compared  to  those 
generated  from  the  actual  data  using  the  noncompartmental  algorithm  in  Eq.  3-9. 
Endogenous  Cortisol 

The  pharmacokinetics  of  endogenous  Cortisol  were  also  investigated  for  future 
incorporation  into  the  PK/PD  model  describing  the  effects  of  exogenous  corticosteroids 
on  blood  lymphocytes.  The  circadian  rhythm  of  Cortisol  was  described  using  a  linear 
release  model  (Rohatagi,  Bye  et  al.  1996;  Derendorf,  Mollmann  et  al.  1997), 
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where  the  release  rate  of  Cortisol  is  described  by  two  linear  functions  (Re)  between  the 
maximum  Cortisol  release  (Rmax)  at  tmax  (6-10  AM)  and  a  minimum  release  at  tmin  (8  PM- 
2  AM),  ke  ort  is  the  elimination  rate  of  Cortisol,  and  Vd  is  the  volume  of  distribution 
(Rohatagi,  Hochhaus  et  al.  1995).  Cortisol  concentration-time  data  after  placebo  (n=12) 
were  fitted  to  a  modified  linear  release  model  with  gradient  based  distribution  and 
nonlinear  elimination.  This  modified  model  was  based  on  Eq.  3-16  through  Eq.  3-18  with 
the  assumption  that  Cortisol  is  released  into  the  central  compartment  (Vd=Vc=3.22  L)  and 
eliminated  by  hepatic  metabolism  (keCort=ki0)  according  to  Eq.  3-13  through  Eq.  3-15. 
For  comparison,  the  same  data  set  was  fitted  to  the  linear  release  model  (one- 
compartment)  used  in  previous  publications  with  a  fixed  volume  of  distribution  (33.7  L) 
and  a  first-order  elimination  rate  constant,  keCort  (Rohatagi,  Bye  et  al.  1996). 
Data  Fitting 

All  data  fitting  procedures  were  performed  using  SCIENTIST  (Micromath,  Salt 
Lake  City,  Utah),  a  program  that  allows  models  to  be  entered  in  either  analytic  or 
differential  form.  The  integrity  of  the  models  was  compared  using  the  Model  Selection 
Criteria  (MSC),  a  modified  Akaike  information  criterion  generated  in  the  statistics  option 
in  this  program.  In  general,  a  larger  MSC  value  is  indicative  of  a  better  description  of  the 
data.  The  complexity  of  the  model  and  the  number  of  parameters  needed  to  construct  the 
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model  are  also  factors  in  determining  the  MSC. 

Results  and  Discussion 
Preliminary  Calculations 

The  results  of  the  preliminary  calculations  of  the  unbound  fraction  in  plasma  (ful) 
for  hydrocortisone  (HC)  data  after  20  mg  (single  IV  bolus  doses  to  eight  volunteers)  are 
shown  in  Figure  3-6. 


Free  fraction  (fu)  after  20mg  HC  (IV) 


0.25 


Figure  3-6.  Free  fraction  (f„i)  of  hydrocortisone  in  plasma. 


From  comparing  these  values  to  those  in  Figure  3-1,  it  is  apparent  that  complete 
saturation  of  transcortin  does  not  occur  after  a  20  mg  dose,  fu]  not  having  reached  the 
plateau  of  0.25.  At  roughly  five  hours  after  injection,  a  sufficient  amount  of 
hydrocortisone  has  been  eliminated  to  allow  a  dynamic  equilibrium  to  be  established 
between  free  and  bound  drug  and  fui  becomes  relatively  constant  at  just  under  0.05. 

The  volume  of  distribution  at  various  time  points  is  shown  in  Figure  3-7.  In  linear 
kinetics,  the  volume  of  distribution  rises  from  the  initial  value,  Vc,  to  that  in  the 
elimination  phase,  Vdp,  with  the  steady  state  value  situated  somewhere  between  Vc  and 


61 


Vdp.  This  is  clearly  not  the  case  for  this  dose  of  hydrocortisone  wherein  Vd  rises  and  then 
declines  toVdp.  The  same  trend  was  observed  when  Vd  was  simulated  using  the  gradient- 
based  distribution  model  and  parameters  obtained  from  the  analysis  of  clinical  data  (20 
and  100  mg  hydrocortisone,  IV  bolus  administration).  The  rise  in  Vd  may  be  explained  by 
the  fact  that  the  large  initial  fu,  allows  for  rapid  distribution  of  free  drug  into  the 
periphery  where  tissue  binding  occurs.  Based  on  Eq.  3-8,  there  is  a  disproportion^ 
decline  in  plasma  concentration  in  comparison  to  the  amount  of  drug  remaining  in  the 
body,  which  may  be  described  only  by  an  increase  in  the  volume  of  distribution. 
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Simulated  Vd  after  20mg  HC  (IV) 
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Figure  3-7.  Calculated  volume  of  distribution  after  20  me  hydrocortisone. 
For  the  plot  on  the  left,  a  noncompartmental  algorithm  (Eq.  3-9)  was  used  to  calculate  the 
volume  of  distribution  (Vd)  at  each  data  time  point  after  the  IV  bolus  administration  of  20 
mg  hydrocortisone  (n=8).  These  were  averaged  and  plotted  along  with  the  standard 
deviation.  For  the  plot  on  the  right,  the  volume  of  distribution  was  calculated  using  the 
same  algorithm  for  a  simulated  concentration-time  profile.  This  simulation  was 
performed  using  the  gradient-based  distribution  model  (Eq.  3-13)  with  parameters 
obtained  from  fitting  clinical  data  to  the  model  (20  mg  doses,  n=8). 

A  large  unbound  fraction  also  yields  a  rapid  hepatic  elimination  at  early  time 

points.  The  influence  of  fui  on  CL  was  illustrated  in  Figure  3-2.  After  reaching  steady 

state,  the  time  at  which  free  levels  in  the  central  and  peripheral  compartments  are  equal, 
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drug  begins  to  flux  back  into  the  central  compartment  where  it  is  tightly  bound  to 
transcortin.  The  amount  of  drug  removed  from  the  central  compartment  is 
counterbalanced  by  the  amount  of  drug  returning  from  the  peripheral  compartment  and, 
for  a  brief  period,  plasma  concentrations  are  relatively  constant.  In  fact,  expanding  the 
fui-time  profile  logarithmically  reveals  a  slight  plateau  corresponding  to  the  time  at  which 
the  peak  Vd  is  observed.  This  plot  is  provided  in  Figure  3-8. 


Log-Log  plot  of  the  free  fraction  (fu) 
after  20mg  hydrocortisone 
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Figure  3-8.  Free  fraction  of  hydrocortisone  in  plasma  (f,,,)  expanded  logarithmically. 
The  unbound  fraction  (fui)  of  hydrocortisone  was  calculated  from  clinical  data  (n=8)  and 
population  average  values  of  protein  concentrations  and  binding  affinities.  These  were 
averaged  for  each  time  point  and  plotted  along  with  the  standard  deviation.  Values  of  ful 
were  relatively  constant  at  time  points  around  0.5  h. 


In  this  scenario,  there  is  a  disproportional  decline  in  the  amount  of  drug  remaining 
in  the  body  and  the  plasma  concentration  leading  to  a  smaller  apparent  volume  of 
distribution.  It  may  then  be  hypothesized  that  the  binding  of  hydrocortisone  to  plasma 
proteins,  specifically  transcortin,  is  much  greater  than  its  affinity  to  tissues.  Similar  trends 
were  observed  when  the  100  mg  data  was  subjected  to  this  type  of  analysis.  Again,  there 
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was  a  slight  plateau  at  approximately  0.5  h  with  very  little  change  in  ful  from  0.3  to  0.7  h 
Compartmental  Models  with  Linear  and  Nonlinear  Pharmacokinetics 

The  results  of  simultaneously  fitting  hydrocortisone  data  after  a  20  mg  or  100  mg 
IV  bolus  dose  to  a  two-compartment  model  are  shown  in  Table  3-1 .  This  model  is 
described  in  Eq.  3-6  in  terms  of  macroconstants.  These  results  may  be  compared  to  those 
obtained  from  fitting  the  same  data  sets  to  the  new  two-compartment  model  described  in 
Eq.  3-13  and  Eq.  3-14.  Parameter  values  obtained  from  the  gradient-based  distribution 
model  are  shown  in  Table  3-2. 


Table  3-1.  Parameters  obtained  from  simultaneously  fitting  20  mg  andlOO  mp 

hydrocortisone  data  to  a  two-compartment  model  with  macroconstants. 
Parameter        Value         SD* 


A  (ng/ml) 

24.9 

2.6 

B  (ng/ml) 

16.0 

0.7 

a(l/h) 

6.2 

0.8 

P(l/h) 

0.35 

0.02 

MSC 

2.7 

♦Standard  deviations  were  generated  using  the  statistics  option  in  SCIENTIST.  Clinical 
data  after  20  and  100  mg  hydrocortisone  (IV  bolus  administration)  were  fitted 
simultaneously  to  a  two-compartment  model  (Eq.  3-6).  A  single  set  of  parameters 
obtained  since  the  dose  was  entered  as  a  separate  constant,  20  or  100  mg. 


was 


Table  3-2.  Parameters  obtained  from  simultaneously  fitting  20  mg  and  100  mp 
hydrocortisone  data  to  a  two-compartment  model  with 
gradient-based  diffusion- 
Parameter  Value  SD* 


K(l/h) 

223.7 

10.4 

CLint(L/h) 

93.6 

4.1 

fU2 

0.30 

0.01 

♦Standard  deviations  were  generated  using  the  statistics  option  in  SCIENTIST.  Clinical 
data  after  20  and  100  mg  hydrocortisone  (IV  bolus  administration)  were  fitted 
simultaneously.  A  single  set  of  parameters  was  obtained  since  the  dose  was  entered  as  a 
separate  constant,  20  or  100  mg. 


A  larger  MSC  value  was  obtained  using  the  gradient-based  distribution  model  4.2 
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compared  to  2.7  for  the  normal  two-compartment  model.  When  the  20  mg  and  100  mg 
data  were  analyzed  separately  (simultaneously  fitting  8  data  sets  after  20  mg  and  5  data 
sets  after  100  mg),  the  differences  in  MSC  values  between  the  two  models  were  even 
greater,  at  6.1  and  5.2,  respectively,  for  the  gradient-based  distribution  model  compared 
to  2.4  and  2.3  for  a  conventional  two-compartment  model. 

Visibly  better  fits  were  obtained  using  the  gradient-based  distribution  model, 
particularly  at  early  time  points.  This  is  shown  in  both  Figure  3-9  and  Figure  3-10.  The 
latter  is  an  enlargement  of  the  fitted  100  mg  data  illustrating  that  a  brief  plateau,  between 
0.5  and  1.5  h,  precedes  the  elimination  phase. 
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Figure  3-9.  Fitted  hydrocortisone  data  after  20  mg  and  1 00  mg. 
Hydrocortisone  data  were  fitted  to  the  conventional  two-compartment  model  Eq  3-6  and 
a  new  distribution  model,  Eq.  3-13  through  Eq.  3-15,  based  on  the  flux  of  free  drug  down 
a  concentration  gradient. 
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The  conventional  two-compartment  model  was  able  to  accurately  describe  part  of 
the  concentration-time  profile  of  the  100  mg  data  set.  However,  as  seen  in  Figure  3-10, 
the  predicted  profile  deviates  from  the  actual  data  values  in  the  post-distribution  phase. 


Conventional  2-Compartment  Model 
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Figure  3-10.  Fitted  hydrocortisone  data  after  100  mg  IV  using  a  two-compartment  model 

with  gradient-based  distribution. 
Hydrocortisone  data  were  fitted  to  the  conventional  two-compartment  model  Eq.  3-6  and 
a  new  distribution  model,  Eq.  3-13  through  Eq.  3-15,  based  on  the  flux  of  free  drug  down 
a  concentration  gradient. 


The  deviation  observed  using  the  conventional  two-compartment  model  is  not 
observed  in  the  fitted  profile  using  the  gradient-based  distribution  model.  In  general, 
either  the  20  or  1 00  mg  data  could  be  described  using  a  conventional  two-compartment 
model  using  weight  factors  but  both  could  not  be  described  simultaneously  with  the  same 
set  of  parameters.  In  addition  to  describing  the  data,  this  model  provided  an  estimate  of 
the  intrinsic  clearance  of  hydrocortisone  and  the  overall  tissue  binding,  principle 
objectives  in  this  study.  The  unbound  fraction  in  the  peripheral  compartment  (fu2)  of  0.3 
indicates  that  the  affinity  of  hydrocortisone  to  tissues  is  roughly  the  same  as  that  to 
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albumin.  Interestingly,  this  value  had  been  predicted  in  the  preliminary  calculations  using 
the  previously  reported  volume  of  distribution  at  steady-state  (Vdss)  of  33.7  L  and  a 
fraction  unbound  in  plasma  of  0.05  (for  low  concentrations).  Using  this  value  for  fu2  and 
0.26  for  ful  (the  unbound  fraction  of  hydrocortisone  in  plasma  after  total  saturation  of 
transcortin),  the  predicted  Vdss  is  36.7  L  (calculated  from  Eq.  3-5).  This  is  very  close  to 
the  maximum  observed  in  the  time-dependent  volume  of  distribution  in  Figure  3-6.  It 
may  be  hypothesized  that  the  maximum  Vd  in  Figure  3-6  corresponds  to  the  steady-state 
volume  of  distribution  (formally  defined  as  the  Vd  at  the  time  when  free  concentrations  in 
the  central  and  peripheral  compartments  are  equal). 

Although  diffusion  coefficients  are  defined  in  different  terms  depending  on  the 
application,  the  value  of  K  determined  in  this  pharmacokinetic  model  is  certainly  related 
to  the  membrane  permeability  of  hydrocortisone  and  logP  (Chapter  2).  As 
discussed  in  Appendix  A,  the  diffusion  coefficient  in  the  gradient-based  distribution 
model  has  units  of  volume-per-time.  For  this  reason,  K  may  be  thought  of  as  the 
clearance  of  the  drug  from  the  peripheral  compartment,  assuming  drug  metabolism 
occurs  only  in  the  central  compartment.  With  knowledge  of  the  physiological  parameters 
determined  in  the  analysis  of  clinical  data  (diffusion  coefficients,  enzyme  activity,  and 
tissue  binding)  and  the  affinity  constants  for  plasma  proteins,  it  should  be  possible  to 
accurately  predict  the  pharmacokinetic  disposition  of  hydrocortisone  all  doses. 
Simulations  with  the  Gradient-Based  Distribution  Model 

To  illustrate  the  usefulness  of  the  gradient-based  distribution  model,  simulations 
were  performed  for  two  different  scenarios.  These  are  shown  in  Figure  3-1 1  and  3-12  and 
the  user  interface  of  the  Excel™  spreadsheet  is  shown  in  Figure  3-13. 
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Figure  3-11.  Simulated  concentration-time  profiles  for  different  values  of  K. 
These  profiles  were  generated  using  the  gradient-based  distribution  model.  The  diffusion 
coefficient  for  Casel  was  1  and  that  for  Case2  was  10.  All  other  input  parameters  and 
plasma  protein  binding  characteristics  are  identical  and  equal  to  those  for  hydrocortisone. 


The  first  scenario  is  a  rather  straightforward  example  illustrating  how  different 
diffusion  coefficients  (based  on  the  physicochemical  properties  of  the  drug  as  discussed 
in  Chapter  2)  modulate  the  disposition  of  the  drug.  In  Figure  3-1 1,  all  input  parameters 
were  equal  to  those  obtained  in  the  analysis  of  clinical  data  for  hydrocortisone  (20  and 
100  mg  IV,  Table  3-2)  except  for  the  diffusion  coefficients.  The  values  for  K  differ  by  a 
factor  of  10  (1  for  Casel  and  10  for  Case2).  The  smaller  diffusion  coefficient  for  Casel 
caused  a  slow  partitioning  of  drug  into  the  peripheral  compartment  after  IV  bolus 
administration  and,  as  a  result,  drug  concentrations  were  low  (but  rather  constant)  in  the 
peripheral  compartment  while  the  drug  was  rapidly  eliminated  from  the  central 
compartment.  The  concentration  of  drug  in  the  peripheral  compartment  for  Case2, 
however,  rose  rapidly  and  showed  a  rapid  decline  as  well  due  to  the  relatively  free 
partitioning  between  the  central  and  peripheral  compartments.  This  is  a  general  example 
for  drugs  exhibiting  distribution  phases  and  the  same  would  be  predicted  whether  or  not 
nonlinear  distribution  or  elimination  occurs. 
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Figure  3-12.  Simulated  concentration-time  profiles  for  different  values  of  fig. 
These  profiles  were  generated  using  the  gradient-based  distribution  model.  The  unbound 
fraction  in  tissues  for  Casel  was  0.1  and  that  for  Case2  was  0.9.  All  other  input 
parameters  and  plasma  protein  binding  characteristics  are  identical  to  those  for 
hydrocortisone. 


In  Figure  3-12,  a  somewhat  more  complex  scenario  is  presented.  Here,  the  only 
difference  is  found  in  the  fraction  unbound  in  the  tissues  (or  peripheral  compartment) 
with  all  other  input  parameters  equal  to  those  obtained  from  the  analysis  of 
hydrocortisone  data  after  20  and  100  mg  IV  bolus  doses.  For  Casel  fu2  equals  0.1  and  for 
Case2,  0.9.  As  predicted  from  general  pharmacokinetic  principles,  the  greater  degree  of 
tissue  binding  in  Casel  (having  a  lower  unbound  fraction)  in  comparison  to  Case2  results 
in  larger  concentrations  in  the  peripheral  compartment  and  a  larger  volume  of 
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distribution.  Likewise,  this  drug  exhibits  a  longer  terminal  half-life  (observed  in  the 
terminal  slope  of  the  concentration-time  profile  for  the  central  compartment),  which  may 
be  explained  by  the  high  degree  of  tissue  binding  and  the  large  volume  of  distribution. 
Although  the  nonlinear  protein  binding  affects  the  volume  of  distribution  to  some  degree 
(in  that  Vd  for  Casel  rises  and  then  declines  as  was  observed  for  hydrocortisone),  the 
nonlinear  kinetics  are  much  more  evident  in  the  clearance. 


Two-Compartment  Model  with  Gradient-Based  Distribution 
An  Application  of  Fick's  1"  Law  of  Diffusion 

Hepatic  and/or  Renal  Clearance  -  Linear  and/or  Nonlinear  Binding 
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FiRure  3-13.  Interface  for  the  gradient-based  distribution  model  spreadsheet 
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In  this  scenario,  a  change  in  the  tissue  binding  indirectly  affected  the  clearance  of 
the  drug  from  the  central  compartment  by  prolonging  the  time  during  which  there  is 
nonlinear  binding  to  plasma  proteins.  Of  course,  clearance  is  concentration-dependent 
rather  than  time-dependent  and  the  clearance  values  of  both  drugs  (Casel  and  Case2)  are 
equal  for  a  given  concentration.  This  is  evident  from  viewing  the  concentration-time 
profiles  for  the  central  compartment  and  the  clearance  profile  (immediately  below  that  for 
concentration  at  roughly  2.5  h).  For  Case2,  there  is  pronounced  nonlinearity  in  the 
clearance  at  early  time  points,  becoming  linear  at  just  over  4h.  For  Casel,  however,  a 
larger  fraction  of  the  administered  dose  partitions  into  the  tissues,  reentering  the  central 
compartment  over  a  long  period  of  time,  and  clearance  does  not  become  constant  until 
roughly  12  h  after  the  dose  was  given. 
Endogenous  Cortisol 

The  results  of  fitting  endogenous  Cortisol  data  to  the  linear  release  model  with  and 
without  the  introduction  of  a  second  compartment  with  nonlinear  distribution  and 
elimination  are  shown  in  Table  3-3  and  plots  are  contained  in  Figure  3-14. 

Table  3-3.  Parameters  obtained  from  modeling  endogenous  Cortisol  after  placebo. 
Parameter         Linear  Release  Model '      Linear  Release  Model2 

(normal)  (modified^) 

Rmax(ug/h)  3309.2  2288.1 

tmax(h)  20.6  21.9 

tmin(h)  18.5  191 

ke^Ch-1)  0.7  _ 

fut  —  0.7 

K(L/h)  ...  4  j  391 

,       MSC  2.1  2.4 

The  normal  linear  release  model  is  a  one  compartment  model  with  first-order 
elimination.  See  Eq.  3-16  through  Eq.  3-18. 
2The  linear  release  model  was  modified  to  include  gradient-based  distribution 
into  the  periphery  and  nonlinear  hepatic  clearance.  See  Eq.  3-13  and  Eq.  3-15. 
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Despite  having  one  additional  fitted  parameter,  slightly  better  fits  were  obtained 
with  the  gradient-based  distribution  model  compared  to  the  normal  linear  release  model. 
However,  it  may  not  be  necessary  to  use  this  model  considering  that  the  normal 
physiological  concentrations  of  Cortisol  are  much  lower  than  the  concentrations  of 
hydrocortisone  analyzed  in  this  study  after  the  20  and  100  mg  IV  bolus  doses.  In 
addition,  since  exogenous  corticosteroids  suppress  the  release  of  Cortisol,  it  is  predicted 
that  the  concentrations  of  Cortisol  are  sufficiently  low  as  to  render  the  inclusion  of 
nonlinear  processes  an  unnecessary  complication  in  modeling  procedures. 
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Figure  3-14.  Fitted  endogenous  Cortisol  data  after  placebo. 
In  these  fitting  procedures,  average  Cortisol  concentrations  after  placebo  were  used  (n=6). 
In  the  gradient  based  distribution  model,  it  was  assumed  that  Cortisol  was  released  into 
and  eliminated  from  the  central  compartment. 


Conclusions 
Better  model  selection  criteria  (MSC)  values  were  obtained  when  the  two- 
compartment  model  was  redefined  to  allow  nonlinear  distribution  and  elimination. 
Accordingly,  the  nonlinear  model  provides  better  fits  of  the  data,  particularly  at  early 
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time  points.  The  fui  values  for  hydrocortisone  calculated  from  the  total  concentrations 
and  predicted  free  levels  range  from  0.04  to  0.19  for  a  20  mg  dose  and  up  to  0.24  for  a 
100  mg  dose.  The  results  obtained  from  the  preliminary  calculations  and  from  the 
analysis  of  clinical  data  indicate  that  the  nonlinear  plasma  protein  of  hydrocortisone  is 
responsible  for  nonlinear  pharmacokinetics. 

In  regard  to  comparisons  between  the  newly  developed  model  and  the  familiar 
two-compartment  system,  it  should  be  noted  that  the  commonly  used  two-compartment 
body  model  (Eq.  3-6,  Figure  3-5)  provides  a  distribution  phase  only  via  the  rate  constants 
a  and  (3.  The  pre-exponential  factors  A  and  B  are  defined  in  terms  of  Vc  and  no  change 
in  Vd  is  taken  into  consideration  per  se.  However,  in  the  nonlinear  model  with  gradient- 
based  distribution,  Eq.  3-13  and  Eq.  3-14,  the  apparent  volume  of  distribution  changes 
over  time  and  may  be  viewed  as  an  inherent  component  of  the  model  since  fu],  fu2,  Vi, 
and  V2  determine  flux  from  one  compartment  to  the  other.  Previous  models  for 
describing  changes  in  volume  of  distribution  were  based  on  Eq.  3-5,  including  either 
nonlinear  binding  in  plasma  or  tissue  (Gibaldi  and  Perrier  1982).  As  found  in  the 
preliminary  calculations  and  data  analysis  in  the  present  study,  such  models  are  an 
improvement  over  conventional  models  with  constant  values  for  CL  and  Vd  but  fail  to 
describe  the  distribution  phase  after  IV  bolus  administration.  The  manner  in  which 
nonlinear  distribution  is  described  in  the  gradient-based  distribution  model  proposed  in 
this  chapter  is  unique.  Furthermore,  this  model  is  unique  in  its  ability  to  describe  both 
nonlinear  elimination  and  distribution  simultaneously. 

Assuming  that  fui  may  be  calculated  using  affinity  constants  for  drug  binding  to 
plasma  proteins,  this  model  has  the  distinct  advantage  of  describing  nonlinear  processes 
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as  well  as  predicting  overall  tissue  binding  (or  V,/fu2)  and  enzyme  activity,  should  hepatic 
metabolism  be  a  major  elimination  pathway.  In  the  event  that  the  drug  under 
consideration  exhibits  linear  binding  to  plasma  proteins  such  that  ful  is  constant,  the 
integrated  form  of  this  model  may  be  used  and  fui  treated  a  fitted  parameter.  Thus,  the 
new  model  allows  both  distribution  into  the  periphery  and  elimination  to  be  defined  in 
terms  of  physiological  properties  and  established  physical  principles. 


CHAPTER  4 

PHARMACOKINETIC/PHARMACODYNAMIC  MODELING  OF  TOTAL 

LYMPHOCYTES  AND  SELECTED  SUBTYPES  AFTER  BUDESONIDE 

Introduction 

As  stated  in  previous  chapters,  pharmacokinetic/pharmacodynamic  (PK/PD) 
models  enable  the  effects  of  a  drug,  whether  local  or  systemic,  desired  or  undesired,  to  be 
predicted  as  a  function  of  time.  Systemic  effects  monitored  in  previous  clinical  trials 
include  decreased  bone  density  and  growth  (Prummel,  Wiersinga  et  al.  1991;  Doull, 
Freezer  et  al.  1996;  Allen  1998;  Efthimiou  and  Barnes  1998),  blood  eosinophil  counts 
(Pincus,  Humeston  et  al.  1997),  hypertension  and  effects  on  the  cardiovascular  system 
(Sholter  and  Armstrong  2000;  Maxwell,  Moots  et  al.  1994;  Yunis,  Bitar  et  al.  1999), 
psychotic  symptoms  (Curtis,  Fogel  et  al.  1970;  Seifritz,  Hemmeter  et  al.  1994),  and 
cushingoid  appearance  (Wilson,  Blumsohn  et  al.  2000).  Since  many  of  these  effects  are 
often  difficult  to  assess  and  quantify,  the  application  of  PK/PD  models  in  clinical  studies 
often  relies  on  easily  measured  surrogate  markers,  or,  in  broader  terms,  biomarkers  in  the 
blood.  Two  commonly  used  surrogate  markers  for  the  systemic  effects  of  coricosteroids 
are  Cortisol  suppression  (Derendorf,  Hochhaus  et  al.  1993)  and  reduction  of  blood 
lymphocytes  resulting  from  their  redistribution  into  peripheral  tissues  (Fauci  1979). 

In  the  present  study,  the  systemic  effects  of  budesonide  (BUD)  were  monitored 
after  oral  administration.  The  two  formulations  considered  were  Budenofaulk®  (BF)  and 
Entocort®  (EC).  BF  is  a  pH-modified  delayed  release  formulation  that,  in  effect,  achieves 
local  delivery  of  the  drug  after  oral  administration  for  the  treatment  of  inflammatory 
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bowel  diseases  and  EC  is  as  sustained  release  formulation.  The  decline  in  blood 
lymphocytes  was  chosen  as  the  surrogate  marker  for  the  systemic  activity  of  BUD  since 
this  is  directly  related  to  the  suppression  of  the  immune  system  often  observed  in  chronic 
corticosteroid  therapy  and  is,  therefore,  clinically  relevant.  Lymphocytopenia  has  been 
monitored  as  both  a  necessary  outcome  in  the  treatment  of  allergic  inflammation  (Oneda 
1999)  and  an  indication  of  undesired  systemic  activity  in  the  treatment  of  asthma 
(Fokkens,  van  de  Merwe  et  al.  1999)  or  colitis  (Van  Gossum,  Schmit  et  al.  1998).  There 
is  some  indication  that  successful  therapy  of  asthma  results  in  changes  in  the  ratios  of 
helper  and  suppressor  T-lymphocytes  in  both  systemic  circulation  and  bronchoalveolar 
lavage  fluid  (Milgrom  1991)  although  no  direct  correlation  has  been  found  for  T  cell 
downregulation  and  disease  state  indices  (Majori,  Piccoli  et  al.  1997). 

From  Chapter  1 ,  subclasses  of  lymphocytes  play  different  roles  in  regulating  the 
immune  system.  Knowledge  of  how  corticosteroids  influence  the  number  of  cells  in 
circulation  could  allow  immunosuppression  to  be  monitored  and  predicted  more 
effectively.  Pharmacokinetic-phamacodynamic  (PK/PD)  modeling  has  been  used  to  fit 
and  predict  Cortisol  levels  and  total  lymphocyte  counts  following  corticosteroid 
administration  (Derendorf,  Hochhaus  et  al.  1993).  Although  lymphocyte  subsets  have 
been  monitored  in  various  clinical  trials  (Boss,  Neeck  et  al.  1999;  Fokkens,  van  de 
Merwe  et  al.  1999;  Milner,  Kent  et  al.  1999),  previous  PK/PD  studies  of  lymphocyte 
subpopulations  were  limited  to  CD4  and  CD8  after  methylprednisolone  (Milad,  Ludwig 
et  al.  1994)  and  CD3  and  CD8  after  prednisone  (Imani,  Jusko  et  al.  1999).  Since 
lymphocytes  are  modulated  by  both  endogenous  and  exogenous  corticosteroids,  the 
present  PK/PD  study  was  conducted  to  assess  the  combined  effects  of  budesonide  and 
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Cortisol  on  total  lymphocytes  and  all  relevant  subclasses.  Representative  subpopulations 
were  chosen  for  total  T  cells  (CD3),  helper  T  cells  (CD4),  suppressor  T  cells  (CD8),  B 
cells  (CD  19),  and  natural  killer  or  NK  cells  (CD56  and/or  CD  16). 

Clinical  Procedures 

Single  doses  of  budesonide,  Budenofalk  (BF)  and  Entocort  (EC),  and  placebo 
were  administered  at  8  AM  to  the  same  group  of  volunteers.  Five  subjects  received  3  mg 
BF,  four  received  9  mg  BF,  and  five  received  9  mg  EC.  Placebo  capsules  were 
administered  to  six  subjects  in  the  study.  A  commercially  available  radioimmunoassay 
(RIA)  was  used  to  determine  Cortisol  concentrations  and  BUD  concentrations  were 
determined  by  HPLC/RIA  analysis  with  a  sensitivity  of  0.133  ng/mL  (Hochhaus, 
Froehlich  et  al.  1998).  Cortisol  levels  and  lymphocyte  counts  (described  below)  were 
determined  at  the  Medical  Clinic  Bergmannsheil,  University  of  Bochum,  Bochum, 
Germany  and  quantitative  analysis  of  BUD  was  performed  at  the  Department  of 
Pharmaceutics,  University  of  Florida,  Gainesville,  FL. 

The  Blood  samples  were  collected  into  tubes  containing  EDTA  for  total 
lymphocyte  counts  and  containing  heparin  for  subset  analysis.  Total  leukocyte  counts  and 
percentages  were  measured  using  an  automated  cell  counter  (Coulter  STKR). 
Lymphocyte  subsets  were  analyzed  by  flow  cytometry  (FACSCalibur  and  Simulset 
software,  Becton  Dickinson).  Each  blood  sample  had  a  concurrent  sample  containing 
CD54  FITC  and  CD  14  PE  analyzed  to  insure  proper  gating  techniques.  An  isotypic 
control  was  utilized  to  identify  any  non-specific  binding  attributed  to  the  monoclonal 
antibodies  used.  Total  T  cells  were  measured  using  anti-CD3.  The  CD4+  cells  were 
defined  by  dual  staining  with  anti-CD3  and  anti-CD4  and  CD8+  cells  were  defined  using 
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anti-CD3  with  anti-CD8.  The  number  of  B  cells  was  measured  using  anti-CD19.  Natural 
killer  cells  were  CD3 -expressing  CD56/CD16.  The  CD56  and  CD  16  subsets  were 
measured  collectively.  Aliquots  of  blood  were  incubated  with  antibody  pairs  for  dual 
staining  for  15  min  at  room  temperature.  Erythrocytes  were  lysed  and  fixed  using  FACS 
Lysing  solution  and  cells  were  analyzed  within  4  h  after  fixation.  All  monoclonal 
antibodies  were  obtained  from  Becton  Dickinson. 

Pharmacokinetic/Pharmacodynamic  Modeling  of  Blood  Lymphocytes 

Various  aspects  of  the  PK/PD  model  describing  the  combined  effects  of 
endogenous  and  exogenous  corticosteroids  on  lymphocytes  were  discussed  in  preceding 
chapters.  The  complete  model  is  presented  here.  A  modular  assembly  concept  was  used 
to  construct  the  PK/PD  model  and  each  of  the  procedures  is  discussed  below. 
Budesonide  Concentrations 

The  modeling  of  the  kinetic  data  for  the  single  oral  doses  of  3  mg  BF,  9  mg  BF, 
and  9  mg  EC  was  performed  using  a  two-compartment  body  model  with  first-order 
absorption  and  first  order  elimination.  A  lag  time  (Tiag)  was  used  to  describe  the  time 
delays  in  absorption  exhibited  by  the  two  formulations.  This  pharmacokinetic  model  is 
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Eq.  4-1 
where  ka  is  the  absorption  rate  constant,  Tiag  is  the  lag  time  between  the  time  of 
administration  and  the  start  of  GI  absorption,  a  and  p  are  macroconstants  defined  in 
terms  of  the  microconstants  (first-order  rate  constants  for  drug  transfer  between  the 
central  and  peripheral  compartment  and  elimination  from  the  central  compartment),  and 
A  and  B  are  pre-exponential  factors  defined  in  terms  of  dose,  bioavailability,  and  the 
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volume  of  distribution  of  the  central  compartment  (discussed  in  Chapter  3).  Fitted  values 
of  ka  and  (3  were  converted  into  the  absorption  half-life  (ti/2abs)  and  elimination  half-life 
(ti/2eiim),  respectively,  using 

l\l2abs   -      ,  Mud   M/2elim    ~        n       ' 

Individual  concentration-time  profiles  (after  administration  of  a  given  dose  and 
formulation)  for  each  volunteer  were  fitted  to  this  pharmacokinetic  model. 
Pharmacokinetic  parameters  for  a  given  dose/formulation  (3  mg  BF,  9  mg  BF,  or  9  mg 
EC)  were  averaged  prior  to  subsequent  PK/PD  modeling  procedures.  Average  parameter 
values  and  standard  deviations  are  reported.  Although  the  program  used  in  the  fitting 
procedures  contains  a  statistics  option,  averages  and  standard  deviations  of  all 
pharmacokinetic  parameters  were  calculated  manually  from  the  results  obtained  from 
fitting  individual  data  sets  to  Eq.  4-1. 
Circadian  Rhythm  of  Cortisol 

The  circadian  rhythm  of  endogenous  Cortisol  after  placebo  was  described  by  a 
linear  release  model  (Rohatagi,  Bye  et  al.  1996;  Derendorf,  Mollmann  et  al.  1997), 

dCr 


dt 


=  Rc-kr-CCorl  Eq.4-2 
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where  the  release  rate  of  Cortisol  is  described  by  two  linear  functions  (Re)  between  the 
maximum  Cortisol  release  (Rmax)  at  tmax  (6-10  AM)  and  a  minimum  release  at  tmjn  (8  PM- 
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2  AM),  keCort  is  the  elimination  rate  of  Cortisol,  and  Vd  is  the  volume  of  distribution  fixed 
at  a  predetermined  value  of  33.7  L  (Derendorf,  Mollmann  et  al.  1991;  Rohatagi, 
Hochhaus  et  al.  1995).  Placebo  Cortisol  data  was  averaged  (n=6)  and  fitted  to  this  model 
in  order  to  determine  values  of  Rmax,  W,  tmin,  and  k<?on  for  this  particular  group  of 
subjects.  Standard  deviations  were  obtained  from  the  statistics  option  in  the  program  used 
for  the  fitting  procedures  (SCIENTIST,  discussed  below).  The  values  of  these  parameters 
were  fixed  prior  to  subsequent  modeling  procedures. 
Cortisol  Suppression 

The  expression  for  baseline  (placebo)  Cortisol  in  Eq.  4-2  may  be  modified  to 
include  the  suppression  of  Cortisol  by  budesonide.  By  incorporating  an  Emax-based 
expression,  the  Cortisol  release  after  budesonide  administration  may  be  described  by 


dCCo" 
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dt 
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C"  Eq.  4-3 


where  Cf  u    is  the  free  budesonide  concentration  in  plasma,  EC50B^C  is  the  free 
budesonide  concentration  that  reduces  the  release  of  Cortisol  by  50%,  CtotCort  is  total 
Cortisol  concentration,  and  keCort  is  the  first-order  elimination  rate  of  Cortisol  (Rohatagi, 
Bye  et  al.  1996;  Derendorf,  Mollmann  et  al.  1997).  Total  budesonide  concentrations  were 
converted  into  unbound  concentrations  using  a  fraction  unbound  of  0.12  (Ryrfeldt, 
Andersson  et  al.  1982).  Since  most  of  the  parameters  in  this  model  had  been  obtained 
previously  from  the  analysis  of  placebo  data,  EC50B_>C  was  the  only  parameter  obtained  in 
the  PK/PD  analysis  of  treatment  data  (Cortisol  levels  after  the  administration  of  BUD). 
Average  values  and  standard  deviations  are  reported  (obtained  from  the  statistics  option 
in  the  program  used  for  the  modeling  procedures).  Values  of  EC50B^C  for  a  given 
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dose/formulation  (3  mg  BF,  9  mg  BF9,  and  9  mg  EC)  were  determined  and  reported 
separately  in  terms  of  the  free  levels. 
Nonlinear  Binding  to  Plasma  Proteins 

Since  Cortisol  exhibits  nonlinear  binding  to  transcortin  in  the  physiological 
concentration  range  (Ballard  1979b),  a  nonlinear  binding  algorithm  (Meibohm, 
Derendorf  et  al.  1999)  was  used  to  determine  unbound  concentrations.  This  was 
necessary  because  the  pharmacologically  active  unbound  fraction  is  not  constant  but, 
rather,  concentration-dependent.  Cortisol  binds  to  albumin  and  transcortin,  and  total 
Cortisol  concentration  (CtotCort)  can  be  described  according  to 

rCor,         ^Tc     Vtc  •*-/  *Alb  ' V Alb     <-/  rc.orl  p         ,, 

""  i  ,  jf      rCon  \,V       rcon     +  LV  tq.  4-3 

where  KTc  is  the  affinity  constant  between  Cortisol  and  transcortin  (3  x  107  L/mol),  KAib 
the  affinity  constant  between  Cortisol  and  albumin  (5  x  1 04  L/mol),  QTc  the  total 
transcortin  concentration  (0.7  umol/L),  QAib  the  total  albumin  concentration  (550 
umol/L)  (Ballard  1979b),  and  CfCort  the  unbound  Cortisol  concentration.  Since  binding  to 
albumin  in  the  physiological  Cortisol  concentration  range  is  approximately  linear 
(K-Aib-Cf  ort  «  1),  a  simplified  form  may  be  used 


tCort 

c2" =  '7.:,cZn  +ka»-  qm  ■  ccr  +  ccr  Eq.  4-6 
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Solving  for  CfCort  yields  the  unbound  concentration  responsible  for  the  fluctuations  in 
lymphocyte  counts  observed  for  baseline  (placebo)  conditions.  Free  Cortisol  levels  were 
calculated  from  average  Cortisol  data  prior  to  modeling  lymphocyte-time  profiles,  both 
placebo  and  treatment. 
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Lymphocytes  after  Placebo 

An  indirect  response  model  was  used  to  characterize  the  transient  depletion  of  the 
lymphocytes  from  blood  induced  by  the  circadian  rhythm  of  Cortisol  (Wald,  Law  et  al. 
1992;  Mollmann,  Wagner  et  al.  1998).  In  this  model,  the  number  of  lymphocytes  N  in  the 
blood  is  described  by  a  zero-order  influx  of  cells  (rate  constant  k,n)  and  a  first-order 
efflux  of  cells  (rate  constant  koUt), 
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Eq.  4-7 

where  Emax  is  the  maximum  effect  of  Cortisol  (or  corticosteroids  in  general)  and  EC50C^L 
is  the  Cortisol  concentration  producing  50%  of  Emax.  Average  lymphocyte-time  profiles 
(n=6)  for  total  lymphocytes  and  the  CD3,  CD4,  CD8,  CD  19,  and  CD56/16  subtypes  after 
placebo  were  fitted  to  this  model.  As  discussed  below,  placebo  and  all  treatment  data  for 
a  given  lymphocyte  type  were  fitted  simultaneously.  This  was  done  in  order  to  minimize 
the  error  likely  to  occur  when  EC50C^L  and  Emax  are  determined  from  fitting  only  the 
small  effects  of  endogenous  Cortisol  on  lymphocytes. 
Lymphocytes  after  Budesonide 

The  combined  pharmacodynamic  effects  (Ariens  1954;  Meibohm,  Derendorf  et 
al.  1999)  of  budesonide  and  Cortisol  on  lymphocytes  was  described  by 
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Eq.  4-8 


where  Emax  is  the  maximum  effect  of  Cortisol  and  budesonide  on  the  influx  of 


82 

lymphocytes,  CfCortis  the  unbound  Cortisol  concentration,  and  EC50Bh>l  and  EC50C^L  are 
the  unbound  concentrations  of  budesonide  and  Cortisol,  respectively,  which  produce  50% 
of  Emax.  Average  lymphocyte-time  profiles  for  each  lymphocyte  type  after  a  given 
treatment  were  fitted  to  this  model  simultaneously  with  average  placebo  data.  For 
example,  average  total  lymphocyte  data  after  placebo  (n=6),  average  total  lymphocyte 
data  after  3  mg  BF  (n=5),  average  total  lymphocyte  data  after  9  mg  BF  (n=4),  and 
average  total  lymphocyte  data  after  9  mg  EC  (n=5),  four  average  data  sets  in  all,  were 
fitted  simultaneously  to  the  models  in  Eq.  4-7  and  Eq.  4-8  in  order  to  determine  kin,  koUt, 
Emax,  EC50B^L,  and  EC50C^L  for  total  lymphocytes.  Both  EC50B^L  and  EC50C^L  values  are 
reported  as  free  levels.  This  same  procedure  was  used  to  determine  these  parameters  for 
the  selected  subtypes.  Parameter  values  were  reported  along  with  the  standard  deviations 
obtained  from  the  statistics  option  in  the  SCIENTIST  program. 

Multi-Step  Modeling  Procedure 

Absolute  lymphocyte  counts  were  converted  to  a  percentage  of  the  initial  count  at 
8  AM.  All  fitting  procedures  were  performed  using  the  SCIENTIST  software  package 
(SCIENTIST  1993).  Analysis  was  performed  step- wise,  beginning  with  the  kinetic  data 
and  Cortisol  data  after  placebo.  Cortisol  data  after  budesonide  was  then  analyzed  to 
determine  the  degree  to  which  Cortisol  was  suppressed  for  a  particular  dose  and 
formulation.  Due  to  the  small  overall  effects  on  lymphocytes,  lymphocyte  data  after 
budesonide  (3  mg  BF,  9  mg  BF,  and  9  mg  EC)  for  a  given  type  (total  or  subpopulation) 
were  analyzed  simultaneously  with  placebo  data  for  that  type.  Average  values  and 
standard  deviations  for  the  lymphocyte  module  in  the  PK/PD  model  were  obtained  from 
the  statistics  option  in  the  SCIENTIST  program  and  compared  using  Tuckey's  multiple 
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comparison  test  with  a  95%  confidence  interval. 

Noncompartmental  Analysis 
The  area  under  the  curve  (AUC  for  the  concentration-time  or  effect-time  profile) 
was  calculated  for  all  pharmacokinetic  and  pharmacodynamic  data.  For  the 
pharmacokinetic  data,  AUCs  were  calculated  using  the  trapezoidal  rule  for  a  24-hour 
period  (the  duration  of  the  clinical  study)  and  the  terminal  AUC  was  determined  by 
extrapolating  the  concentration-time  profile  to  infinity  using  the  terminal  slope,  obtained 
from  a  semi-logarithmic  plot  (natural  log  of  the  concentration  versus  time)  of  the  data. 
Total  AUCs  were  calculated,  averaged,  and  reported  with  the  standard  deviation  for  each 
dose/formulation  of  BUD.  The  AUCs  of  the  pharmacodynamic  data  (cortisol-time  and 
lymphocyte-time  profiles)  were  determined  for  a  24-hour  period  (the  duration  of  the 
clinical  study)  for  both  placebo  and  treatment  data  using  the  trapezoidal  rule.  No  attempt 
was  made  to  extrapolate  treatment  profiles  back  to  baseline  (placebo)  levels/counts.  The 
AUCs  of  the  dynamic  data  were  then  used  to  compute  the  net  effects  of  BUD  on  both 
Cortisol  and  lymphocytes  over  a  24-hour  period.  These  cumulative  effects  were  calculated 
as  follows: 

The  cumulative  Cortisol  suppression  (CCS)  was  calculated  according  to 

CCS  =  1 00  .aucp>°^-auc<™< 

AUCplacebo  b^-4-9 

where  AUCp|acebo  is  the  24-hour  AUC  for  placebo  Cortisol  data  and  AUCtreat  is  the  24- 
hour  AUC  of  the  cortisol-time  profile  after  budesonide  administration. 

The  cumulative  lymphocyte  reduction  (CLR)  was  calculated  in  a  similar  manner, 

CLR  =  lOQ.AUCp,™bo~AUC'™<.  Fn410 
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where  AUCpiacebo  is  the  24-hour  AUC  for  placebo  lymphocyte  data  and  AUCtreat  is  the  24- 
hour  AUC  of  the  lymphocyte-time  profile  after  budesonide  administration. 

Average  values  for  CCS  and  CLR  are  reported  with  standard  deviations.  Since 
these  parameters  are  based  on  both  placebo  and  treatment  data,  which  differ  in  sample 
size,  the  standard  deviations  (SD)  were  calculated  using  a  weighted  average  of  the  sample 
variances, 


SD  = 


(k.-O-^+C^-D-s' 


Eq.4-11 


O,  +  n2  -  2) 

where  S]  and  S2  are  the  standard  deviations  (the  square  root  of  the  sample  variance)  of 
the  two  quantities  (in  this  case,  the  AUC  of  placebo  data  and  the  AUC  of  treatment  data) 
and  m  and  n2  are  the  respective  sample  sizes,  and  expressed  as  a  percentage  of  the  values 
ofCCSandCLR. 

Results 

Pharmacokinetic  parameters  (lag  times,  rate  constants,  and  pre-exponential 
factors)  obtained  from  fitting  budesonide-time  profiles  to  a  two-compartment  model  in 
macroconstants  are  listed  in  Table  4-1  along  with  the  total  areas  under  the  concentration- 
time  profiles  extrapolated  to  infinity.  Plots  of  the  fitted  concentration-time  profiles 
(average  data  with  standard  deviation)  after  3  mg  BF,  9  mg  BF,  and  9  mg  EC  are 
provided  in  Figure  4-1 . 

From  the  analysis  of  placebo  Cortisol  data  using  Eq.  4-2,  keCort  was  0.61  ±  0.07  h"1 
and  the  AUC24  (the  area  under  the  cortisol-time  profile  over  a  24-hour  period)  was  1521 
±  366  ng/mL/h.  The  maximum  release  of  Cortisol  (Rmax)  was  2853  ±  280  ng/mL/h.  The 
times  of  maximum  and  minimum  Cortisol  release,  tmax  and  tmin,  were  5  AM  ±  1 .5  h  and 
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1:30  AM  ±  1 .4  h,  respectively. 

Values  of  EC50B^C  obtained  from  fitting  Cortisol  data  after  budesonide 
administration  (3  mg  BF,  9  mg  BF,  and  9  mg  EC)  to  the  PK/PD  model  in  Eq.  4-3  are 
provided  in  Table  4-2  along  with  the  areas  under  the  cortisol-time  profiles  (for  a  24-hour 
period)  and  the  cumulative  Cortisol  suppressions  (calculated  according  to  Eq.  4-9).  For 
comparison,  parameters  from  the  analysis  of  placebo  data  are  included  in  Table  4-1.  Plots 
of  the  fitted  data  are  shown  in  Figure  4-2. 

Table  4- 1 .  Pharmacokinetic  parameters  of  budesonide. 
Parameter  3  mg  BF*  9  mg  BF*  9  mg  EC* 

Tlag  (h)  3.47  ±0.55        3.24  ±0.65        0.78  ±  0.25 

A  (ng/mL)  9.15  ±  6.17        5.78±2.16        7.12±  3.95 

B  (ng/mL)  1.98  ±0.74        2.98  ±2.19        2.55  ±1.57 

ti/2abs  (h)  0.12  ±0.05        0.09  ±0.05        0.36  ±0.16 

ti/2eiim  (h)  3.15  ±0.92        2.48  ±0.31         4.56  ±1.17 

a  (h"1)  3.09  ±  2.33        0.89  ±  0.74        2.96  ±  1 .49 

MSC  (average)  4.11  4.93  2.71 

AUCcc  (ng/mL/h)        12.48  ±6.13       16.22  ±4.84       15.49  ±8.19 
*Budenofaulk  (BF),  Entocort  (EC). 

Lag  times  and  absorption/elimination  half-lives  were  obtained  from  individual  fitting  data 
to  a  two-compartment,  Eq.  4-1.  The  area  under  the  curve  (AUC)  of  the  budesonide-time 
profile  was  calculated  using  the  trapezoidal  rule  and  extrapolating  to  infinity.  Standard 
deviations  were  calculated  manually. 

Table  4-2.  Parameters  from  the  analysis  of  Cortisol  data  after  oral  budesonide. 
Parameter  Placebo  3  mg  BF  9  mg  BF  9  mg  EC 

EC50B^C  (ng/mL)  —  0.09  ±  0.04        0.05  ±  0.04        0.08  ±  0.03 

AUC24(ng/mL/h)  1521  ±366  1448  ±167        1313  ±126        1093  ±224 

CCS(%)  —  4.8  ±0.9  13.7  ±2.5  28  2  ±6  3 

Rmax  (ng/mL/h)  2853  ±  280 

tmax(h)*  0500  ±1.5 

tmin(h)*  01 30  ±1.4 

ke^Gf1)  0.61  ±0.07 

*Military  time. 

PK/PD  parameters  of  Cortisol  after  placebo  and  budesonide  are  described  in  Eq.  4-2  and 
Eq.  4-3,  respectively.  EC50B^C  is  expressed  in  terms  of  free  BUD  concentrations.  The 
area  under  the  curve  (AUC)  was  determined  using  the  trapezoidal  rule  for  a  24-hour 
period.  The  cumulative  Cortisol  suppression  (CCS)  was  calculated  according  to  Eq  4-9 
Standard  deviations  for  AUC24  and  CCS  were  calculated  manually  and  those  for  model 
parameters  were  obtained  from  the  statistics  option  in  SCIENTIST. 
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Budesonide  after  3mg  BF 


Budesonide  after  9mg  BF 


Budesonide  after  9mg  EC 


Time  (hr) 


Figure  4-1 .  Fitted  budesonide-time  profiles  for  3  mg  and  9  mg 


BF  and  9  me  EC. 


Two  formulations  of  budesonide,  Budenofalk  (BF)  and  Entocort  (EC),  were 
administered.  Sample  sizes  were  3  mg  BF  (n=5),  9  mg  BF  (n=4),  9  mg  EC  (n=5). 
Pharmacokinetic  data  were  fitted  to  a  two-compartment  model  with  first-order  absorption 
and  elimination,  Eq.  4-1 .  A  lag  time  was  included  to  describe  the  delay  in  the  release  of 
budesonide  from  the  two  formulations. 
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Cortisol  concentration  after  placebo 


Cortisol  concentration  after  3mg  BF 


Time  (hr) 


Cortisol  concentration  after  9mg  BF 


Cortisol  concentration  after  9mg  EC 


Figure  4-2.  Fitted  cortisol-time  profiles  after  placebo  and  oral  budesonide. 
Average  Cortisol  data  after  placebo  and  oral  budesonide  were  fitted  to  the  linear  release 
model,  Eq.  4-2  and  Eq.  4-3,  respectively.  Two  formulations  of  budesonide,  Budenofalk 
(BF)  and  Entocort  (EC),  were  administered.  Sample  sizes  were  3  mg  BF  (n=5),  9  mg  BF 
(n=4),  9  mg  EC  (n=5).  Dotted  lines  in  the  3  mg  BF,  9  mg  BF,  and  9  mg  EC  plots 
represent  the  cortisol-time  profile  after  placebo  (simulated  using  parameters  obtained 
from  fitting  Cortisol  data  to  Eq.  4-2). 
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The  NK  subpopulation  data  (CD56  and/or  CD  16)  did  not  fit  the  present  PK/PD 
model.  All  other  subtypes  and  total  lymphocytes  were  well-described  by  the  model. 
Parameters  obtained  from  fitting  baseline  lymphocyte-time  profiles  (after  placebo)  to  the 
PK/PD  in  Eq.  4-7  are  contained  in  Table  4-3  along  with  the  areas  under  the  lymphocyte- 
time  profiles.  Fitted  data  are  shown  in  Figure  4-3. 

The  EC50  ~*  values  obtained  from  fitting  lymphocyte-time  profiles  after  budesonide 
administration  to  the  PK/PD  model  in  Eq.  4-8  are  shown  in  Table  4-4.  The  fitted  data  for 
total  lymphocytes  and  subtypes  (CD3,  CD4,  CD8,  and  CD  19)  are  shown  in  Figure  4-4 
and  Figure  4-5,  respectively.  A  general  trend  was  observed  in  that  the  EC50B^L  value  for 
the  CD8  subtype  was  larger  and  those  for  CD4  and  CD  19  were  somewhat  smaller. 
However,  due  to  the  large  standard  deviations,  these  values  are  not  statistically  different. 

Results  for  the  areas  under  the  lymphocyte-time  profiles  and  the  cumulative 
lymphocyte  reductions  are  provided  in  Table  4-5. 

Conclusions  and  Discussion 
Budesonide 

As  shown  in  Figure  4-1,  the  BF  formulation  exhibits  a  considerable  lag  time 
before  the  release  of  budesonide,  roughly  3.5  h.  However,  once  absorption  begins, 
concentrations  increase  rapidly.  The  sustained  release  characteristics  of  the  EC 
formulation  are  represented  in  the  slow  rise  in  BUD  concentrations.  The  differences  in 
these  two  formulations  are  also  seen  in  the  fitted  model  parameters  in  Table  4-1 .  The  EC 
formulation  has  a  longer  absorption  half-life  as  well  as  a  longer  elimination  half-life.  It  is 
likely  that  the  sustained  release  of  BUD  influenced  the  elimination  half-life  calculated 
from  the  terminal  slope. 
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Total  after  placebo 


CDS  after  placebo 


CD3  after  placebo 


IB  24 


CD19  after  placebo 
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Figure  4-3.  Fitted  lymphocyte  data  (total  and  subtypes)  after  placebo. 
Average  lymphocyte  data  (n=6)  were  fitted  to  the  PK/PD  model  in  Eq.  4-7  describing  the 
effect  of  the  circadian  rhythm  of  endogenous  Cortisol  on  lymphocyte  counts.  The 
CD56/16  subtype  did  not  fit  the  present  PK/PD  model.  The  data  are  shown  here  and  were 
subjected  only  to  noncompartmental  analysis. 
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Table  4-4.  ECso  ~*L  values  lymphocyte  data  after  budesonide  administration- 
Lymphocytes       EC^L  (ng/mU* 
Total  0.22  ±0.14 

CD3  0.17  ±0.10 

CD4  0.08  ±0.03 

CD8  0.25  ±  0.23 

CD19  0.10  ±0.05 

*EC50  "    values  for  budesonide  acting  on  lymphocytes  were  obtained  from  fitting 
lymphocyte  data  to  Eq.  4-8.  All  data  sets  (placebo,  3  mg  BF,  9  mg  BF,  and  9  mg  EC)  for 
a  given  subtype  were  analyzed  simultaneously.  Other  model  parameters  are  provided  in 
Table  4-3.  Standard  deviations  were  obtained  from  the  statistical  analysis  option  in 
SCIENTIST. 
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Total  lymphocytes  after  3mg  BF 


Total  lymphocytes  after  9mg  BF 
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Figure  4-4.  Fitted  data  for  total  lymphocytes  after  oral  budesonide. 
Average  data  for  total  lymphocyte-time  profiles  after  oral  budesonide  (3  mg  BF,  9  mg 
BF,  and  9  mg  EC)  and  placebo  were  fitted  simultaneously  to  the  model  in  Eq.  4-8.  Dotted 
lines  represent  lymphocyte-time  profiles  after  placebo  (simulated  using  parameters 
obtained  from  fitting  placebo  data  to  Eq.  4-7). 
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Figure  4-5.  Fitted  data  for  the  CD3,  CD4,  CD8.  and  CD  19  lymphocyte  subtypes. 
Average  data  were  fitted  to  Eq.  4-8.  Dotted  lines  represent  the  lymphocyte-time  profile 
after  placebo.  The  CD56/16  subtype  did  not  fit  the  present  PK/PD  model. 
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The  pharmacokinetic  profiles  of  the  different  budesonide  formulations  were  not 
linear  with  respect  to  the  dose  administered.  As  seen  in  Table  4-1 ,  the  area  under  the 
budesonide-time  profile  (AUG*,  extrapolated  to  infinity  using  the  terminal  slope)  was 
somewhat  smaller  for  the  3  mg  dose  of  the  BF  formulation  in  comparison  to  the  other 
doses/formulations.  However,  it  was  not  three  times  smaller  than  that  for  9  mg  BF  as 
might  be  expected  based  on  the  principle  that  the  AUC  is  proportional  to  the  dose  for 
linear  kinetics, 

ATir       DF 

AUC-=~^~  Eq.4-12 

where  D  is  the  dose,  F  is  the  oral  bioavailability,  and  CL  is  total  body  clearance.  This  is 
not  an  indication  that  BUD  exhibits  nonlinear  pharmacominetics  but,  rather,  may  be 
attributed  to  the  small  number  of  volunteers  in  the  study  and  the  typically  high  variability 
in  bioavailability  (F)  after  oral  administration  and,  perhaps,  the  lack  of  sufficient  data 
points  during  the  release  time  of  the  BF  formulation.  The  objective  of  this  study  was  to 
ascertain  relationships  between  corticosteroid  concentrations  and  effects  on  both  total 
lymphocytes  and  selected  subtypes.  This  task  was  not  thwarted  by  the  above  observation. 
Cortisol  after  Placebo 

Baseline  lymphocyte  counts  after  placebo  are  influenced  by  the  circadian 
rhythm  of  endogenous  Cortisol  (Meibohm,  Derendorf  et  al.  1999).  From  viewing  the 
results  in  Table  4-3  of  fitting  baseline  lymphocyte  data  to  the  PK/PD  model  in  Eq.  4-7, 
the  EC50  ~"  values  for  endogenous  Cortisol  acting  on  lymphocytes  are  somewhat  larger 
than  those  previously  reported  (Milad,  Ludwig  et  al.  1994).  In  this  previous  study 
conducted  by  other  researchers,  total  Cortisol  levels  were  used  in  the  modeling  procedures 
(rather  than  the  pharmacologically  active  free  concentrations)  and  the  circadian  rhythm  of 


% 


Cortisol  was  described  in  terms  of  harmonic  oscillations  (a  sum  of  sine  and  cosine 
functions).  This  group  reported  EC50C^L  values  for  total  Cortisol  of  56.0-66.3  ng/mL. 
Converting  these  into  free  levels  using  the  nonlinear  binding  algorithm  applied  in  the 
present  study,  Eq.  4-6,  gives  2.62-3.20  ng/mL,  which  do  not  compare  well  to  the  free 
levels  of  6.2-13.6  ng/mL  found  in  the  present  study.  It  is  unclear  whether  these 
discrepancies  are  due  to  the  performance  of  the  models  or  the  nonlinear  binding  of 
Cortisol.  Comparisons  between  these  results  and  the  referenced  values  are  difficult  since 
it  is  not  known  whether  the  nonlinear  binding  to  plasma  proteins  allows  free  EC50C^L 
values  to  be  accurately  calculated  retrospectively. 
Cortisol  after  Budesonide 

For  the  budesonide  dosages  considered  in  this  study,  Cortisol  suppression  was 
small,  ranging  from  4.8%  to  28.3%.  However,  the  CCS  after  EC  is  considerably  greater 
than  that  after  BF  for  the  same  dose  (9  mg).  This  may  be  explained  on  the  basis  that  the 
sustained  release  of  budesonide  from  the  EC  capsule  suppresses  Cortisol  over  a  longer 
time  period.  As  seen  in  Figure  4-2,  the  Cortisol  profile  after  9  mg  EC  is  still  suppressed 
below  the  baseline  (placebo)  level  at  24  h.  On  the  other  hand,  although  BF  is  a  delayed 
release  formulation,  absorption  of  budesonide  is  rapid  once  the  drug  is  released  and,  due 
to  the  large  clearance  of  BUD,  the  drug  is  rapidly  eliminated  from  systemic  circulation 
(Derendorf,  Hochhaus  et  al.  1998;  Derendorf,  Hochhaus  et  al.  2000).  Cortisol  declines 
roughly  three  hours  after  administration  of  3  mg  and  9  mg  BF,  the  time  at  which 
budesonide  absorption  begins  from  the  pH-modified  delayed  release  capsules.  The  effect 
of  BUD  on  Cortisol  from  the  EC  formulation  is  observed  almost  immediately  after 
administration. 
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Lymphocytes  after  Placebo 

Although  influx  rates  (k,n)  of  B  cells  (CD  19)  in  Table  4-3  are  somewhat  larger 
than  those  for  the  other  subtypes,  efflux  rates  (lw)  are  very  similar  as  are  EmaxL  values. 
The  maximal  effects  of  corticosteroids  (Cortisol  and  BUD)  on  lymphocytes  are  0.8  to  0.9 
or  between  80%  and  90%.  These  values  are  practically  identical  to  the  EmaxL  value  of  0.8 
for  total  lymphocytes  used  or  obtained  in  previous  PK/PD  studies  (Rohatagi,  Hochhaus  et 
al.  1995;  Meibohm,  Derendorf  et  al.  1999).  Thus,  corticosteroids  can  not  completely 
inhibit  the  influx  of  cells  into  the  blood  from  the  extravascular  space.  Likewise,  EC50C"*L 
values  for  the  action  of  Cortisol  on  lymphocytes  are  very  similar  for  total  lymphocytes 
and  all  subtypes  except  CD19.  The  lower  EC50C^L  value  for  CD19  indicates  that  B  cells 
are  somewhat  more  sensitive  to  the  circadian  rhythm  of  Cortisol  than  are  the  other 
lymphocyte  subtypes.  This  is  also  shown  in  Figure  4-3  wherein  the  CD19  subtype  shows 
more  pronounced  diurnal  variations  (inversely  proportional  to  Cortisol,  Figure  4-2)  than 
the  other  subtypes;  Total  lymphocytes  and  CD8  exhibit  only  a  slight  circadian  rhythm. 
Lymphocytes  after  Oral  Budesonide 

The  EC50  ~    values  for  BUD  acting  on  lymphocytes  from  Table  4-4  range 
between  0.10  and  0.25  ng/mL.  These  are  considerably  lower  than  the  EC50  values  of  4.95 
to  6.53  ng/mL  reported  for  methylprednisolone  (Milad,  Ludwig  et  al.  1994)  acting  on 
lymphocytes,  as  might  be  predicted  based  on  the  difference  in  potency.  The  CD8  subtype 
(suppressor  T  cells)  had  the  largest  EC50B"L  value,  indicating  that  this  subtype  is  less 
sensitive  to  the  presence  of  corticosteroids  than  are  the  other  subtypes.  The  larger 
EC50B^L  for  CD8  is  consistent  with  suppressor  T  cells  traits  reported  in  other  studies 
(Milad,  Ludwig  et  al.  1994).  As  seen  in  Table  4-5,  CDS  exhibited  the  smallest  cumulative 
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effects,  as  determined  by  the  CLR.  Negative  CLR  for  some  CDS  data  indicate  that  the 
CDS  count  rose  above  the  placebo  profile  at  some  point  after  BUD  administration.  This 
overshooting  of  the  baseline  value  has  been  previously  documented  (Milad,  Ludwig  et  al. 
1994).  Just  the  opposite  is  observed  for  CD4  (helper  T  cells)  and  CD19  (B  cells),  each 
having  somewhat  smaller  EC50^L  values  and  larger  CLRs  in  comparison  to  the  other 
subtypes.  Thus,  there  is  some  consistency  between  the  PK/PD  analysis  (determination  of 
EC50B-*L)  and  the  purely  noncompartmental  analysis  performed  on  the  raw  data 
(calculation  of  CLR). 

Despite  the  differences  in  lymphocyte  profiles  between  the  BF  and  EC 
formulations,  CLR  values  are  similar  for  a  given  subset  and  no  apparent  dose-response 
relationships  are  observed.  For  example,  from  Table  4-5,  the  CLR  values  for  total 
lymphocytes  after  3  mg  BF,  9  mg  BF,  and  9  mg  EC  are  1 1 .2,  9.9,  and  9. 1 ,  respectively, 
indicating  that  total  lymphocytes  are  reduced  by  roughly  10%  over  a  24-hour  period. 
Such  similarities  were  not  observed  for  CD56/16.  The  small  effects  combined  with  the 
small  study  size  prohibit  definitive  conclusions  to  be  made  concerning  the  effects  of 
budesonide  on  this  subtype. 

Similar  trends  regarding  the  relative  numbers  of  helper  and  suppressor  T  cells  and 
B  cells  have  been  previously  reported,  with  corticosteroids  reducing  CD4  and  CD19 
lymphocytes  more  so  than  CDS  (Cupps,  Edgar  et  al.  1984;  Walzer,  LaBine  et  al.  1984; 
Ozbek,  Yetgin  et  al.  1998).  In  most  other  studies,  comparisons  between  subtypes  were 
based  on  a  limited  number  of  data  points.  The  present  study,  however,  was  more 
comprehensive  in  that  total  lymphocytes  and  all  relevant  subtypes  of  lymphocytes  were 
observed  over  a  24-hour  period  and  compared  to  baseline  (placebo)  counts. 
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Another  important  feature  of  the  present  study  is  that  data  analysis  was 
approached  from  two  different  perspectives.  The  application  of  the  PK/PD  model  allowed 
the  combined  effects  of  BUD  and  Cortisol  on  lymphocytes  to  be  described  over  time  and 
pharmacodynamic  parameters  were  obtained  for  the  effects  of  both  BUD  and  Cortisol. 
These  parameters  may  later  be  used  for  predicting  effects  under  different  conditions  such 
as  variable  doses  and  administration  times.  Of  course,  the  ability  to  describe  and/or 
predict  these  effects  is  dependent  on  the  performance  of  the  model.  For  this  reason,  a 
noncompartmental  approach  (based  only  on  the  area  under  the  lymphocyte-time  profiles) 
was  also  used.  Analyzing  the  data  in  this  manner  allows  unbiased  results  to  be  reported 
without  any  dependency  on  the  construction  or  performance  of  a  PK/PD  model.  A 
correlation  between  these  two  approaches  could  serve  as  a  test  for  the  robustness  of  the 
PK/PD  model.  As  discussed  above,  such  a  correlation  was  observed  between  the 
noncompartmental  approach  and  the  PK/PD  model  used  in  this  study. 

In  general,  the  results  of  the  present  PK/PD  study  suggest  that  total  lymphocyte 
counts  adequately  characterize  the  systemic  effects  of  corticosteroids  and  there  is  no 
apparent  advantage  in  monitoring  specific  lymphocyte  subtypes  solely  for  the  purpose  of 
obtaining  a  suitable  surrogate  marker.  Although  small  overall  effects  on  lymphocytes  are 
observed  after  oral  administration  of  3  mg  and  9  mg  budesonide,  results  from  the  subset 
analysis  suggest  that  corticosteroids  may  have  a  two-fold  impact  on  the  immune  system. 
First,  the  number  of  helper  T  cells  (CD4)  and  B  cells  (CD  19)  in  circulation  is  reduced. 
Thus,  the  ability  to  initiate  an  immune  response  is  impeded.  Second,  the  number  of 
suppressor  T  cells  (CD8)  in  circulation  is  not  reduced  to  the  same  degree  as  the  other 
subpopulations.  With  CD8  cells,  the  subtype  that  suppresses  the  action  of  other  cells  in 
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the  immune  system,  present  in  a  greater  ratio,  sufficient  antibody  production  could  be 
hindered. 


CHAPTER  5 
FACTORS  INFLUENCING  LYMPHOCYTE  REDUCTION 

Introduction 


After  the  administration  of  an  exogenous  corticosteroid  (ECS),  Cortisol  is 
suppressed  and  there  is,  therefore,  less  total  corticosteroid  (endogenous  plus  exogenous) 
acting  on  lymphocytes.  It  may  be  said,  then,  that  Cortisol  suppression  buffers  the  effects 
on  lymphocytes  upon  the  administration  of  an  ECS.  The  results  of  previous  simulations 
using  a  pharmacokinetic/pharmacodynamic  (PK/PD)  describing  the  combined  effects  of 
an  ECS  and  Cortisol  on  blood  lymphocytes  suggested  that  Cortisol  suppression  buffers  the 
effects  observed  on  lymphocytes  and  that  this  buffering  capacity  is  dose-dependent 
(Hochhaus,  Meibohm  et  al.  1998).  This  being  the  case,  there  may  be  an  optimal  dose  for 
which  the  concentrations  of  an  ECS  and  the  decrease  in  endogenous  Cortisol  levels 
counterbalance  one  another.  Theoretically,  such  a  dose  would  produce  no  net  effect  on 
lymphocyte  counts.  In  addition  to  dose  and  drug  properties,  it  has  been  shown  that 
administration  time  and  route  of  administration  also  influence  the  degree  of  lymphocyte 
reduction.  The  dependency  on  the  administration  time  may  be  correlated  to  the  circadian 
rhythm  of  Cortisol  since  baseline  lymphocyte  counts  exhibit  a  circadian  rhythm  inversely 
proportional  to  that  of  endogenous  Cortisol  (Meibohm,  Derendorf  et  al.  1999). 

The  present  study  was  conducted  in  order  to  determine  the  factors  that  influence 
the  lymphocyte  reduction  observed  after  the  systemic  administration  of  corticosteroids. 
Of  major  interest  were  factors  such  as  dose,  administration  time,  dosing  frequency,  the 
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circadian  rhythm  of  Cortisol,  and  corticosteroid  potency.  Pharmacokinetic  properties  were 
investigated  to  a  lesser  degree. 

Methods 

A  PK/PD  approach  was  used  to  assess  the  dependency  of  systemic  effects  on 
factors  such  as  dose,  administration  time,  etc.  The  PK/PD  model  applied  in  this  study 
describes  the  combined  effects  of  Cortisol  and  exogenous  corticosteroids  on  blood 
lymphocytes.  Initial  input  parameters  were  obtained  from  a  clinical  study  of  budesonide 
(Stark,  Werner  et  al.  2001).  This  model  and  the  parameters  obtained  from  the  analysis  of 
clinical  data  are  described  in  Chapter  4.  PK/PD  simulations  were  performed  in  order  to 
generate  lymphocyte-time  profiles.  Systemic  effects  were  quantified  and  compared  using 
the  areas  under  the  lymphocyte-time  profiles  generated  in  the  simulations.  All 
simulations  were  performed  using  the  SCIENTIST  program.  Although  initial  parameters 
were  obtained  from  a  clinical  study  of  single  doses,  all  simulations  were  performed  for 
multiple  doses  at  steady-state. 
Lymphocyte-Time  Profiles 

In  order  to  determine  the  degree  to  which  Cortisol  buffers  the  effects  on 
lymphocytes,  simulations  were  performed  with  and  without  Cortisol  suppression  by  the 
ECS.  Cortisol  suppression  was  removed  from  the  model  simply  by  replacing  the  Cortisol 
concentrations  by  the  baseline  values  observed  after  placebo.  In  order  to  determine  the 
degree  to  which  the  circadian  rhythm  of  Cortisol  influences  the  net  effects  observed  on 
lymphocytes  after  the  administration  of  an  ECS,  some  simulations  were  repeated  without 
circadian  rhythm  in  the  PK/PD  model,  replacing  the  normal  cortisol-time  profile  by  an 
average  Cortisol  level  that  was  suppressed  by  the  ECS  via  an  Emax-type  function. 
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Noncompartmental  Pharmacokinetic/Pharmacodynamic  Parameters 

In  order  to  quantify  the  effect  of  corticosteroids  on  lymphocytes,  a  new  parameter 
was  defined  based  on  the  area  under  the  curve  of  lymphocytes  (AUCL)  after  the 
administration  of  an  ECS  and  placebo.  The  cumulative  lymphocyte  reduction  (CLR)  was 
defined  as  the  percent  difference  between  the  lymphocyte-time  profile  after  the 
administration  of  an  ECS,  AUCL(treatment),  and  that  after  placebo,  AUCL(baseline),  over 
a  24-hour  period.  This  is  shown  in  Eq.  5-1 . 

r,„_.nn  A UCL (baseline)  -  A UCL (treatment) 

LLK  —  1UU Eq.  5-1 

A  UCL  (baseline) 

CLR  is  analogous  to  the  24-hour  cumulative  Cortisol  suppression  (CCS)  calculated  in 
previous  PK/PD  studies  (Derendorf,  Hochhaus  et  al.  1998). 

Since  Cortisol  suppression  may  buffer  the  effects  of  an  ECS  on  lymphocyte 
counts,  another  new  parameter  was  proposed  in  order  to  quantify  this  buffering  effect  of 
Cortisol.  The  Cortisol  buffering  capacity  (CBC)  was  defined  as  the  difference  between  the 
CLR  with  and  without  allowing  Cortisol  to  be  suppressed  by  an  ECS.  The  CBC  is  defined 
as 

CBC  =  CLRW/Ocorlsup  -  CLRw/consup  Eq.  5-2 

where  CLRw/o_coltSup  is  the  cumulative  lymphocyte  reduction  calculated  from  lymphocyte- 
time  profiles  in  which  Cortisol  levels  were  fixed  as  the  baseline  (placebo)  values  and 
CLRw/cortsup  is  the  cumulative  lymphocyte  reduction  calculated  from  lymphocyte-time 
profiles  in  which  Cortisol  levels  are  suppressed  by  an  ECS. 
Input  Parameters 

Initial  parameter  values  were  based  on  those  determined  from  the  PK/PD  analysis 
of  clinical  data  (budesonide  concentrations,  Cortisol  levels,  and  lymphocyte  counts)  after 
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the  administration  of  9  mg  Budenofalk  (BF)  and  placebo  at  8  AM  (Chapter  4).  Although 

Budenofalk  is  a  delayed  release  formulation  of  budesonide  (BUD),  no  lag  times  were 

used  in  these  simulations;  Unless  stated  to  the  contrary,  all  doses  were  treated  as  IV  bolus 

injections  (or  a  completely  bioavailable,  rapidly  absorbed  oral  dose).  The  values  of  the 

parameters  for  the  kinetic  models  are  shown  in  Table  5-1 .  Those  for  the  Cortisol  and 

lymphocyte  models  are  shown  in  Table  5-2. 

Table  5-  1.  Parameters  for  the  kinetic  models- 
Parameter  Value(s) 
D  (me)          0.001,  0.005,  0.04,  0.05,  0.1,  0.5,  1.0*,  5.0 
ke  (hT)  0.25 
Vd(L)  183 
x  (hr)                          24,  12,  constant  infusion 
♦Initial  values  obtained  from  the  analysis  of  clinical  data.  An  IV  bolus  dose  of  1  mg 
produces  concentrations  comparable  to  9  mg  oral  budesonide  with  rapid  absorption. 

Table  5-2.  Parameters  for  the  dynamic  models. 

Parameter  Value(s) 

ECsoecs^c  (ng/mL)       0.00051,  0.00102,  0.0051,  0.0255, 

0.051*,  0.102,  0.51,  2.55,  5.1 
Rmax  (ng/mL/h)  2853 

k^Of1)  0.61 

tmax(h)  5:00  AM 

tmin(h)  1:30  AM 

kin%  (h"1)  88.46 

kout  (h'1)  0.72 

EmaxL  0.75 

EC50C^L  (ng/mL)  15.3* 

ECsoecs^l  (ng/mL)       0.00153,  0.00306,  0.0153,  0.0765, 

0.153*,  0.306,  1.53,7.65,  15.3 
initial  values  obtained  from  the  PK/PD  analysis  of  clinical  data. 

Dose,  Frequency,  and  Disposition 

The  dose  was  varied  over  a  large  range  starting  with  1  mg  IV,  comparable  to  the 
dose  reaching  systemic  circulation  after  9  mg  oral  budesonide.  Administration  time  was 
varied  in  a  systematic  manner  with  the  dose  being  given  once  a  day  at  8  AM,  10  AM, 
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noon,  and  so  on.  To  determine  how  dosing  frequency  effects  the  degree  of  systemic 
activity  of  an  ECS,  the  dosing  interval  was  varied  as  well.  Once-a-day,  twice-a-day,  and 
constant  rate  infusion  regimens  were  considered.  In  order  to  assess  the  influence  of 
administration  time  on  lymphocytes  for  drugs  with  different  kinetic  profiles  (specifically 
different  half-lives,  ti/2),  simulations  were  repeated  for  hypothetical  high-,  middle-,  and 
low-extraction  drugs  while  the  volume  of  distribution  was  held.  These  simulations  were 
repeated  yet  another  time  with  different  volumes  of  distribution  while  total  clearance  was 
held  constant,  producing  the  same  ti/2  as  in  the  simulations  with  various  extraction  ratios. 
Corticosteroid  Potency 

To  determine  how  steroid  potency  modulates  lymphocyte  reduction,  the  EC50 
values  for  an  ECS  acting  on  both  Cortisol  and  lymphocyes  were  varied.  The  EC50  value 
for  Cortisol  acting  on  lymphocytes  was  considered  to  be  a  constant  physiological 
parameter  and  was  not  modified.  From  the  results  of  PK/PD  analysis  of  clinical  data  in 
Chapter  4,  the  initial  ratio  of  EC50ECS~*C:  EC50ECS^L:  EC50C^L  was  approximately  1 :3:300. 
However,  different  PK/PD  studies  (unpublished  results  after  single  and  multiple  doses) 
yielded  different  ratios  due  to  an  increase  in  EC50C^L.  Values  for  EC5oECS^c  and 
EC50       '  were  increased  and  decreased  1 000-fold  from  the  starting  values.  The  results 
on  lymphocyte  reduction  are  reported  and  discussed  in  conjunction  with  the  EC50  ratios. 

Results  and  Discussion 
Dose 

CLR  and  CBC  values  obtained  for  a  large  range  of  doses  of  BUD  (administered  at 
8  AM)  are  shown  in  Figure  5-1.  The  CLR  plot,  analogous  to  the  familiar  dose-response 
curve,  exhibits  a  typical  Emax  profile  due  to  the  manner  by  which  the  influx  of 
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lymphocytes  is  modulated  in  the  PK/PD  model.  The  effect  does  not  approach  100%  with 
large  doses  since  EmaXL  was  75%  and  it  is  not  possible  to  completely  inhibit  the  influx  of 
lymphocytes  from  the  extravascular  space  (Mollmann,  Hochhaus  et  al.  1995;  Rohatagi, 
Hochhaus  et  al.  1995;  Mollmann,  Wagner  et  al.  1998).  From  this  profile  and  that  for  CBC 
(Figure  5-1),  it  is  apparent  that  Cortisol  buffers  the  effects  observed  on  lymphocytes  for 
doses  between  O.OlandlO.O  mg.  Outside  of  this  dose  range,  the  CLR  profiles  with  and 
without  Cortisol  suppression  in  the  model  are  superimposed.  With  the  BUD  parameter  set 
used  in  these  simulations,  the  maximum  buffering  capacity  of  Cortisol  is  found  for  doses 
of  just  less  than  1  mg. 
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Figure  5-1.  Cumulative  lymphocyte  reduction  and  Cortisol  buffering  capacity 
for  a  large  range  of  doses  of  budesonide. 


Administration  Time 


Representative  lymphocyte-time  profiles  (simulated)  for  various  administration 
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times  of  the  same  daily  dose  (lmg  BUD  IV,  roughly  equivalent  to  9  mg  BF  with  rapid 
absorption  and  no  lag  time)  at  steady-state  are  shown  in  Figure  5-2.  In  some  of  the  plots 
(e.g.,  administration  at  4  PM  or  8  PM),  the  lymphocyte  counts  rise  slightly  before 
returning  to  the  baseline  (placebo)  value.  This  is  often  observed  in  clinical  data  and  is  due 
to  the  longer  duration  of  Cortisol  suppression  (Milad,  Ludwig  et  al.  1994).  Since  Cortisol 
is  typically  a  more  sensitive  marker  for  the  presence  of  an  ECS  than  are  blood 
lymphocytes,  the  effects  on  Cortisol  are  observed  over  a  longer  time  period. 
Corresponding  CLR  and  CBC  for  the  same  dosing  schedules  are  shown  in  Figure  5-3. 
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Figure  5-2.  Lymphocyte-time  profiles  after  multiple  doses  of  1  mg  BUD 
at  steady-state  administered  at  various  times. 


The  plots  in  Figure  5-3  illustrate  the  reduction  of  lymphocytes  induced  by  an  ECS 
is  dependent  on  the  administration  time.  Minimal  systemic  effects  are  observed  for  early 
morning  dosing  (4  AM)  while  maximum  effects  are  observed  for  afternoon 
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administration.  Despite  the  small  effects  found  in  these  simulations,  the  systemic  effects 
for  afternoon  administration  were  42%  greater  than  for  morning  administration  of  the 
same  dose.  Hence,  the  buffering  effect  of  Cortisol  suppression  is  more  pronounced  when 
an  ECS  is  administered  in  the  morning;  The  decreased  concentration  of  endogenous 
Cortisol  during  this  time  counterbalances  the  concentration  of  the  ECS  entering  systemic 
circulation. 
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Figure  5-3.  Cumulative  lymphocyte  reduction  ("CLR)  and  Cortisol  buffering  capacity 
(CBC)  observed  after  1  mg  BUD  IV  at  various  administration  times. 


Furthermore,  Cortisol  suppression  has  no  impact  on  lymphocyte  reduction 
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observed  when  an  ECS  is  administered  in  the  afternoon  or  early  evening  as  shown  in  the 
top  graph  in  Figure  5-3.  These  results  correlate  well  to  those  reported  in  previous  PK/PD 
studies  in  which  Cortisol  suppression  was  the  sole  marker  for  the  systemic  activity  of 
ECSs  (Meibohm,  Hochhaus  et  al.  1997).  In  these  studies,  greater  Cortisol  suppression  was 
observed  for  administration  in  the  early  morning  (4  AM).  However,  different  conclusions 
may  be  reached  depending  on  whether  Cortisol  suppression  or  the  reduction  of  blood 
lymphocytes  is  chosen  as  the  surrogate  marker  for  the  systemic  activity  of  ECSs.  The 
buffering  effects  of  Cortisol  are  further  discussed  in  Appendix  B. 

Drug  concentrations,  Cortisol  levels,  and  lymphocyte  counts  for  the  simulations  of 
drugs  having  different  extraction  ratios  are  shown  in  Figure  5-4.  The  CLR  values 
calculated  after  various  administration  times  of  these  hypothetical  ECSs  are  plotted  in 
Figure  5-5.  These  results  indicate  that  although  the  degree  of  systemic  effects  depends  on 
the  administration  time  of  a  typical  ECS  having  a  high  systemic  clearance,  such  is  not  the 
case  for  low  extraction  drugs.  For  low  extraction  drugs,  there  is  greater  accumulation  of 
drug  in  the  body  leading  to  greater  systemic  effects  and,  due  to  the  longer  half-life, 
administration  time  is  no  longer  an  important  factor.  Similar  results  are  observed  for 
drugs  with  long  half-lives  due  to  an  increased  volume  of  distribution  (as  in  Figure  5-5 
with  ti/2  of  7  h).  In  the  latter  case,  effects  were  smaller  since  clearance  was  held  constant 
and  the  drug  did  not  accumulate  appreciably.  The  results  of  these  simulations  may  have 
been  based  on  the  fact  that  during  steady  state  dosing  of  a  drug  with  a  long  half-life,  there 
is  little  fluctuation  between  peak  and  trough  concentrations. 
Dosing  Frequency 

BUD  concentration-,  Cortisol  level-,  and  lymphocyte  count-time  profiles  after 
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Figure  5-5.  The  influence  of  administration  time  on  the  cumulative  lymphocyte 

reduction  (CLR)  observed  after  corticosteroids  with  different  kinetic  profiles. 


different  dosing  regimens  (once  a  day  at  8  AM,  twice  a  day  at  8  AM  and  8  PM,  constant 
infusion)  of  the  same  daily  dose  (1  mg)  of  BUD  are  shown  in  Figure  5-6. 
Noncompartmental  parameters,  CLR  and  CBC,  for  these  dosing  regimens  are  tabulated  in 
Figure  5-7.  CLR  and  CBC  after  various  daily  doses  and  dosing  frequencies  are  shown  in 
Figure  5-8. 

During  a  constant  infusion  of  BUD,  the  lymphocyte-time  profile  is  almost  flat  and 
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no  longer  exhibits  the  fluctuations  induced  by  the  circadian  rhythm  of  Cortisol.  With  more 
frequent  administration  of  the  same  daily  dose,  greater  systemic  effects  were  observed: 
CLR  and  CCS  both  increase.  For  once-a  day  dosing  of  1  mg,  lymphocytes  and  Cortisol 
were  reduced  14%  and  55%  below  baseline  levels,  respectively. 


70 
60 
50 
40 
30 
20 
10 
0 


■  1mg  8am 
D0.5mg  8am, 8pm 
01mg/day  Infusion 


CLR 


CBC 


CCS 


Figure  5-7.  Cumulative  lymphocyte  reduction  (CLR),  Cortisol  buffering  capacity  (CBC), 
and  cumulative  Cortisol  suppression  (CCS)  after  lmg  BUD/day  administered 
with  various  dosing  frequencies. 

When  the  same  dose  was  administered  via  a  constant  rate  infusion,  CLR  and  CCS 
increased  to  22%  and  68%,  respectively,  corresponding  to  a  57%  increase  in  CLR  and  a 
24%  increase  in  CCS  compared  to  the  values  observed  for  once-a-day  dosing.  Increased 
Cortisol  suppression  did  lead  to  a  slightly  larger  buffering  capacity  (CBC)  when  BUD 
was  administered  more  frequently  but  this  increase  was  not  significant  (<1%).  Additional 
concepts  involving  dosing  regimens  and  PK/PD  approaches  for  dose  optimization  are 
discussed  in  Appendix  C. 
Circadian  Rhythm  of  Cortisol 

Since  baseline  lymphocyte  counts  are  modulated  by  the  circadian  rhythm  of 
Cortisol,  simulations  were  performed  with  and  without  including  circadian  rhythm. 
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The  CLR,  CBC,  and  CCS  values  obtained  from  the  simulations  with  and  without  diurnal 
variations  in  Cortisol  are  shown  in  Figure  5-9.  These  results  suggest  that  circadian  rhythm 
does  modulate  the  ability  of  Cortisol  suppression  to  buffer  the  effects  of  an  ECS  on 
lymphocytes.  This  is  most  obvious  in  the  CBC-dose  profiles  wherein  the  CBC  is  still 
increasing  for  a  5  mg  dose.  (It  should  be  noted  that  the  same  results  were  not  observed 
when  simulations  were  conducted  for  drug  administered  via  constant  infusion.  During  a 
constant  infusion,  Cortisol  is  suppressed  and  the  circadian  rhythm  is  not  as  pronounced. 
Thus,  excluding  circadian  rhythm  from  the  PK/PD  model  with  constant  rate  infusion 
input  produced  minimal  changes  in  the  CLR  profiles). 
Potency 

CLR-  and  CBC-dose  profiles  for  hypothetical  ECSs  producing  different  degrees 
of  effects  (various  values  of  EC50ECS"L)  on  lymphocytes  are  shown  in  Figure  5-10.  Similar 
profiles  for  ECSs  with  differing  effects  on  endogenous  Cortisol  are  shown  in  Figure  5-11. 

From  viewing  the  results  in  Figure  5-10,  as  the  EC50ECS"L  increases  such  that  the 
ECS  and  endogenous  Cortisol  are  equipotent  (acting  on  lymphocytes),  the  CLR  decreases, 
indicating  less  overall  systemic  effects.  For  some  doses,  "negative  effects"  are  observed, 
stemming  from  the  longer  duration  of  Cortisol  suppression.  Under  these  conditions,  the 
Cortisol  buffering  capacity  (CBC)  increases  and  maximum  buffering  is  observed  for 
larger  doses.  Note  how  the  ratio  of  EC5oECS"C:EC5oECS"L  increases  from  1:3  to  1 :300  in 
Figure  5-10. 

In  Figure  5-1 1,  the  same  ratios  appear.  In  the  simulations  for  these  CLR/CBC 
profiles,  these  ratios  were  produced  by  increasing  the  EC50ECS"C  values  (i.e.  the  ECS 
exerts  smaller  effects  on  Cortisol)  rather  than  by  decreasing  EC5oECS"L.  Again,  overall 
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Fifiure  5-10.  Cumulative  lymphocyte  reduction  and  Cortisol  buffering  capacity  with 


varying  exogenous  corticosteroid  potency  on  lymphocytes  (ECsn        ) 

♦Ratios  of  EC5oECS"C:EC5oECS"L  are  shown  in  parentheses. 
Starting  value  obtained  from  a  PK/PD  study  of  budesonide. 
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Figure  5-11.  Cumulative  lymphocyte  reduction  and  Cortisol  buffering  capacity  with 
varying  exogenous  corticosteroid  potency  on  Cortisol  (ECg;      c).  " 
♦Ratios  of  EC5oECSC:EC5oECS"L  are  shown  in  parentheses.  " 
T  Starting  value  obtained  from  a  PK/PD  study  of  budesonide. 
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effects  on  lymphocytes  decline,  "negative  effects"  are  observed  for  certain  doses,  and  the 
CBC  values  increase.  However,  under  these  conditions,  maximum  buffering  is  observed 
for  smaller  doses.  Figure  5-10  and  Figure  5-1 1  also  contain  CLR/CBC  profiles  for  which 
the  ratio  of  EC5oECS"c:EC5oECS"L  is  decreased  by  either  decreasing  EC50ECS"L  or  increasing 
EC50     "  ,  respectively.  In  either  case,  when  this  ratio  is  smaller  than  1 :3,  CBC  values 
decline.  Thus,  under  these  conditions,  Cortisol  suppression  has  very  little  buffering  effect 
on  the  reduction  of  lymphocytes  and  would  not  be  relevant  from  a  clinical  viewpoint. 
CLR-  and  CBC-dose  profiles  for  ECSs  with  increasing  and  decreasing  overall 
potency  (such  that  the  initial  ratio  of  EC5oECS"C:EC5oECS"L  values  remained  constant)  are 
shown  in  Figure  5-12.  Obviously,  as  the  potency  of  the  ECS  increases,  greater  effects  are 
observed  on  lymphocytes  (as  seen  in  the  CLR  profiles).  As  expected,  with  larger  EC50 
values,  the  CLR  is  not  as  pronounced.  Providing  the  ratio  of  EC50ECS"c:EC5oECS"L  remains 
constant  (1:3  in  these  simulations),  Cortisol  suppression  buffers  the  effects  observed  on 
lymphocytes  regardless  of  the  ECS  potency.  However,  the  maximum  CBC  is  observed  at 
different  doses.  In  Figure  5-12,  as  potency  increases,  the  maximum  CBC  (roughly  3%  for 
parameter  set)  appears  at  progressively  smaller  doses.  Alternatively,  as  potency 
decreases,  the  maximum  buffering  of  Cortisol  suppression  is  observed  for  larger  doses. 
Simulations  for  EC50  ratios  of  1 :3,  1 :5,  and  1 :20  are  of  particular  interest  as  these  ratios 
have  been  observed  in  the  PK/PD  analysis  of  clinical  data  after  multiple  doses 
(unpublished  results). 

Conclusions 
Although  Cortisol  is  a  sensitive  marker  for  the  presence  of  corticosteroids  in  systemic 
circulation,  it  is  not  necessarily  a  preferred  surrogate  marker  for  systemic  effects.  After 
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Figure  5-12.  Cumulative  lymphocyte  reduction  and  Cortisol  buffering  capacity  with 
varying  overall  exogenous  corticosteroid  potency 
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*Ratio  of  EC50  values  remains  1:3. 
f  Starting  values  obtained  from  a  PK/PD  study  of  budesonide. 
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the  administration  of  an  ECS,  endogenous  Cortisol  levels  decline.  It  may  be  said,  then, 
that  the  Cortisol  in  systemic  circulation  is  being  replaced  by  another  corticosteroid,  albeit 
a  more  potent  one  in  most  circumstances.  However^both  the  present  lymphocyte  study 
and  previous  Cortisol  studies  revealed  that  minimal  systemic  effects  areobserved  for 
once-a-day  dosing  in  chronic  corticosteroid  therapy  and  that  administration  time  is  an 
important  factor  (Meibohm,  Hochhaus  et  al.  1997).  When  Cortisol  is  monitored  as  a 


surrogate  markerfor  the  systemic  effects^ojfanECS^maller  effects  arepredicted  for 
afternoon  administration.  Just  the  opposite  is  true  for  the  present  study  whej 


lymphocyte  counts  were  chosen  a  surrogate  marker.  The  results  of  this  study  indicate  that 


smaller  systemic  effects  will  be  achieved  when  an  ECS  is  administered  in  the  morning. 

In  comparing  the  results  of  PK/PD  modeling  of  various  corticosteroids, 
particularly  those  used  in  the  treatment  of  asthma,  it  was  noted  that  the  ratio  of  EC50 
values  for  the  ECS  on  Cortisol  and  on  lymphocytes  was  often  1 :3  (Rohdewald,  Rehder  et 
al.  1987;  Mollmann,  Hochhaus  et  al.  1995;  Rohatagi,  Hochhaus  et  al.  1995;  Mollmann, 
Derendorf  et  al.  1997;  Meibohm,  Derendorf  et  al.  1999).  This  was  true  despite  the  wide 
range  in  potencies  of  these  drugs  and  no  explanation  has  been  proposed  for  this  finding 
thus  far.  All  that  may  be  concluded  is  that  Cortisol  is  normally  three  times  more  sensitive 
to  the  presence  of  an  ECS  than  are  lymphocytes.  However,  analysis  of  some  clinical  data 
revealed  different  EC50  values  with  EC50ECS"c:EC5oECS"L  ratios  between  1 :3  and  1 :20.  The 
simulations  in  the  present  study  were  conducted  over  a  wide  range  of  potencies  to 
account  for  these  differences  and  enable  prediction  of  the  effects  on  lymphocytes  should 
the  EC50  values  vary  from  the  1 :3  ratio  observed  in  Chapter  4.  The  different  ratios  from 
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other  studies  were  due  to  changes  in  EC50ECS"C  For  multiple  dosing,  an  ECS  may  have 
diminished  effects  on  Cortisol. 

For  drugs  that  differ  only  in  their  receptor  binding  affinities,  the  lymphocyte-time 
profiles  for  various  doses  would  be  predicted  to  be  superimposable  when  the  dose  range 
is  shifted.  This  was  observed  in  the  present  study  as  illustrated  in  the  CLR/CBC  profiles 
in  Figures  5-12.  It  should  be  noted  that  the  mechanism  is  not  well  understood  for  either 
the  ECS-cortisol  interaction  or  the  ECS-lymphocyte  interaction.  However,  these  appear 
to  be  receptor  mediated  and  Emax-based  functions  adequately  describe  the  decrease  in 
both  the  release  of  Cortisol  and  the  influx  of  lymphocytes  observed  after  the 
administration  of  ECSs.  As  EC5oECS"C:EC50ECS"L  increases  above  1:3,  the  buffering 
capacity  of  Cortisol  plays  a  larger  role  in  modulating  the  reduction  of  blood  lymphocytes. 
If  this  increase  is  due  to  an  increase  in  EC50ECS"L,  the  maximum  buffering  is  found  at 
larger  doses.  The  opposite  is  observed  if  the  increase  in  the  EC50  ratio  is  due  to  a  decrease 
in  EC5o     "  .  The  buffering  capacity  of  Cortisol  has  virtually  no  influence  on  lymphocyte 
reduction  when  this  ratio  is  smaller  than  1:3.  It  may  be  concluded,  therefore,  that  if  the 
ratio  of  EC50     "  :EC50      '  is  1 :3  or  greater,  Cortisol  suppression  buffers  the  systemic 
effects  of  ECS  as  determined  by  lymphocyte  counts.  This  buffering  capacity  is  of  greater 
clinical  relevance  as  this  ratio  increases.  For  more  potent  ECSs,  the  maximum  buffering 
capacity  will  be  seen  for  smaller  doses. 

The  use  of  noncompartmental  parameters  in  this  study  facilitated  the 
quantification  and  comparison  of  the  simulated  lymphocyte-time  profiles.  These 
parameters  were  based  on  the  area  under  the  lymphocyte-time  profile  (treatment,  placebo, 
and  with  Cortisol  fixed  at  the  baseline  levels)  and  allowed  the  effects  on  lymphocytes  to 
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be  described  over  a  given  time  period  with  a  single  value  (the  cumulative  lymphocyte 
reduction,  CLR,  or  the  Cortisol  buffering  capacity,  CBC).  This  was  an  important 
consideration  due  to  the  number  of  simulations  conducted  and  the  number  of  input 
parameters  in  the  PK/PD  model. 

Other  noncompartmental  parameters  were  proposed  during  the  course  of  this 
study.  However,  due  to  the  complexity  of  the  PK/PD  model  describing  the  combined 
effects  of  Cortisol  and  an  ECS  on  blood  lymphocytes,  these  parameters  were  not  readily 
applicable  to  actual  or  simulated  lymphocyte  data.  The  parameters  quantify  both  the 
duration  (the  PK/PD  equivalent  to  mean  residence  time)  and  the  magnitude  of  a  given 
effect  and  are  discussed  in  Appendix  D. 


CHAPTER  6 
CONCLUSIONS 


In  the  present  study,  PK/PD  modeling  of  the  systemic  effects  observed  after  the 
administration  of  corticosteroids  was  reviewed  from  various  perspectives.  The  first 
chapter  dealt  with  basic  principles  of  PK/PD  modeling  and  the  application  of  surrogate  or 
biomarkers.  Although  Cortisol  has  been  described  as  a  surrogate  marker  for  the  systemic 
activity  of  corticosteroids  in  previous  studies,  it  may  better  be  termed  a  biomarker.  In 
order  to  be  considered  a  surrogate  marker,  a  direct  correlation  must  be  shown  between 
the  effect  on  the  marker  and  the  observed  clinical  outcome.  Although  Cortisol  is  a  very 
sensitive  marker  for  the  presence  of  an  exogenous  corticosteroid  in  systemic  circulation, 
it  is  not  necessarily  a  marker  for  systemic  effects.  Upon  the  administration  of  an 
exogenous  corticosteroid,  Cortisol  is  suppressed,  often  dramatically.  However,  in  this 
scenario,  Cortisol  is  being  replaced  by  another  corticosteroid,  albeit  a  more  potent  one  in 
most  cases.  Thus,  although  the  exogenous  drug  administered  may  produce  systemic 
effects,  such  effects  may  not  be  attributed  directly  to  suppressed  Cortisol  levels. 

The  other  surrogate  or  biomarker  for  the  systemic  activity  of  corticosteroids 
considered  in  this  study  is  the  reduction  of  blood  lymphocytes.  This  surrogate  is  more 
closely  linked  to  the  clinical  effects  observed  during  chronic  corticosteroid  therapy, 
namely  immunosuppression.  As  discussed  in  the  review  of  the  literature  and  illustrated  in 
the  analysis  of  clinical  data,  corticosteroids  reduce  the  number  of  lymphocytes  in 
circulation.  This  reduction  occurs  rapidly  and  it  has  been  hypothesized  that  it  is  due  to  a 
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decline  in  the  influx  of  cells  into  the  blood.  As  discussed,  the  reduction  in  lymphocytes  is 
directly  related  to  the  suppression  of  the  immune  system.  Several  lymphocyte  subtypes, 
components  of  the  immune  system  that  play  different  roles  in  generating  as  immune 
response,  were  considered  in  this  study. 

In  Chapter  2,  various  binding  properties  of  corticosteroids  were  considered.  This 
was  done  for  two  main  reasons.  The  first  is  that  since  only  unbound  drug  is 
pharmacologically  active,  knowledge  of  plasma  protein  binding  is  critical  for  describing 
the  effects  of  a  drug  using  PK/PD  models.  Although  protein  binding  may  not  be 
considered  a  pharmacodynamic  parameter,  it  certainly  influences  the  degree  of  the 
response  observed  for  a  given  dose  of  a  given  corticosteroid.  The  second  reason  for 
studying  binding  properties  is  that  these  strongly  impact  the  pharmacokinetic  disposition 
of  a  drug.  In  general,  protein  binding  may  control  both  clearance  and  the  volume  of 
distribution.  Most  new  therapeutic  corticosteroids  are  characterized  by  a  high  systemic 
clearance,  and,  thus,  protein  binding  has  minimal  effects  on  the  rate  of  metabolism. 
However,  these  drugs  differ  in  their  volumes  of  distribution,  determined  by  the  relative 
degree  of  binding  to  plasma  proteins  and  to  tissues.  Both  of  these  binding  phenomena 
were  studied  and  correlated  to  the  physicochemical  properties  of  the  drug. 

The  degree  of  binding  to  plasma  proteins  and  tissues  was  again  addressed  in 
Chapter  3,  specifically  in  regard  to  the  nonlinear  binding  to  plasma  proteins  exhibited  by 
hydrocortisone.  Both  exogenous  hydrocortisone  and  endogenous  Cortisol  were  studied  to 
assess  the  degree  to  which  nonlinear  binding  affects  the  pharmacokinetic  disposition. 
Unlike  many  other  corticosteroids  currently  used,  hydrocortisone  is  a  low  extraction  drug 
with  a  clearance  considerably  lower  than  liver  blood  flow.  It  was  found  that  the  nonlinear 
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binding  influences  both  elimination  and  distribution  and  a  new  pharmacokinetic  model 
was  developed  to  describe  these  nonlinear  processes  simultaneously. 

Chapters  4  and  5  are  linked  by  the  application  of  the  PK/PD  model  describing  the 
combined  effects  of  exogenous  corticosteroids  and  endogenous  Cortisol  on  total 
lymphocytes  and  selected  subtypes  in  circulation.  Clinical  data  (kinetic  and  dynamic) 
were  fitted  to  this  model.  The  data  were  also  analyzed  using  noncompartmental  methods 
based  on  the  area  under  the  effect-time  profile.  It  was  found  that  the  degree  of  systemic 
effects  is  small  upon  single  doses  of  9  mg  oral  budesonide,  with  total  lymphocytes  being 
reduced  by  roughly  10%  over  a  24-hour  period.  The  analysis  of  lymphocyte  subtypes 
revealed  that  both  B  cells  and  helper  T  cells,  the  cells  initiating  an  immune  response,  are 
reduced  to  a  greater  degree  (up  to  20%)  while  suppressor  T  cells,  the  subtype  that 
terminates  the  immune  response,  are  not  appreciably  reduced  during  the  same  time 
interval.  The  EC50  values  obtained  from  the  PK/PD  fitting  procedures  correlated  well 
with  the  24-hour  cumulative  lymphocyte  reduction.  The  subtypes  exhibiting  the  greatest 
decline  (B  cells  and  helper  T  cells)  also  had  the  smallest  EC50  values,  indicating  that 
these  subtypes  are  more  sensitive  to  the  presence  of  an  exogenous  corticosteroid  than  are 
other  subtypes.  Likewise,  the  subtype  exhibiting  the  least  decline  (suppressor  T  cells)  had 
the  largest  EC50  value. 

The  parameters  obtained  in  the  PK/PD  fitting  procedures  were  then  used  for 
numerous  simulations  for  predicting  the  systemic  effects  of  corticosteroids  under 
different  conditions,  one  of  the  key  features  in  PK/PD  modeling  and  one  that  is  gaining 
more  exposure  in  drug  development.  The  complex  interplay  of  endogenous  and 
exogenous  corticosteroids  was  studied  and  it  was  found  that  Cortisol  suppression  buffers 
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the  effects  of  exogenous  corticosteroids  on  lymphocytes.  One  of  the  most  clinically 
relevant  findings  is  that  the  degree  of  systemic  effects  observed  during  corticosteroid 
therapy  is  dependent  on  the  administration  time.  Additional  simulations  revealed  that  this 
time-dependency  is  mainly  due  to  the  underlying  circadian  rhythm  of  Cortisol, 

Throughout  the  present  study,  other  concepts  were  introduced  for  designing 
dosing  regimens  based  on  both  pharmacokinetic  and  pharmacodynamic  parameters.  With 
the  large  variability  in  pharmacological  response  often  observed  during  drug  therapy,  it 
was  believed  that  dosing  regimens  could  be  better  customized  by  considering  both  kinetic 
and  dynamic  factors.  Due  to  the  complex  interaction  between  exogenous  and  endogenous 
corticosteroids  on  blood  lymphocytes,  many  of  the  ideas  proposed  can  not  readily  be 
applied  when  lymphocyte  counts  are  used  as  a  surrogate  marker  for  the  systemic  effects 
of  corticosteroids.  Less  complex  models  were  considered  and  PK/PD  approaches  to 
designing  dosing  regimens  are  presented  in  Appendix  C. 

Novel  noncompartmental  methods  for  PK/PD  data  analysis  for  quantifying  both 
the  duration  and  the  magnitude  of  the  effects  observed  after  drug  administration  are 
presented  in  Appendix  D.  New  parameters  were  defined  mathematically  and,  when 
applicable,  derived  from  basic  concepts.  These  new  parameters  allow  data  to  be  analyzed 
without  the  assumption  of  a  specific  pharmacokinetic  or  pharmacodynamic  model. 
Application  of  these  methods  to  the  study  of  corticosteroids  on  blood  lymphocytes  was 
hindered  by  the  diurnal  baseline  values  of  lymphocytes.  The  theoretical  framework  for 
statistical  moments  in  PK/PD  analysis  was  developed  and  applied  with  less  complex 
effect  models.  In  addition,  special  conditions  were  proposed  to  allow  these  parameters  to 
be  determined  in  cases  where  the  required  integrals  do  not  exist  in  closed  form. 
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The  primary  motive  for  the  continuance  of  pharmacokinetic/pharmacodynamic 
(PK/PD)  development  is  to  construct  meaningful  physiological  models  that  may 
ultimately  enable  a  clinician  to  determine  the  right  dose  for  the  right  patient  and 
administer  it  in  such  a  way  that  the  desired  therapeutic  effects  are  maintained  while  all 
others  are  minimized.  It  is  hoped  that  the  ideas  conceived,  the  understanding  gained,  and 
the  conclusions  drawn  before,  during,  and  after  this  investigation,  respectively,  supported 
by  a  sound  theoretical  framework,  sufficient  experimental  data,  and  meticulous 
examination  of  the  facts,  will  sustain  this  mission. 


APPENDIX  A 
A  TWO-COMPARTMENT  MODEL  WITH  GRADIENT-BASED  DISTRIBUTION 


Most  drugs  exhibit  a  distribution  phase  after  an  intravenous  (IV)  bolus  dose. 
Thus,  two-compartment  systems  and  bi-exponential  curves  may  be  necessary  to 
adequately  describe  the  pharmacokinetics  of  most  drugs.  If  distribution  occurs  by  passive 
diffusion,  the  concentration  gradient  between  the  central  and  peripheral  compartments  is 
the  driving  force  for  drug  distribution.  Surprisingly,  this  concept  is  not  applied  in  the 
traditional  construction  of  the  two-compartment  body  model  in  pharmacokinetics  and, 
rather,  drug  distribution  is  described  using  two  independent  first-order  rate  constants.  The 
two-compartment  model  presented  here  is  based  on  passive  diffusion  following  the 
concentration  gradient  between  the  two  compartments.  With  the  use  of  free 
concentrations,  this  model  is  unique  in  its  ability  to  describe  changes  in  the  rate  of 
distribution  due  to  nonlinear  binding  to  either  plasma  proteins  or  peripheral  tissues. 

Fick's  first  law  of  diffusion  states  that  the  net  flux  of  matter  across  a  semi- 
permeable membrane  is  determined  by  the  concentration  gradient  on  either  side  of  the 
membrane  and  a  proportionality  constant.  Since  precipitated  solutes  or  bound  material  is 
not  part  of  this  gradient,  only  free  particles  should  be  considered.  This  is  illustrated  in 
Figure  A-l  and  the  net  flux  across  the  membrane  is  given  by 

dA 


dt 


--K-VC/  Eq.  A-l 
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where  A  is  the  amount  of  material,  K  is  the  diffusion  coefficient  (a  proportionality 
constant  with  units  of  volume/time),  and  VCf  is  the  gradient  in  free  concentration. 


Compartment  1 

O      o 

Dissolved/Unbound 

o 


Compartment  2 

o 

Dissolved/Unbound 

o 


Figure  A-l .  Diffusion  of  free  matter  across  a  semi-permeable  membrane 


In  the  context  of  pharmacokinetics,  elimination  of  the  drug  often  occurs  from  the  central 
compartment  (compartment  1)  by  hepatic  metabolism  or  renal  processes  as  illustrated  in 
Figure  A-2. 


Compartment 


Compartment  2 


K 


Passive  Diffusion 


(  1st  order  elimination  rate  constant) 


Figure  A-2.  A  two-compartment  pharmacokinetic  model  with  passive  diffusion 
and  first-order  elimination. 


The  differential  equation  describing  the  change  in  amount  of  drug  in  the  central 
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compartment  (Ai)  is  then 
dA 


L\ 


,   =-K-(crf^-C2-fu2)~k]0-A]  Eq.  A-2 

where  k)0  is  the  first  order  elimination  rate  constant  from  the  central  compartment  (with 
units  of  time" )  and  Ci  and  C2  are  the  total  concentrations  in  central  and  peripheral 
compartments,  respectively,  which  have  been  converted  into  free  levels  by  multiplying 
by  the  respective  unbound  fractions,  ful  and  fu2.  Dividing  this  equation  through  by  Vu  the 
volume  of  the  central  compartment,  gives  the  rate  expression  for  the  concentration  in  the 
central  compartment, 

dC.         K-(C,-f.-C,-f,) 

L  = v    l    Ju\      ^-2    Jul)      h-       r  n       a    i 

dt  v  *io-ci-  Eq.A-3 

The  change  in  amount  of  drug  in  the  peripheral  compartment  my  be  expressed  as 

dA 

—[^  =  -K-(C2-fu2-C]-fu])  Eq.  A-4 

where  K  is  the  same  diffusion  coefficient  as  given  in  Eq.  A-2  since  the  membrane  is 
assumed  to  be  equally  permeable  from  both  sides.  Dividing  Eq.  A-4  through  by  V2,  the 
volume  of  the  peripheral  compartment,  gives  the  rate  expression  in  terms  of 
concentration, 

dC2_    K-(C2-fu2-CrfJ 

dt  "  v2  •  Eq-A'5 

After  collecting  constants,  Eq.  A-3  and  Eq.  A-5  may  be  rewritten  as 


dCx 

~d7 

and 


^^  +  k   ) 
V  10 


Ci  +     77^ '  Q  Eq.  A-6 
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dC2      K-f.    „      K-f, 


dt 


l\ 


Eq.  A-7 


Eq.  A-6  and  Eq.  A-7  comprise  a  system  of  separable  first-order  linear  differential 
equations  that  must  be  solved  in  order  to  obtain  integrated  functions  for  Ci  and  C2. 

Laplace  transforms  may  be  used  to  transform  the  rate  expressions  from  Eq.  A-6 
and  Eq.  A-7  in  the  time  domain  to  a  system  of  linear  differential  equations  in  the  Laplace 
domain, 


s-C,-C,(0)  =  - 


X'ful 


+  k 


10 


C,+ 


K-ful 


V, 


c, 


s-C2-C2(0)  = 


K'fu\        Q  K-  fu2        Q 

v'1      v''2 


Eq.  A-8 


Eq.  A-9 


where  s  is  the  Laplace  operator  and  Ci(0)  and  C2(0)  are  the  initial  concentrations  of  drug 
in  the  central  and  peripheral  compartments,  respectively.  After  an  IV  bolus  dose,  the 
initial  concentration  in  the  peripheral  compartment  is  zero  while  that  in  the  central 
compartment  is  dose  (D)  divided  by  volume, 

D 


C,(0)  = 


V, 


Eq.  A- 10 


The  expression  for  Q  is  therefore  nonhomogeneous  (can  not  be  set  equal  to  zero).  After 
collecting  like  terms,  this  system  of  linear  equations  may  be  written  as 


s  +  £10  + 


K-f* 


C, 


i     j 


■C2=C,(0) 


Eq.A-11 


and 


V2 


■c,+ 


s  + 


K-fu, 


C2  =0. 


Eq.  A-12 
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The  method  of  determinants  (Cramer's  Rule)  may  be  used  to  solve  the  equations  for  each 
compartment. 


C 


C,(0) 


K-/.j 


'.+  *•'■' 


v 


V. 


2        J 


s  +  kl0  + 


*•/., 


V, 


i     J 


*■/., 


s  + 


*-/„2 


2        J 


Eq.A-13 
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i     y 
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v    ^     y 


v     P.     y 
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2     y 
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The  denominator  in  Eq.  A- 14  may  be  factored  such  that 

(s  +  a)(s  +  ^)  =  s2  +(a  +  /?)-s  +  cr/? 
and  rewritten  as 


Eq.  A- 14 


C,(0) 


'.+*•'■' 


C, 


I' 


2     y 


(s  +  a)(s  +  /?) 


Eq.A-15 


or 


C, 


P(s) 

e(s) 


Eq.A-16 


The  roots  9?  of  the  polynomial  Q(s)  are  (-a)  and  (-p).  The  Heaviside  Expansion  Theorem 
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states 

In  the  present  example,  the  polynomial  Q(s)  is 
and  the  first  derivative  is  then 


Q'(s)  =  2-s  +  a  +  0. 


Thus, 


c,(O  =  c,(0) 


K'f«2     J  (K-fu2 


a 


J  e-'  + 


V. 


{p-a) 


-fi 

(a-0) 


-pi 


Eq.A-17 


or,  in  simpler  notation, 


C,{t)  =  A-e'a'  +B-e-p' 


where 


D 


fK-fu2        ^ 
-a 


v    V2 


J 


Vr{fi-a) 


and 


D 


B  = 


K-fu2 

V2 


-P 


Vr(a-fi) 


The  values  of  a  and  P  from  Eq.  A- 15  may  be  determined  from  the  relationships 
(i)         a  +  /?  =  £10+^^  +  ^LAl 


Vy  V2 


and 
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(ii)        a-p  = 


Rearranging  (i)  gives 


P  =  K  + 


*■/„  ,  K-fu7 


K  v2 


-a . 


Substituting  this  expression  for  p  into  (ii)  and  rearranging  gives 


a   -a 


k]0  + 


£■/„>   ,  K-f* 


V,  V2     , 


+ 


K-k]0-ful 


=  0. 


After  solving  this  quadratic  equation,  a  and  (3  may  be  written  as 


a  =  -■{ 
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The  integrated  form  of  the  free  concentration  gradient-based  distribution  model  is 
applicable  only  to  linear  pharmacokinetics.  If,  for  example,  fu]  is  not  constant,  the 
differential  form  must  be  used. 


APPENDIX  B 
BUFFERING  EFFECT  OF  CORTISOL  SUPPRESSION  ON  THE  REDUCTION  OF 

BLOOD  LYMPHOCYTES 

It  has  been  hypothesized  that  Cortisol  suppression  upon  the  administration  of 
certain  doses  of  an  exogenous  corticosteroid  "buffers"  the  systemic  effects  as  determined 
by  the  decline  of  blood  lymphocytes  (Hochhaus,  Meibohm  et  al.  1998).  This  hypothesis 
is  based  on  the  fact  that  lymphocyte  counts  are  influenced  by  both  endogenous  and 
exogenous  corticosteroids.  If  the  exogenous  drug  administered  suppresses  Cortisol,  there 
is  less  total  corticosteroid  present.  The  decline  in  lymphocyte  counts  is  then  less 
pronounced.  Theoretically,  it  may  be  possible  to  determine  a  dose  for  which  no  net 
effects  on  lymphocytes  are  observed. 

This  buffering  phenomenon  is  complex,  as  it  is  dependent  on  the  drug-cortisol, 
drug-lymphocyte,  and  cortisol-lymphocyte  interactions.  In  a  series  of  PK/PD  simulations 
of  lymphocyte -time  profiles,  it  was  found  that  least  effects  on  lymphocytes  are  observed 
for  the  4  AM  administration  of  a  1  mg  dose  of  budesonide  while  greatest  effects  for  the 
same  dose  and  same  route  of  administration  are  observed  upon  4  PM  administration. 
Lymphocytes  exhibit  a  circadian  rhythm  influenced  by  and  inversely  proportional  to  that 
of  Cortisol  (Meibohm,  Derendorf  et  al.  1999).  Thus,  it  was  predicted  that  the  underlying 
fluctuations  in  Cortisol  were  responsible  for  the  distinct  effects  on  lymphocytes  observed 
for  different  administration  times.  An  additional  series  of  simulations  and  calculations 
proved  this  to  be  the  case.Concentration-time  profiles  were  simulated  for  budesonide  and 
Cortisol  after  multiple  IV  bolus  doses  at  steady  state.  A  0. 1  mg  dose  was  used  rather  than 
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1.0  mg  for  scaling  purposes  (discussed  later).  Doses  were  administered  at  either  4  AM  or 
4  PM.  The  total  concentration-time  profiles  were  then  converted  to  free  concentration- 
time  profiles  using  a  free  fraction  of  budesonide  of  0.12  (Ryrfeldt,  Andersson  et  al.  1982) 
and  a  nonlinear  binding  algorithm  for  Cortisol  (Meibohm,  Derendorf  et  al.  1999).  This 
was  done  since  any  effects  observed  on  lymphocytes  are  due  to  unbound, 
pharmacologically  active  molecules.  Since  budesonide  is  more  potent  than  Cortisol, 
converting  free  Cortisol  into  free  budesonide  equivalents  normalized  the  effective 
concentrations  of  Cortisol  and  budesonide.  This  was  done  by  multiplying  the  Cortisol 
concentrations  by  the  ratio  of  the  EC50  value  for  the  budesonide-lymphocyte  interaction 
to  the  EC50  value  for  the  cortisol-lymphocyte  interaction, 


cjr(^/v)=cjr-3b-  Eq.B-i 


'50 


This  conversion  could  have  also  been  performed  using  the  ratio  of  relative  receptor 
binding  affinities  since  the  ratio  of  binding  affinities  is  inversely  proportional  to  that  of 
the  EC50  values, 

R  P  A  Cnr' 

C%(equiv)  =  C%.f^.  Eq.B-2 

The  free  budesonide  concentrations  were  then  added  to  the  normalized  free  Cortisol  levels 
to  yield  the  total  effective  free  corticosteroid  concentration  at  any  given  time  points 
during  the  24-hour  dosing  interval.  These  plots  are  shown  in  Figure  B-l  for  4  AM 
administration  and  Figure  B-2  for  4  PM  administration. 

Two  important  facts  were  drawn  from  these  simulations.  The  first  is  that 
budesonide  administration  at  4  AM  prevents  Cortisol  from  increasing  as  it  does  under 
baseline  conditions  (after  placebo).  Since  lymphocyte  counts  decline  as  Cortisol 
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concentrations  rise,  the  administration  of  budesonide  indirectly  blocked  the  decline  in 
lymphocytes  that  is  due  to  the  presence  of  Cortisol. 


Free  Cp(Cort,BUD)  after  0.1mg  BUD(IV)  -  4am 


CP(EFFECTIVE-TOTAL)  vs  T 
BUD  vs  T 

CORTISOL(EQUIV)  vs  T 
CORTISOL(BASELINE)  vs  T 


Figure  B-l .  Effective  corticosteroid  concentrations  after  administration  at  4  AM. 


Free  Cp(Cort.BUD)  after  0.1  mg  BUD(IV)  -  4pm 


CP(EFFECTIVE-TOTAL)  vs  T 
BUD  vs  T 

CORTISOL(EQUIV)  vs  T 
CORTISOL(BASELINE)  vs  T 


Figure  B-2.  Effective  corticosteroid  concentrations  after  administration  at  4  PM. 


The  second  is  that  at  upon  administration  at  4  AM,  Cortisol  levels  return  to  basline  levels 
only  during  the  acrophase.  Thus,  during  the  24-hour  dosing  interval,  Cortisol  can  not 
induce  a  decline  in  lymphocytes.  The  changes  in  Cortisol  levels  are  shown  in  Figure  B-3 
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and  Figure  B-4  in  which  Cortisol  levels  were  not  converted  to  budesonide  equivalents. 


Cortisol  and  BUD  after  1  mg  IV  -  4am 


CORTISOL(BASELINE)  vs  T 
CORTISOL  vs  T 
BUD  vs  T 


Figure  B-3.  Corticosteroid  concentrations  after  administration  at  4  AM. 


Cortisol  and  BUD  after  1  mg  IV  -  4pm 


Time  (hr) 


CORTISOL(BASELINE)  vs  T 
CORTISOL  vs  T 
BUD  vs  T 


Figure  B-4.  Corticosteroid  concentrations  after  administration  at  4  PM. 


These  simulations  illustrate  that  Cortisol  suppression  is  the  basis  for  the  negligible 
effects  on  lymphocytes  upon  the  administration  of  certain  doses  of  corticosteroids.  No  net 
effects  are  observed  when  the  free  concentration  of  the  exogenous  corticosteroid  is 
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proportional  to  the  decline  in  free  Cortisol  equivalents.  Although  Cortisol  levels  are 
suppressed,  they  are,  in  effect,  replaced  by  the  exogenous  corticosteroid.  The  result  is  an 
effective  corticosteroid  concentration-time  profile  that  is  quite  similar  to  the  baseline 
circadian  rhythm  of  Cortisol  and  baseline  lymphocyte  counts  are  sustained. 


APPENDIX  C 

PHARMACOKINETIC/PHARMACODYNAMIC  CONSIDERATIONS  IN 

DESIGNING  A  DOSING  REGIMEN 


Most  often,  patient-specific  dosing  recommendations  are  implemented  using 
clinically  determined  pharmacokinetic  (PK)  parameters.  However,  since 
pharmacodynamic  (PD)  variability  within  a  population  normally  exceeds 
pharmacokinetic  variability,  a  major  source  of  interpatient  diversity  is  often  omitted, 
namely  the  dynamic  parameters.  As  illustrated  below,  PK/PD  models  may  be  used  to 
design  a  dosing  regimen.  In  addition  to  the  common  PK  parameters,  the  PK/PD  approach 
requires  knowledge  of  EC50,  the  drug  concentration  that  produces  50%  of  a  given  effect. 
With  this  information,  however,  both  the  appropriate  dose  and  dosing  interval  may  be 
easily  determined  to  maintain  the  desired  effect  for  a  given  patient  (Stark  and  Hochhaus 
1998).  Due  to  the  complexity  of  the  PK/PD  model  used  to  describe  the  decline  in 
lymphocytes  after  corticosteroid  administration,  this  approach  could  not  be  readily 
applied.  Simpler  models  were  investigated  in  the  development  of  this  concept  and  the 
PK/PD  approach  to  dosing  recommendations  is  presented  here  using  the  receptor- 
mediated  response  model. 

Before  the  appropriate  dose  may  be  determined,  it  is  necessary  to  state  the  desired 
effect  as  fractions  of  the  maximum  possible  effect  (Emax).  In  the  expressions  below,  a 
refers  to  the  peak  effect  just  after  IV  bolus  administration  and  p  refers  to  the  trough 
effect.  To  obtain  a  desired  effect  within  a  given  effect  window,  a*Emax  decreasing  to 
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p*Emax  before  the  next  dose,  the  dosing  interval  (t)  for  multiple  IV  bolus  doses  (D)  and 
may  be  calculated  as 

"«•(!-/?)' 


t  =  —  In 
k. 


Eq.  C-l 


where  ke  is  the  first-order  elimination  rate  constant  of  the  drug.  The  recommended  dose  is 
then 

D  = 2— ^ '-  Eq.  C-2 

(1-or) 

where  Vd  is  the  volume  of  distribution  for  the  drug. 

Although  these  expressions  are  useful  in  designing  a  dosing  regimen  to  maintain 

effects  within  a  given  effect  window,  the  possibility  of  undesired  effects  must  also  be 

considered.  Since  theophylline  is  one  common  drug  having  a  rather  narrow  therapeutic 

window,  it  is  used  here  as  an  example.  Population  average  PK  and  PD  parameters  used  in 

this  example  are  provided  in  Table  C-l  (Holford,  Hashimoto  et  al.  1993). 

Table  C-l .  Pharmacokinetic  and  pharmacodynamic  properties  of  theophylline. 

Parameter  Value 


Emax  (%)*  63 

EC50(mg/L)  10 

ke  (h"1)  0.08 

Vd  (L)  30 

*  Emax  is  the  maximum  effect  of  theophylline  in  asthma  patients  expressed  as  a 
percentage  of  normal  forced  expiratory  volume  (FEVi). 


The  published  parameters  were  used  to  calculate  dosing  recommendations  that 
provide  a  wide  range  of  effects  (percentages  of  the  maximum  forced  expiratory  volume). 
A  summary  of  these  calculations  is  provided  in  Table  C-2.  In  viewing  the  results,  the 
normal  therapeutic  concentration  range  for  theophylline  provides  an  effect  between  50- 
67%  of  the  maximum  effect.  Using  the  data  in  Table  C-l,  this  maximum  effect  may  be 
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translated  into  a  32-42%  improvement  in  lung  function  for  an  asthmatic  patient  on  a 

dosing  regimen  of  320  mg  theophylline  every  eight  hours. 

Table  C-2.  Dosing  recommendations  for  theophylline  based  on  PK/PD  methods- 
Effect  Dose        Interval      Peak  Cp^*      TrpughCpss 
(%ofEmax)         (mg)  (h)  (mg/L)  (mg/L)  ~ 
99->l%            30,000          84               990  1.1 
80^20%            1100            34                 40                    2.5 


67->50%  320  8  20  10 


50->30%  170  11  10  4.3 

*Cpss  represents  the  plasma  concentration  of  theophylline  for  multiple  IV  bolus  doses  at 
steady- state. 


Other  calculations  are  provided  to  illustrate  the  limitation  of  the  PK/PD  approach 
when  narrow  therapeutic  window  drugs  are  considered.  In  such  cases,  it  is  not  possible  to 
arbitrarily  choose  a  desired  effect  range.  To  be  of  practical  use  for  such  drugs,  it  is 
necessary  to  consider  the  plasma  concentration-side  effect  relationships.  This  has  not  yet 
been  incorporated  into  the  theoretical  framework  presented  here. 

For  drugs  having  wide  therapeutic  windows,  it  is  possible  to  administer  a  larger 
dose  and  increase  the  dosing  interval  in  order  to  achieve  better  patient  compliance.  For 
maintaining  the  effect  above  some  threshold  value,  P*Emax,  the  appropriate  dosing 
interval  for  a  predetermined  dose  is 

D-(\-p) 


t  =  —  In 
k.. 


+  1 


Eq.  C-3 


J3-Vd-ECso 

This  expression  may  be  rearranged  in  order  to  calculate  an  appropriate  dose  once  a 
convenient  dosing  interval  has  been  chosen, 

fi-ECx.Vd.(e'*--l) 
D~  (!_£)  Eq.C-4 

During  a  multiple  dosing  regimen,  plasma  concentrations  increase  as  drug 
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accumulates  in  the  body.  For  first-order  elimination,  the  peak  plasma  level  at  steady-state 
is  augmented  by  the  same  factor  as  the  trough  level  when  compared  to  the  peak  and 
trough  after  the  first  dose,  respectively.  Effects  do  not  accumulate  and  are  dependent 
solely  on  the  drug  concentration  at  a  given  time.  However,  it  is  possible  to  define 
accumulation  factors  for  the  peak  and  trough  effects  independently.  These  accumulation 
factors,  rPeak  and  rtr0ugh,  are  useful  for  predicting  the  fluctuation  in  the  observed  effect  for 
a  given  dosing  regimen  of  a  drug  with  known  pharmacokinetic  parameters.  These 
accumulation  factors  are 

I  Cp0  ) 

rpeak    ~  {  r„      \  E<T  C-5 


§r)(l-e"'")+1 


and 


_  I CPo  J P    r, 

'trough  /  \  ■  tq.  L.-0 

^   -(l-e-^J-e-^+l 

\CPo  J 

where  Cp0  is  the  peak  plasma  concentration  observed  after  a  single  IV  bolus  dose. 

In  the  dosing  regimen  expressions  and  accumulation  factors  described  above,  the 
EC50  value  is  required.  For  a  receptor-mediated  pharmacological  effect,  the  EC50  value 
may  be  determined  by  fitting  data  to  a  PK/PD  model.  This  requires  several  blood  samples 
to  be  collected  over  a  given  time  period.  Since  this  is  impractical  outside  of  a  clinical  trial 
or  hospital  setting,  an  alternative  method  is  presented  here.  By  obtaining  a  single 
concentration/effect-time  point  after  two  different  doses,  the  patient-specific  EC50  value 
may  be  calculated.  To  prove  this,  let  Ea(t)  and  Ep(t)  be  the  measured  effects  at  time  t  after 
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drug  administration  of  doses  a  and  P,  respectively, 
and 

where  Ca(t)  and  Cp(t)  are  the  plasma  concentrations  at  time  t  after  doses  a  and  p.  The 
times  need  not  be  the  same.  Solving  Eq.  C-7  and  Eq.  C-8  for  Emax  and  setting  them  equal 
to  each  other  gives 

_  Ea  (Q  •  (EC50  +  Ca  (Q)     £,  (Q  ■  (£C50  +  C,  (Q) 

c.(0  c,<r)  ■  bqC"9 

These  equations  may  be  solved  for  EC50  to  yield 

"■"~£.<*C,(i)-£,(0-C.(0-  Eq'C"10 

Substituting  this  expression  for  EC50  into  Eq.  C-7  also  allows  Emax  to  be  determined  from 
the  same  two  data  points, 


E      -E"{t) 


max 


c.<0-c#W.M-^(0)+ 


£a  (0-^(0-^(0-^(0 

Despite  the  limitations  of  some  of  the  PK/PD  concepts  presented  here,  they  are 
potentially  useful  for  optimizing  drug  therapy  for  a  given  patient  and  for  understanding 
the  fluctuations  in  effects  after  drug  administration. 


APPENDIX  D 

NOVEL  NONCOMPARTMENTAL  PARAMETERS  IN 

PHARMACOKINETIC/PHARMACODYNAMIC  ANALYSIS 


Introduction 

Noncompartmental  methods  have  been  used  in  various  respects  throughout  the 
present  pharmacokinetic/pharmacodynamic  (PK/PD)  study.  However,  in  the  analysis  of 
kinetic/dynamic  data,  the  term  "noncompartmental"  is,  perhaps,  a  misnomer  and  a  better 
term  may  be  "model-independent"  analysis.  In  pharmacokinetic  data  analysis, 
noncompartmental  methods  refer  to  the  assessment  of  pharmacokinetic  parameters  such 
as  clearance,  bioavailability,  and  the  mean  residence  time  without  fitting  the  data  to 
conventional  compartment  models;  Calculations  are  based  on  the  area  under  the 
concentration-time  profile  of  the  drug.  For  application  to  PK/PD  studies,  calculations  are 
based  on  the  area  under  the  effect-time  profile  without  the  assumption  of  a  specific 
kinetic-dynamic  model. 

Noncompartmental  methods  provide  a  general  approach  for  pharmacokinetic 
analysis  without  requiring  the  assumption  of  a  specific  compartmental  model  for  the 
drug.  There  has  been  a  distinct  shift  away  from  computer-based  curve-fitting  of 
experimental  data  and  elaboration  of  compartmental  models  and  toward 
noncompartmental  methods  of  analysis  (Gibaldi  and  Perrier  1982).  In  some  cases,  the 
choice  of  a  specific  compartmental  model  is  somewhat  subject.  For  example,  a 
distribution  phase  is  required  to  adequately  describe  the  concentration-time  profile  after 
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the  IV  bolus  administration  of  most  drugs  and,  thus,  a  two-compartment  model  is  used. 
However,  after  oral  administration,  the  distribution  phase  is  not  always  apparent  and, 
thus,  a  one-compartment  model  with  first-order  absorption  may  adequately  describe  the 
concentration-time  profile  despite  the  fact  that  distribution  still  takes  time. 

Similarly,  the  use  of  noncompartmental  parameters  in  PK/PD  analysis  allows 
effects  to  be  quantified  in  a  rather  straightforward  and  unbiased  manner.  This  approach  is 
advantageous  due  to  the  complexity  of  the  PK/PD  models  necessary  to  describe  some 
effects  and,  thus,  the  rationale  for  establishing  noncompartmental  methods  in  PK/PD 
analysis  is  comparable  to  that  for  defining  noncompartmental  parameters  in 
pharmacokinetics.  In  the  present  study,  cumulative  Cortisol  suppression  (CCS)  and 
cumulative  lymphocyte  reduction  (CLR)  were  used  to  quantify  the  effect  of  budesonide 
on  Cortisol  and  lymphocytes,  respectively.  These  parameters  were  based  solely  on  the 
area  under  the  cortiol-time  or  lymphocyte  count-time  profile  and,  although  PK/PD 
modeling  was  also  performed  using  the  same  data  sets,  no  curve-fitting  was  required  to 
generate  them. 

Concepts 

The  time  course  of  drug  concentration  may  be  regarded  as  a  statistical  distribution 
curve.  For  this  reason,  noncompartmental  analysis  is  sometimes  referred  to  as  statistical 
moment  analysis.  The  first  application  of  statistical  moment  analysis  to  biological 
systems  was  to  describe  the  input-output  of  cholesterol  in  the  body  (Perl  and  Samuel 
1969).  These  concepts  were  later  applied  to  pharmacokinetics  and  data  analysis  may  be 
accomplished  by  the  same  approach  regardless  of  the  complexity  of  the  model  that  may 
best  describe  the  data  (Yamaoka,  Nakagawa  et  al.  1978).  One  of  the  basic  parameters  in 
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noncompartmental  analysis  is  the  mean  residence  time  (MRT),  which  refers  to  the  mean 
time  for  the  drug  to  remain  intact  in  the  body  after  administration.  As  will  be  discussed, 
MRT  is  defined  in  terms  of  the  area  under  the  curve  (AUC)  of  the  concentration-time 
profile.  Although  the  use  of  AUC  is  well-established  in  pharmacokinetics,  applications  in 
PK7PD  modeling  have  been  limited.  The  idea  of  mean  effect  time  (MET)  has  been 
introduced  to  describe  the  average  duration  of  receptor  occupancy  for  the  triamcinlone 
acetonide-glucocorticoid  receptor  interaction  in  in-vivo  animal  experiments  (Gonzalez- 
Rothi,  Suarez  et  al.  1996).  Although  no  formal  derivation  of  MET  was  presented,  this 
paper  contained  a  straightforward  application  of  noncompartmental  analysis  for  the 
effect-time  profile. 

Mean  residence  time  and  the  PK/PD  equivalent  mean  effect  time  are  based  on 
statistical  moments  about  the  time-axis  (or  x-axis  on  the  x,y-plane).  This  is  only  one  of 
two  possible  moments  as  a  statistical  moment  may  also  be  defined  about  the  y-axis 
(concentration  or  effect  in  the  present  discussion).  Thus,  another  noncompartmental 
parameter  may  be  defined,  which  is  the  mean  concentration  or  mean  effect.  This 
parameter  is  a  weighted  average  of  all  y-values  in  a  given  interval  on  the  x-axis.  The  two 
moments  about  the  x  and  y  axes  are  shown  in  Figure  D-l  in  a  very  general  format  and  the 
concepts  apply  whether  or  not  baseline  values  are  present.  For  most  pharmacokinetic 
studies,  an  exogenous  drug  is  administered  and  there  is  no  baseline  (placebo)  level.  Thus, 
noncompartmental  parameters  in  pharmacokinetics  are  defined  using  the  area  between 
the  concentration-time  profile  and  the  x-axis.  This  is  not  always  the  case  for  PK/PD 
studies,  depending  on  the  measured  effect.  In  the  present  PK/PD  study,  effects  of 
corticosteroids  were  monitored  using  Cortisol  and  lymphocytes,  both  of  which  have 
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baseline  values.  Thus,  noncompartmental  parameters  in  PK/PD  are  defined  using  the  area 
between  the  observed  effect-time  profile  after  placebo  and  that  after  drug  administration. 
Together,  the  statistical  moments  about  both  axes  define  the  center  (center  of  gravity  in 
mathematics  and  physics)  of  the  area  under  the  data-time  profile  or  the  area  between  the 
time  profiles  after  treatment  and  placebo. 


R  =  The  region  bound  by  f(x)  and  g(x)  on  [a,b] 


H 


"►*(Mx,My) 


a 


f(x)  >  g(x)  for  a  <  x  <  b 


Figure  D-l .  Statistical  moments  about  the  x  and  y  axes. 
In  pharmacokinetics,  My  is  the  statistical  moment  about  the  y-axis  used  in  calculating  the 
mean  residence  time  or  mean  effect.  Mx  is  the  statistical  moment  about  the  x-axis  used  to 
calculate  the  weighted  average  concentration  or  effect  over  the  time  interval  [a,b].  In 
addition  to  the  two  moments  Mx  and  My,  the  area  of  region  R  is  needed  for  the 
calculations.  This  area  is  often  referred  to  as  the  zero  moment.  Together,  the  two 
moments  define  the  center  point  (center  of  gravity)  of  the  area  bound  by  f(x)  and  g(x). 


The  concept  of  using  two  statistical  moments  to  characterize  data  is  presented 
below  in  simple  mathematical  terms  for  functions  on  the  x,y-plane. 
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Let  f(x)  and  g(x)  be  continuous  functions  on  the  interval  [a,b]  with 

f(x)  >  g(x)  for  a  <  x  <  b. 
Let  R  be  the  region  between  f(x)  and  g(x)  on  [a,b]. 
The  moment  Mx  of  R  about  the  x-axis  is 

Mx  =  \\-{f{x)]2-[g{x)]2\dx.  Eq.D-1 

a 

If  g(x)  =  0,  this  reduces  to 

Mx  =  \l-\f(x)]2.dx.  Eq.D-2 

The  moment  My  of  R  about  the  y-axis  is 

b 

My  =  \x  ■  lf(x)  ~  g(*)]  dx.  Eq.  D-3 

a 

If  g(x)  =  0,  this  reduces  to 

b 

My  =  jxf(x)dx.  Eq.  D-4 

a 

These  moments  may  be  used  to  define  the  center  of  gravity,  the  midpoint  of  the  area 
bound  by  f(x)  and  g(x)  on  [a,b],  of  the  region  R.  If  R  has  a  positive  area  (A),  the  center  of 
gravity  of  R  is  the  point  ( x ,  y )  where 

-     My         -     Mx 

x  =  — -  and  y  =  — x-  Eq.  D-5 

A  A  n 

and  the  area  A  is  given  by  the  integral 

b 

A  =  JVO)  -  g(x)}  ■  dx  Eq.  D-6 

a 

or 
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u 

A=jf(x)-dx  ifg(x)  =  0. 

a 

In  pharmacokinetics,  mean  residence  time  (MRT),  the  value  of  x  in  Eq.  D-5,  is 
calculated  using  My  and  A  from  Eq.  D-4  and  Eq.  D-6,  respectively.  In  these  equations, 
the  function  f(x)  is  Cp(t),  the  drug  concentration  in  plasma  at  a  given  time  t,  g(x)  is  equal 
to  zero  (assuming  administration  of  an  exogenous  drug),  and  all  integrals  are  over  [0,oo]. 
The  mean  residence  time  is  then 


\t-Cp{t)-dt 

\AT>T       o  AUMC 

MRT  =  ^ = pa  n_7 

co  AUC  bqU/ 

\CP-dt 


where  AUMC  in  the  area  under  the  first  moment  curve  and  AUC  is  simply  the  area  under 
the  curve  (the  zero  moment  curve).  Although  Gibaldi  refers  to  MRT  as  the  first  moment 
of  the  concentration-time  curve  (Gibaldi  and  Perrier  1982),  it  is,  more  precisely,  the  first 
normal  moment,  being  normalized  by  dividing  through  by  the  zero  moment. 

Although  mean  residence  time  may  be  stated  by  definition  (based  on  elementary 
theorems  in  calculus,  it  may  also  be  approached  conceptually.  A  derivation  of  mean 
residence  time  is  presented  below  to  illustrate  that  the  mathematical  definition  given  in 
Eq.  D-7  describes  the  mean  time  that  a  drug  molecule  resides  (intact)  in  the  body  after 
administration.  The  derivation  is  loosely  based  on  that  given  by  Rowland  and  Tozer 
(Rowland  and  Tozer  1995).  The  means  residence  time  in  pharmacokinetics  is  equal  (or 
proportional)  to  the  mean  effect  time  when  drug  action  may  be  described  by  a  linear 
effect  PK/PD  model  and,  thus,  an  understanding  of  the  pharmacokinetic  principle  may 
lay  the  foundation  for  the  more  complex  application  to  PK/PD. 
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Mean  Residence  Time 

Each  dose  of  a  drug  contains  millions  of  molecules.  As  an  illustration,  consider 
the  9  mg  dose  of  budesonide  administered  in  the  clinical  studies  in  Chapter  4.  This  dose 
corresponds  to  roughly  2.4  x  1030  molecules.  These  molecules  spend  different  amounts  of 
time  in  the  body.  Some  molecules  are  eliminated  quickly  and  some  distribute  into  various 
tissues.  Obviously,  not  all  molecules  will  interact  with  corticosteroid  receptors  to  initiate 
an  effect.  The  drug  course  through  the  body  is,  thus,  a  stochastic  process.  The  term 
stochastic  is  derived  from  the  Greek  word  related  to  a  target.  When  arrows  are  aimed  at  a 
target,  most  will  not  hit  the  exact  center  and  there  is  a  random  scatter  pattern.  It  is  easy  to 
understand  why  this  word  has  become  part  of  the  jargon  in  statistical  analysis. 

Applying  the  concept  of  stochastic  processes  allows  the  concentration-time 
profile  for  a  drug  to  be  treated  as  a  statistical  distribution  curve.  It  is  not  possible  to 
monitor  the  fate  of  each  individual  molecule  in  the  body  and  only  the  overall  properties 
of  large  collections  of  molecules  (that  is,  a  sample  or  population)  can  be  observed.  The 
residence  time  of  the  drug  in  the  body  can,  therefore,  be  regarded  as  a  statistically 
random  variable.  Since  time  points  are  considered  quantitative  data,  it  is  possible  to 
determine  the  mean  (u),  which  is  the  arithmetic  average  of  the  measurements  in  a  data 
set.  It  is  equally  possible  to  determine  the  variance  (ct2)  of  residence  time,  VRT 
(Yamaoka,  Nakagawa  et  al.  1978;  Gibaldi  and  Perrier  1982).  This  involves  the  second 
central  moment  of  the  concentration-time  profile  (in  contrast  to  the  first  moments  defined 
in  Eq.  D-l  through  Eq.  D-4)  and  will  not  be  covered  here. 

Consider  the  IV  bolus  administration  of  a  dose  containing  N  molecules  of  a  drug 
molecule.  The  mean  residence  time  of  all  drug  molecules  in  the  body  can  be  written  as 
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MRT  =  t^  Eq.D-8 

where  tj  is  the  residence  time  of  the  ith  molecule.  As  stated  previously,  individual 
molecules  can  not  be  counted.  However,  it  is  possible,  theoretically,  to  group  together 
molecules  having  the  same  residence  time.  Let  m  be  the  number  of  molecules  having 
residence  time  ti.  The  total  amount  of  residence  time  for  this  group  is  m.ti.  The  average 
residence  time  for  this  group  (MRTO  is,  trivially, 

MRTX  =  5L4-  =  t  .  E    D.9 

"1 

Consider  another  group  of  n2  molecules  that  remain  intact  in  the  body  for  time  t2.  The 

average  residence  time  of  both  groups  (m  and  n2)  of  molecules  together  is  given  by 

MRTi2  =  -!— J L_2. .  E     D.10 

This  procedure  may  be  continued  such  that  each  of  the  N  molecules  in  the  dose  can  be 
placed  into  m  different  groups  representing  every  possible  residence  time  observed.  The 
mean  residence  time  of  all  m  groups  (m,  n2,  . . .,  nm)  together  is  then 

MRTu,.,m  =  ^J — —^  ■  Eq.  D-l  1 


nx  +n2  +•••  +  «. 


Thus, 


"'  in 


*     i 


i 

The  summation  in  Eq.  D-l  2  may  be  extended  to  include  all  values  of  time  t,  multiplying 
each  time  tj  by  the  number  of  occurrences  nj.  Thus, 
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MRT  =  -^-—,  Eq.D-13 

N 

which  may  be  rearranged  to  give 

MRT  =  Z[^\t-  Eq.D-14 


The  bracketed  quantity 


v7Vy 


is  the  probability  of  observing  the  value  t.  Since  N  is  very 


large,  the  probability  of  observing  any  given  t  is  infinitesimal.  The  summation  over 
infinitesimal  probabilities  is  equivalent  to  integration.  Therefore, 

oo 

\tdn 

MRT  =  \—.  Eq.D-15 

N 

The  transition  from  Eq.  D-14  to  Eq.  D-15  may  also  be  accomplished  by  letting  Eq.  D-14 
be  regarded  as  a  Riemann  sum  and  using  the  approximation 

n  b 

£/('*)-Az*  *//(*)•<&•  Eq.D-16 

k=i  a 

In  Eq.  D-15,  dn  is  the  number  of  molecules  from  time-zero  to  infinity,  the  time  at 
which  all  molecules  have  been  eliminated.  The  derivation  now  shifts  from  the  molecular 
to  the  macroscopic  level  by  converting  numbers  of  molecules  into  amounts  (g)  using 
Avagadro's  number  and  the  molecular  weight  (MW)  of  the  drug.  Let  dX  be  the  amount 
of  drug  eliminated, 

,v       dn  ■  MW 

^'J^^  E*D-17 

and  X00  be  the  amount  of  drug  eliminated  at  time-infinity  (the  entire  dose  consisting  of  N 
molecules).  The  expression  for  MRT  in  Eq.  D-15  may  now  be  written  as 
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jt-dX 


MRT  =  JLx^~-  EqD-18 

By  rewriting  the  expression  in  terms  of  the  amount  of  drug  eliminated,  MRT  may  be 
defined  using  common  concepts  and  parameters  in  pharmacokinetics. 

The  rate  of  elimination  (dX/dt)  of  a  drug  may  be  expressed  in  terms  of  total  body 
clearance  CL  (assumed  to  be  constant), 

—  =  CL-Cp(t)  Eq.D-19 

where  Cp(t)  is  the  drug  concentration  in  plasma  at  time  t.  Solving  Eq.  D-19  for  dX, 

dX  =  CL-Cp(t)-dt,  Eq.D-20 

and  substituted  into  Eq.  D-18  to  give 

oo 

\t-CL-Cp(t)-dt 
MRT  =  * — .  Eq.  D-2 1 

At  this  point  in  the  derivation,  it  is  tempting  to  replace  X00  by  the  dose  (D).  However, 
making  this  substitution  leaves  CL  in  the  numerator.  This  factor  does  not  have  to  remain 
if  the  following  considerations  are  made. 
Total  body  clearance  (CL)  may  be  expressed  as 

y*X) 

CL  =  Z •  Eq.  D-22 

\Cp{t)-dt 

0 

Solving  Eq.  D-22  for  X00  gives 

CO 

X~=CL-\Cp(t).dt  =  CL-AUC  Eq.D-23 
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where  AUC  is  the  area  under  the  curve  (concentration-time  curve).  Substituting  this 
expression  into  Eq.  D-21  yields 

CO  CO 

\t-CLCp[t)dt      \t-Cp(t)-dt 

LL  ■  A  UC  A  UC 

Recognizing  the  numerator  as  the  area  under  the  first  moment  curve  (AUMC)  allows  Eq. 
D-24  to  be  rewritten  as 
.,__,     AUMC 

MRT  =  ^Uc->  E*D-25 

which  is  equal  to  Eq.  D-7  for  MRT  stated  by  definition  from  basic  concepts  in  calculus. 

Q.E.D. 
Statistical  Moments  in  Pharmacokinetics 
Mean  Residence  Time 

As  shown,  mean  residence  time  (MRT)  is  the  average  time  that  a  drug  remains 
intact  in  the  body  after  administration.  After  IV  bolus  administration  of  a  drug  following 
a  one-compartment  model,  MRT  is  inversely  proportional  to  the  first  order  elimination 
rate  constant  (ke).  To  prove  this,  consider  the  statistical  moment  about  the  y-axis  or 
concentration-axis  of  the  concentration-time  curve  (Mc  )  as  in  Eq.  D-3, 

MCp=\t-Cp(t)-dt  =  AUMC  =  ^-]t-e-k<'-dt  Eq.  D-26 

o  "d    o 

where  AUMC  is  the  area  under  the  first  moment  curve  (about  the  y-axis),  D  is  the  IV 

bolus  dose,  and  Vd  is  the  volume  of  distribution.  This  integral  may  be  solved  using 

jx"  •e-qx  -dx  =  -^      for      n>\,      q>0.  Eq.  D-27 

o  9 
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Then, 


D    V.         D        C.(0) 

Uc'-y^-TiTd-^r-AUMC 


The  area  under  the  curve  is 


D    __,,    ^         D         C,(0) 


AUC=lCp(t).dt=j^.dt  =  -^  =  ^  Eq.D-29 

0  0  Y  d  Kt'V  d  Ke 

Thus,  the  mean  residence  time  is 

..__,     AUMC     C,(0)      *e  1 

MRT  = =  — e- —  =  —  pn    n  T.0 

AUC  k)      Cp(0)     *e  bq'D"3° 

Another  potentially  useful  relationship  using  MRT  is  that  for  the  volume  of 
distribution  at  steady-state  (Vdss), 

Vdss=CL-MRTIV,  Eq.D-31 

where  CL  in  the  total  body  clearance  of  the  drug,  related  to  the  area  under  the  curve 
(AUC)  as  in  Eq.  D-23  by 

n         D 

CI  =  — •  Eq.D-32 

Weighted  Average  Concentration 

The  moment  about  the  x-axis  may  be  used  to  quantify  the  weighted  average 
concentration  achieved  by  a  certain  dose.  Recall  that  for  a  single  dose,  it  is  not  possible  to 
calculate  the  average  concentration  ( Cp  )  using  the  AUC  and  the  interval  (t), 

jr     AUC 

Cp  = •  Eq.  D-33 

T 

Since  the  time  interval  is  [0,  oo],  the  average  concentration  after  a  single  IV  bolus  dose 
according  to  Eq.  D-33  is  zero, 
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I 


c,-A 


— >  0       <zs       /  — >  oo , 


This  is  not  the  case  when  calculating  a  weighted  average  based  on  statistical  moments 
(C)  according  to  Eq.  D-2.  The  moment  about  the  time-axis  is 


M,       -fjMJ.A-I.^USC-I.jf^ 


-|2 


-*.-' 


oL^ 


£# 


D2 


*-yi-K 


■e 


.-2k -i 


D2 


Cl(P) 


o    4-f;.*.     4-*. 


Eq.  D-34 


The  weighted  average  concentration  is  then 

1 


C  = 


4£/C  4 


Eq.  D-35 


where  AUSC  is  the  area  under  the  squared  curve  (concentration-time  profile-squared). 
Thus,  the  weighted  average  concentration  after  an  IV  bolus  dose  is  one-forth  of  the  initial 
concentration  at  time-zero.  Unlike  MRT,  the  weighted  average  concentration  is, 

necessarily,  dose-dependent.  The  point  (MRT,  C )  defines  the  center  of  the  concentration- 
time  profile. 

Statistical  Moments  in  Pharmacokinetics/Pharmacodynamics 
This  discussion  of  statistical  moments  applied  to  PK/PD  data  will  be  limited  to 
the  linear  effect  model  and  the  sigmoid  Emax  model  for  a  receptor-mediated 
pharmacological  response.  As  in  the  preceding  application  to  pharamacokinetics,  all 
derivations  will  be  for  a  single  IV  bolus  dose.  From  Eq.  D-l  through  Eq.  D-5,  the 
statistical  moments  about  the  time  and  effect  axes  may  be  used  to  calculate  the  mean 
effect  time  (MET)  and  the  weighted  average  effect  (£),  respectively,  over  a  given  time 
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interval.  The  general  expressions  for  these  two  noncompartmental  PK/PD  parameters  are 
AUMCE 


MET  = 


AUCE 


Eq.  D-36 


where  AUMCE  is  the  area  under  the  first  moment  of  the  effect  curve, 

00 

AUMCE  =\t-E{t)dt,  Eq.D-37 


AUCE  is  the  area  under  the  effect  curve, 


AUCE=  JE(t)dt, 


Eq.  D-38 


and 


1 


AUSCE 


E  = 


AUCE 


Eq.  D-39 


where  AUSCE  is  the  area  under  the  squared  curve  of  the  effect, 


AUSCE  =  JE\t)dt 


Eq.  D-40 


These  two  parameters  define  the  center  of  gravity  of  the  area  bound  by  the  effect-time 
profile,  (MET,  E ).  MET  may  be  considered  a  very  useful  parameter  since  the  average 
duration  of  a  pharmacological  effect  may  enable  suitable  dosing  intervals  to  be  predicted 
for  a  given  drug.  Another  method  for  determining  dosing  intervals  based  on  a  PK/PD 
approach  was  presented  in  Appendix  C.  The  E  may  not  seem  as  appealing  as  MET  since 
it  merely  defines  a  point  within  the  region  bound  by  the  effect-time  profile  rather  than  a 
true  observed  effect.  However,  E  may  be  used  to  compare  the  average  effects  observed 
after  the  administration  of  different  drugs  or  different  doses.  Recall  that  the  average  effect 
determined  in  a  more  conventional  manner  using  the  area  under  the  effect-time  profile  is 
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equal  to  zero  for  a  single  dose  over  the  interval  [0,oo], 

\E{t)-dt 

E  = >  0       as       t  — >  oo . 

t 

It  is  important  to  note  that,  unlike  noncompartmental  pharmacokinetic  parameters,  both 

noncompartmental  PK/PD  parameters  (MET,  E)  are  dose-dependent. 

The  lack  of  additional  applications  of  statistical  moments  to  PK/PD  data  is,  in 
part,  due  to  the  rather  complex  mathematical  expressions  needed  to  describe  some 
effects.  This  was  especially  true  of  the  expressions  needed  for  Cortisol  and/or 
lymphocytes,  which  were  used  as  surrogate  markers  for  systemic  effects  in  the  present 
study.  In  order  to  develop  the  theory  of  noncompartmental  PK/PD  analysis,  less  complex 
models  were  investigated  in  an  attempt  to  define  and  utilize  novel  PK/PD  parameters 
based  on  the  area  under  the  effect-time  profile. 
The  Linear  Effect  Model 

For  situations  in  which  the  observed  effect  is  directly  proportional  to  drug 
concentration  at  a  given  time,  the  linear  effect  model  may  be  used.  This  is  the  simplest 
effect  model  and  the  effect  as  a  function  of  time  is 

E(t)  =  a-Cp(t)  Eq.D-41 

where  Cp(t)  is  the  drug  concentration  at  time  t  and  a  is  a  proportionality  constant.  The 

MET  is  then 

MET  =  a  ■  MRT  Eq  D_42 

and  from  Eq.  D-30 

._     AUMCE  AUMC      a 

MET  = =  a =  — .  Fa  D-43 

AUCE  AUC       k  cq.u<u 
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The  derivation  of  MET  for  a  linear  effect  model  is  analogous  to  that  for  MRT.  Likewise, 
from  Eq.  D-35,  the  weighted  average  effect  is 

_     l.AUSCE         \-AUSC     a.C(0) 

E  =  A =  «•- = ^— .  Eq.  D-44 

AUCE  AUC  4 

Although  useful  for  developing  the  idea  of  statistical  moments  in  PK/PD  analysis,  this 
model  has  limited  application  since  few  drugs  exhibit  a  direct  linear  relationship  between 
concentration  and  pharmacological  effect  over  an  extended  concentration  range. 
The  Emax  Effect  Model 

A  potentially  more  useful  effect  model  is  the  so-called  Emax  model  for  a  receptor- 
mediated  response, 

£max-Cp(0  Eq.D-45 

EC50+Cp(t) 

where  Emax  is  the  maximal  effect,  Cp(t)  is  the  plasma  concentration  of  drug  at  a  given 
time,  and  EC50  is  the  concentration  at  which  50%  of  the  maximal  effect  is  observed. 
Since  the  effects  of  corticosteroids  are  initiated  by  interaction  with  the  glucocorticoid 
receptor,  the  majority  of  work  in  the  present  theoretical  study  focused  on  developing 
relevant  PK/PD  parameters  based  on  the  area  under  the  curve  for  an  effect  observed  upon 
the  formation  of  a  drug-receptor  complex.  As  shown  in  Figure  D-2,  in  order  to  calculate 
MET  and  E ,  three  areas  must  be  determined: 

(1)  The  area  under  the  effect  curve  (AUCE).  The  AUCE  may  be  determined  from  actual 
or  simulated  data  using  the  trapezoidal  rule  and  an  integrated  solution  for  AUCE  has  been 
provided  for  an  Emax  model  with  IV  bolus  input  and  first-order  elimination  (Derendorf, 
Hochhaus  et  al.  1993), 
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Effect  (E^  model) 


t*Effect-Time  profile 
after  IV  bolus  administration 


6  12  18 

Time  (hr) 


24 


Effect-squared  profile 
after  IV  bolus  administration 


6  12  18 

Time  (hr) 


24 


Figure  D-2.  Plotted  functions  for  statistical  moment  analysis  in  PK/PD. 
Plots  were  constructed  using  the  Emax  model  for  a  receptor-mediated  pharmacological 
effect  from  Eq.  D-45.  The  shaded  areas  in  the  plots  above  for  the  effect,  effect-times- 
time,  and  the  effect-squared  represent  the  integrals  for  AUCE,  AUMCE,  and  AUSCE 
given  in  Eq.  D-46,  Eq.  D-47,  and  Eq.  D-50,  respectively. 
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AUCE  =  ^^  -In 
k. 


i+.   D 


^50  '  *  d 


Eq.  D-46 


(2)  The  area  under  the  first  moment  curve  of  the  effect  (AUMCE).  From  Eq.  D-37, 
AUMCE  is 

AUMCE  =  )t-E{t)-dt  =  )t-  *™'C>®  -dt.  Eq.D-47 

0  0         ^50  +(~pV) 

However,  regardless  of  the  input  function  (for  concentration),  an  integrated  solution  of 
AUMCE  does  not  exist  in  closed  form.  A  very  simple  solution  for  this  integral  is  the 
Laplace  approximation  wherein  either  the  pre-exponential  factor  or  the  exponential 
function  is  neglected,  giving 

r-v^\  max   '^-p\') 

E(0=       ECP  Eq.D-48 


or 


£(0  *  ^max »  Eq.  D-49 

respectively.  Such  approximations  provide  good  results  only  for  small  effects 
proportional  to  the  concentration  that  disappear  at  a  rate  equal  to  the  elimination  rate  of 
the  drug  (Eq.  D-48)  or  for  very  large  effects  that  remain  near  the  value  of  Emax  for  a  given 
time  interval  (Eq.  D-49).  A  better  approximation  has  been  provided  based  on  theories  of 
asymptotic  functions  and  applied  to  yield  a  good  estimate  of  MET  (Krzyzanski  and  Jusko 
1997).  Another  solution  has  been  proposed  wherein  the  effect-time  profile  is  subdivided 
into  five  linear  or  exponential  regions  in  order  to  obtain  easily  integrated  functions 
(unpublished  algorithms),  providing  results  within  5%  of  those  obtained  by  numerical 
integration  (with  a  data  density  of  20  points  per  hour). 
(3)  The  area  under  the  squared  curve  of  the  effect  from  Eq.  D-40  is 
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AUSCE  = 


•In 


1  + 


Cp0 


EC 


>(.) 


Cp0 


EC50  +  Cp0 


Eq.  D-50 


where  Cp0  is  the  concentration  at  time-zero,  determined  by  the  dose  D  and  the  volume  of 
distribution  Vd.  Thus,  although  MET  can  not  be  calculated  exactly  such  is  not  the  case  for 
E  since  the  integrals  exist  in  closed  form.  This  parameter  is  not  dependent  on  clearance 
but  only  on  the  volume  of  distribution  (Vd),  the  IV  bolus  dose  (D),  Emax  and  the  EC50 
value, 


E  = 


1    E] 

■AUSCE      2'T 


2 
max 


In 


1  + 


Cp0 


Cp0 


EC50  +  Cp0 


AUCE 


In 


1  + 


D 


ECi0-Vd 


or,  in  reduced  form, 


In 


E  = 


1  + 


Cp0 


EC 


50 


Cp0 


EC50  +  Cp0 


2-ln 


1  + 


Cp0 


EC 


50 


Eq.  D-51 


Two  examples  of  the  application  of  statistical  moments  in  PK/PD  analysis  are  shown  in 
Figures  D-3  and  D-4  and  comparisons  are  made  to  quantifying  effects  and  drawing 
conclusions  only  from  the  area  under  the  effect-time  curve,  AUCE. 
A  Special  Case 

Although  no  general  solution  for  the  integral  of  the  effect-times-time  function, 
(AUMCE)  exists  in  closed  form,  such  is  not  the  case  for  one  interesting  condition.  If  an 
IV  bolus  dose  provides  an  initial  plasma  concentration  equal  to  EC50,  all  necessary 
integrals  reduce  to  ones  having  closed  form  solutions,  including  that  for  AUMCE. 
IfCp(0)=EC50, 
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Application:  Different  Pharmacokinetics 
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Figure  D-3.  Application  of  statistical  moments  in  PK/PD  analysis  for  drugs  with 
different  kinetic  paramters. 
D-IV  bolus  dose,  CL-total  body  clearance,  Vd-volume  of  distribution,  ke-elimination  rate 
constant,  Emax-maximum  response  possible,  EC50-concentration  that  yields  50%  of  the 
maximum  response,  Cp(0)-  the  initial  plasma  concentration  after  IV  bolus  injection, 
AUC-  area  under  the  concentration-time  profile,  AUMC-  area  under  the  concentration- 
times-time  profile,  E(0)-  initial  effect  after  IV  bolus  injection,  AUCE-  area  under  the 
effect-time  profile,  AUMCE-  area  under  the  effect-times-time  profile,  AUSCE-  area 
under  the  effect-squared-time  profile,  MRT-  mean  residence  time,  MET-  mean  effect 
time,  ME-  weighted  average  effect. 

Concentration-time  and  effect-time  profiles  were  simulated  for  two  drugs,  A  and  B, 
which  differ  in  volume  of  distribution,  using  a  one-compartment  body  model  with  IV 
bolus  input  and  an  Emax  effect  model,  Eq.  D-45.  To  produce  the  same  initial  concentration 
and  effect,  the  dose  of  B  (the  drug  with  the  larger  volume  of  distribution)  was  doubled. 
This  resulted  in  an  increase  in  both  AUCE  and  MET  by  a  factor  of  2  for  drug  B  compared 
to  drug  A.  However,  ME,  the  weighted  average  effect  over  the  same  time  interval  (zero  to 
infinity)  was  the  same  for  both  drugs.  Thus,  compensating  for  the  larger  volume  of 
distribution  of  drug  B  by  doubling  the  dose  produced  the  same  average  effect  but,  due  to 
the  pharmacokinetic  disposition  of  the  drug,  the  effect  was  observed  over  a  longer  period 
of  time. 
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Application:  Different  Kinetics  and  Dynamics 
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Figure  D-4.  Application  of  statistical  moments  in  PK/PD  analysis  for  drugs  with 
different  kinetic  and  dynamic  paramters. 
D-IV  bolus  dose,  CL-total  body  clearance,  Vd-volume  of  distribution,  ke-elimination  rate 
constant,  Emax-maximum  response  possible,  EC50-concentration  that  yields  50%  of  the 
maximum  response,  Cp(O)-  the  initial  plasma  concentration  after  IV  bolus  injection, 
AUC-  area  under  the  concentration-time  profile,  AUMC-  area  under  the  concentration- 
times-time  profile,  E(0>  initial  effect  after  IV  bolus  injection,  AUCE-  area  under  the 
effect-time  profile,  AUMCE-  area  under  the  effect-times-time  profile,  AUSCE-  area 
under  the  effect-squared-time  profile,  MRT-  mean  residence  time,  MET-  mean  effect 
time,  ME-  weighted  average  effect. 

Concentration-time  and  effect-time  profiles  were  simulated  for  two  drugs,  A  and  B, 
which  differ  in  volume  of  distribution  and  potency  (EC50  value),  using  a  one- 
compartment  body  model  with  IV  bolus  input  and  an  Emax  effect  model,  Eq.  D-45. 
Clearance  values  were  altered  slightly  to  facilitate  comparisons  between  the 
noncompartmental  parameters.  For  all  practical  purposes,  however,  clearance  is  roughly 
equal  for  both  drugs.  In  this  scenario,  the  areas  under  the  effect-time  profiles  were  equal 
and  MRT  and  MET  were  50%  greater  for  drug  B  compared  to  drug  A.  However,  the 
weighted  average  effect  was  50%  greater  for  drug  A  compared  to  drug  B.  This  example 
illustrates  the  usefulness  of  the  parameters  determined  using  statistical  moment  analysis. 
Although  the  doses  were  selected  to  achieve  the  same  initial  concentration,  the  AUC 
values  were  different  due  to  the  pharmacokinetic  disposition  of  the  drugs  yet  the  AUCE, 
the  areas  under  the  effect-time  profiles  were  identical.  Drug  A  produced  a  greater  effect 
with  a  shorter  duration. 
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m  =  __^_  £q  D_52 

AUCE  =  ^-^  =  Emm.txn,  Eq.D.53 


AUMCE  =  ^f£-,  Eq.D_54 


and 


12  -kl 


E2       f    r  i 

max  i      1^*  A 


^WSC£  =  -^.|ln[2]-i|,  Eq.D-55 

where  Emax  is  the  maximum  effect  of  the  drug  and  ke  is  the  first-order  elimination  rate 
constant.  Thus,  for  this  unique  dose,  MET  may  be  calculated  exactly  and  the  average 
duration  of  the  effect  is  directly  proportional  to  the  half-life  of  the  drug, 

AjfPT     AUMCE  n1  1nif% 

W"^ra"l2^ta2",,712'r,/I  Eq-D-56 

where  tm  is  the  half-life  of  the  drug. 

ME  may  be  expressed  as  a  fraction  of  the  maximal  effect, 

-AUSCE  c  1 

ME  =  ^ =  _!"£<_.  J  i L__l  =  oi39.F  Fn  n  S7 

AUCE  2      {      2-ln[2]J     VA^  *™-  Eq.  D-57 
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